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ABSTRACT OF THESIS submitted by Thomas Barford for the Degree of Doctor 
of Philosophy and entitled The RenormaUsation Group and AppUcadons in Few-Body 
Systems 

Month and Year of Submission: 2004 



This thesis is about effective field theories. We study the distorted wave renormali- 
sation group (DWRG), a tool for constructing power-counting schemes in systems where 
some diagrams must be summed to all orders. We solve the DWRG equation for a variety 
of long-range potentials including the Coulomb, Yukawa and inverse square. We derive es- 
tablished results in the case of the Yukawa and Coulomb potentials and new results in the 
case of the inverse square potential. We generalise the DWRG to systems of three bodies 
and use it to find the power-counting for the three-body force in the KSW EFT. This power- 
counting corresponds to perturbations in the renormalisation group about a limit-cycle. We 
also derive equations for the LO and NLO KSW EFT distorted waves for three bosons and 
for three nucleons in the i channel. The physical results in the nuclear system agree 
well, once the LO three-body force is determined, with predictions of more sophisticated 
potential models. 
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Chapter 1 



Introduction 

From an historical perspective, the hierarchy of scales to be found in the universe has con- 
trolled the development of physics since the middle ages. Time and again the theories of the 
past have been found to be only effective theories that approximate the physical world. 

The first great post-Renaissance physics theory, Newtonian mechanics, is still used 
extensively in many areas, yet it has been 'known' for a long time that it is only an approx- 
imation. The advent of relativistic mechanics brought fresh insight into the nature of the 
universe as well as explaining some mysterious behaviour in the orbit of the planet Mer- 
cury, yet its use in explaining the path of a cricket ball struck through mid-wicket has never 
been advocated. Relativistic mechanics is a far more complicated beast than its 'lesser' 
cousin, in many cases the pay-off for the extra work involved in a calculation is simply 
not worth it. The difference in the path of that cricket ball as explained by Newtonian and 
relativistic mechanics would be unmeasurable. 

Even though Newtonian mechanics is not fundamental to the nature of the universe as 
we know it now, it is still an effective theory; it explains what we see at speeds far below the 
speed of light, c. In quantifying this statement we might say that the difference between any 
classical or relativistic calculation, in which the typical speed is of order v, is of order v^/c^. 
We may then consider any problem, and knowing how accurately we want the result, chose 
between a full relativistic calculation, a simple classical calculation or maybe some middle 
ground in which the classical result is corrected to some order in v^/c^. 

This kind of idea pervades physics today but nowhere is its importance as pronounced 
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as in particle physics. Although the domain of particle physics may ofF-handedly be de- 
scribed as the study of the very small, the range of energies that are scrutinised, from the 
neutrino mass to the unification scale, far outstrips any other area of physics. 

The first great triumph of quantum mechanics, the theory of the very small, was the 
calculation of the spectra of the hydrogen atom. However, the need for corrections to this 
result was soon apparent and there followed calculations for the fine structure. Lamb shift 
and hyperfine structure, each of which relied upon some previously unconsidered higher 
scale physics: spin-orbit coupling, relativistic corrections, loop diagrams in QED and dipole 
interactions. Despite all the corrections made in this calculation there are yet more that can 
be made, such as corrections due to the finite charge radius of the proton. An estimate 
of the magnitude of such a correction can be made in a similar way to the estimate of the 
relativistic corrections above. The corrections are expected to be of the order of ?"^/<3q where 
rp ~ 10~^^m is the charge radius of the proton and uq = 5.29 X 10~^^m is the Bohr radius 
and an estimate of the typical distance of a ground state electron from the proton. Hence we 
expect a correction of the order AE ~ lO'^E. The key point here is the seemingly innocuous 
term 'of the order of. What do we mean by 'of the order of? 

Arguments, like the one above, that involve using a separation of scales and expressions 
like 'of the order of are implicitly assuming 'naturalness'. The assumption of naturalness in 
these arguments is absolutely essential, without it nothing can be said about corrections such 
as relativistic corrections. In a nutshell one might say that the assumption of naturalness 
means that the coefficient, b, of for example the relativistic correction 6 = bv^/c^, is of 
order unity. That is, the number b, which is governed by the dynamics and geometry of the 
problem, is not so large or small to significantly change the size of the correction so that 
expressions like 'of the order of apply. 

Naturalness is not something that can necessarily be taken for granted, and indeed some 
fluke of the parameters may change the face of the problem, resulting in unexpected effects 
such as resonances. 

It is only comparatively recently that physicists have sought to use the separations of 
scales in particle physics quantitatively to produce what have become known as efi'ective 
field theories (EFTs) [1-6]. In an EFT we seek to write down a low-energy equivalent 
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of some 'true' higher energy theory. The concept was first introduced by Weinberg [7] 
in an attempt to use the near chiral symmetry of quantum chromodynamics (QCD) and 
the resulting separation of scales between the pion mass, m„, and all other QCD scales to 
construct a low energy theory equivalent to QCD. Since this initial spark this idea has grown 
to become known as chiral perturbation theory (ChPT). 

ChPT is defined by an eff'ective Lagrangian that contains only those degrees of freedom 
explicit at the energies of interest, Q ~ m;^, i.e. pion fields and depending upon the partic- 
ular problem, nucleon fields, hyper-nucleon fields etc. The fields are coupled by effective 
vertexes that contain the effects of all higher energy physics, of order Aq, not included in the 
Lagrangian, which is said to be 'integrated out'. To include all of the integrated-out physics, 
the Lagrangian must include all terms not forbidden by the symmetries of the higher energy 
physics, invariably resulting in an infinite number of terms. To handle such a formidable 
Lagrangian the terms must be organised according to some 'power-counting' that assigns 
each term an order in an expansion in powers of Q/Aq. Any observable can be calculated 
to arbitrary accuracy by using the power-counting to include all the necessary terms. 

In principle the effective couplings in the ChPT Lagrangian can be found by solving 
QCD, however, in practice they are determined by matching theory to a few experimental 
values, making the theory predictive in other areas. 

In particle physics, the expression EFT was once synonymous with ChPT but it has 
since become an umbrella term for any number of theories in particle and atomic physics 
that subscribe to the same philosophy. Although EFTs are predominantly found in nuclear 
physics, where QCD is non-perturbative and hence almost impossible, they have been ap- 
plied to areas as diverse as quantum gravity and string theory. 

In studying nuclear systems at even lower energies than the pion mass, Q m„, we 
may integrate out the pionic degrees of freedom and deal with the pionless EFT in which 
the only degrees of freedom are the nucleon fields themselves. In this theory the underlying 
physics is that of the pion and the higher energy scale for the theory is the pion mass. The 
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non-relativistic pionless EFT Lagrangian for nucleons of mass M is given by, 

J: = N^ ^ido + ^ j + ^\Nf + y {N^(y - VfN^^N^ + H.c.) + 

^1^1" + ?7 - ^?nAn'^N^ + H.c.) + . . . . 

.36 36 V / 



(1.1) 



It contains four point couplings denoted by C2„ where 2n is the order of the field derivative 
at the vertex. Similarly, since they are not forbidden it contains six point couplings, denoted 
by D2n, and eight point vertexes, ten point vertexes, etc. Relativistic corrections may be 
included by introducing higher order derivatives in the kinetic term [10]. 

Ignoring the six and higher point vertexes, which are not relevant to a two nucleon 
problem, we may immediately construct a power-counting to organise the terms in the ef- 
fective Lagrangian using naive dimensional analysis, which inspired the power-counting in 
ChPT. The four point vertex couplings, C2„, have dimension -2 - 2n. Since they describe 
physics that occurs at the order of the pion mass we may expect. 

In addition each loop between the vertexes contributes a term ~ MQIAn. The order of 
any diagram involving any number of interactions can be evaluated by simply combining n 
vertex terms with the n - I loops in between. A diagram in which the vertex contains v,- 
derivatives will occur at order (Q/m„Y where, 

n 

d = n-l + J^Vi. (1.3) 

!=1 

This power-counting is known as the Weinberg scheme [8]. Unfortunately, when this or- 
ganisation of the terms is used to describe the two-nucleon system it fails. The problem is 
the implicit assumption of naturalness. 

Empirically, we observe in the two nucleon system a low-energy bound state, the 
deuteron, which has binding momentum = 45.7Mev<s; m„. The existence of this and 
the low energy resonances in the other isospin channels alerts us at once to some fine tuning 
of the parameters in the EFT. This fine tuning means that the naive dimensional analysis 
that results in the Weinberg scheme is no longer applicable. Some other power-counting 
must be found if we are to handle the infinite number of terms in the effective Lagrangian 
and produce the low-energy resonances [8-13]. 
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A new power-counting that resolved the problem was first introduced by Kaplan, Sav- 
age and Wise (KSW) [9] using the power divergence subtraction scheme (PDS) but found 
independently by van Kolck [10]. The PDS scheme introduces a new scale, yu, which may 
be chosen to produce the fine tuning eff'ects seen in the two nucleon system. Using the PDS 
scheme the scaling of the coefficients in the effective Lagrangian may be given as [9] 



By choosing /i ~ 2 we introduce the new low-energy scale required. Subsequently the term 
Co occurs at order and so all diagrams containing only the Co vertex must be summed 
to all orders. All other interactions scale with the negative powers of m„ and so may be dealt 
with perturbatively provided each diagram is dressed with Co interaction bubbles. (Since an 
arbitrary number of Co interactions before or after the diagram does not change its order.) 
The subsequent power-counting is given by. 



This KSW power-counting is able to produce the low-energy resonances of order fi seen in 
nuclear systems. 

The origin of these two different power-counting schemes can, in one way, be under- 
stood by using the renormalisation group (RG) [14]. The effective Lagrangian contains, 
by its very definition, an infinite number of non-renormalisable couplings. These result in 
divergent loop integrals which must be regularised, with a sharp momentum cut-off. A, for 
example, and then renormalised by absorbing all cut-off dependence into the effective cou- 
plings. In the context of EFTs the cut-ofF, A, is far more than a UV regulator, it 'floats' 
between the low-energy physics and the high-energy physics and allows us to control the 
introduction of high-energy physics and construct power-counting schemes. After rescal- 
ing all low-energy scales in terms of A, we are led to the concept of Wilson's continuous 
RG [15] initially conceived for use in condensed matter theory. The RG describes how the 
couplings in the theory change as A is varied. 

The use of a sharp momentum cut-off leads us to a very intuitive idea of the RG flow. 
As A gets smaller, more and more of the high energy physics is 'integrated out' and ab- 
sorbed into the couplings of the EFT. Finally as A reaches zero, all high energy physics 



n 




(1.5) 



1=1 
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is removed and all that remains is the rescaled low-energy physics, embodied in an infra- 
red fixed point of the RG. By perturbing about the IR fixed point with powers of A we 
introduce dimensioned constants of integration that must scale with the high-energy scale 
m„. In this way we can re-introduce high-energy physics into the couplings and create a 
power-counting scheme associated with that fixed point. 

In chapter |21 we will study the RG for the two-body pionless EFT and reproduce the 
results of Birse et aZ [14] that reveal two fixed point solutions, a so-called trivial fixed point, 
which is associated with the Weinberg power-counting scheme and a non-trivial fixed point, 
which is associated with the KSW power-counting scheme. The constants of integration in 
the energy dependent perturbations around the trivial fixed point are in one-to-one corre- 
spondence with the terms in the Born expansion. Those around the non-trivial fixed point 
are in one-to-one correspondence with the effective range expansion [16]. 

The two-nucleon pionless EFT is a simple example that demonstrates the use of the 
RG in determining power-counting schemes. However, it is clearly not applicable for the 
interaction of two protons because of their electromagnetic interaction. Photon exchange 
between two protons has the characteristic energy scale, 

aM 

= — = 3.42MeV, (1.6) 

where a is the fine structure constant and M is the mass of the proton. In an EFT for 
protons below the pion mass k must be considered a low energy scale. All single photon 
exchange diagrams must, therefore, be included to all orders. However, as with the KSW 
power-counting, summing some diagrams to all orders and dressing others with 'Coulomb 
bubbles' may completely change the power-counting. 

In chapter 3 we will introduce the distorted wave renormalisation group (DWRG) [17] 
and use it to study, as a first example, proton-proton scattering. The DWRG allows us to 
sum some low energy physics and absorb it into the fixed points of the RG. In particular we 
study a system of two particles interacting via some known long-range potential, equivalent 
to the sum of all known non-perturbative diagrams in the EFT, and a short-range potential, 
equivalent to all the shorter range interactions. By working in the basis of the distorted 
waves (DWs) of the long-range potential we show how to separate the eff"ects of the long 
and short range potentials and then to construct a power-counting for the shorter range 
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interactions. 

The long range potential for the first example in chapter 3 will be the Coulomb poten- 
tial. We will again find two fixed points, a trivial and a non-trivial associated with Weinberg 
and KSW -like schemes respectively. The integration constants in the perturbative series 
about these fixed points are found to be in a one-to-one correspondence with the terms in 
the distorted wave Bom [18] and Coulomb modified eff'ective range [19-21] expansions. 

In the remainder of chapter 3 we will consider two more examples. Firstly, the repulsive 
inverse square potential, which is important to the extension of the pionless KSW EFT to 
three bodies and study of the power-counting in higher partial waves. This example is also 
interesting because it produces novel power-counting schemes that are very difi'erent to the 
schemes that are seen in the Coulomb DWRG and the ordinary RG. 

Secondly, we shall consider a general example of "well-behaved" potentials, which 
includes among many others the Yukawa potential that may be used to include one pion 
exchange (OPE) in a nucleon-nucleon EFT. The solution of this DWRG has much in com- 
mon with the solution of the Coulomb DWRG. This general analysis allows us to construct 
a general method for solving DWRG equations that will prove extremely useful in later ap- 
plications. We conclude the chapter with the discussion of how to introduce explicit pions 
into an EFT for two nucleons. 

After establishing the DWRG we shall look towards using it in the KSW EFT for three 
bodies. As outlined above, the KSW EFT is appropriate for systems with a low energy 
bound state or resonance that constitutes a new low energy scale in the problem. This low 
energy scale is characterised by an unusually large scattering length 1/a ~ 2 «; Aq. Since 
this is a low energy scale its eff'ects must be summed to all orders [22-24] and results in our 
need to use the DWRG. 

Efimov [25] studied systems of three bodies with pairwise interactions characterised 
by a large scattering length, a, and zero eff'ective range. He showed that these systems 
are similar to a two dimensional problem with an inverse square potential. This similarity 
occurs when the characteristic distance between the three particles R ^ a and the problem 
essentially becomes scale free. In the limit of a ^ oo the system is described exactly by an 
inverse square potential. The strength of the inverse square potential is determined purely 
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by the statistics in the system. In the case of three s-wave Bosons and three s-wave nucleons 
in the channel the potential is attractive and singular. The singular nature of these 
potentials has a couple of interesting effects. Firstly, the Thomas effect [26]: the system 
will be no ground state, something first noted by Thomas as long ago as 1935. Secondly, in 
the limit of infinite scattering length, there will be an accumulation of geometrically spaced 
bound states at zero energy, known as Efimov states [25]. 

As the first step towards understanding the KSW EFT for three bodies, in chapter 4 we 
will look at the DWRG for the attractive inverse square potential. The singular nature of 
this potential means that without special measures the Hamiltonian is not self-adjoint. To 
construct a complete set of DWs we must construct a self-adjoint extension to the Hamil- 
tonian [27-29] equivalent to defining some boundary condition near the origin. The con- 
nection to the three-body KSW EFT means that the singular nature of this potential has 
once more come under scrutiny [30-32], in particular with respect to understanding the link 
between the required self-adjoint extension and the three-body force. 

The solution of the DWRG equation for the singular inverse square potential reveals 
that the self-adjoint extension is in fact equivalent to a marginal^ perturbation in the DWRG. 
The general method for solving the DWRG outlined at the end of chapter 3, augmented with 
some special considerations for the bound states, will reveal that the RG flow is controlled 
by limit-cycle solutions that evolve logarithmically in A. 

In chapter 5 we will look at the extension of the DWRG to three body forces (3BD- 
WRG). The 3BDWRG equation is more complex than the standard DWRG because of the 
multi-channelled nature of the system. The 3BDWRG equation has to describe the coupling 
of the three-body force to each of the channels. Fortunately the lessons learnt in chapters 3 
and 4 allow us to find solutions. 

We will demonstrate the solution of the 3BDWRG in an example of "well-behaved" 
pairwise forces and briefly show how the fixed points correspond to Weinberg and an equiv- 
alent KSW counting for three body forces. More interestingly we will look to the derivation 
of the power-counting for the three-body forces in the KSW EFT. 

^The perturbations around fixed points in the RG are described as either stable, scaling with 
positive powers of A, unstable, scaling with negative powers of A, or marginal which do not scale 
with any power of A. 
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The nature of the three-body force in this system is still open to debate. The similarity 
to the singular inverse square potential is clear and the need for a self-adjoint extension is 
almost universally supported [30,32]. It was hoped by some that effective range corrections 
of the two-body force would resolve the singular behaviour, however, results in this direction 
have failed to match experimental data [67]. 

Bedaque et al [22, 24] provide arguments in which they explicitly use a three-body 
force to define a self-adjoint extension and then continue to construct a power-counting for 
it. Phillips and Afnan [33] have shown how the inclusion of one piece of three-body data 
is enough to define all the physical variables without explicitly using a three-body force. 
These related approaches give results that are consistent with each other. Although the 
latter do not specify the actual nature of the three-body force it is clear that the inclusion 
of a piece of three-body data defines a self-adjoint extension of the Hamiltonian and is 
therefore equivalent to the leading-order (LO) three-body force given by Bedaque et al. 

A third view is given by Gegalia and Blankleider [34] who argue that a unitarity con- 
straint is sufficient to define physical quantities with no need for any three-body data at all. 
However, they have failed to produce any results that can be used to check the validity of 
their approach. 

Since three-body force terms are not forbidden by the observed symmetries they must 
be included in the EFT Lagrangian. The 3BDWRG solution for the three-body force shows 
that the self -adjoint extension defining term is equivalent to a LO marginal three-body force 
and, consistent with the results of Bedaque et al [22,24], acts as limit-cycle. Our results 
also support the power-counting suggested by Bedaque et al. 

Having established the power-counting for the three-body force in the pionless KSW 
EFT, in chapter 6 we derive expressions for the DWs in this system at LO and next-to- 
leading order (NLO). The equations for the DWs based upon the two-body interactions 
alone show the singular behaviour anticipated by Efimov. The insertion of the LO three- 
body force is most easily achieved in this equations by demanding some boundary condition 
on the DWs close to the origin. Initially we will derive expressions for a three Boson system 
but in order to produce physically interesting results we will extend the equations to deal 
with a three nucleon system, in particular a system of two neutrons and a proton which 
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contains the two- and three- body bound states, the deuteron and the triton. The LO three- 
body force can be fixed by matching to the neutron-deuteron scattering length. The EFT 
can then be used to construct the triton wavefunction and neutron-deuteron scattering states 
below threshold. 



Chapter 2 

The Renormalisation Group for Short 
Range Forces 

2.1 Introduction 

In this chapter we will introduce the renormalisation group method for constructing EFTs 
[14]. As already observed, the two key ingredients of an EFT are a power-counting scheme 
and a separation of scales. The first of these allows organisation of the infinite number of 
terms that naturally occur in an effective field theory, the second ensures that these terms 
may be truncated to achieve the desired level of precision. 

The separation occurs between the scale of the physics of interest, Q, and that of the 
underlying high-energy physics, Aq. The existence of a separation of scales is usually as- 
sumed with a particular physical system in mind. For example, in nucleon-nucleon scatter- 
ing at energies well below the pion mass, there is a natural separation between the momenta 
of the incoming particles and the lowest energy component of the strong interaction, one 
pion exchange. Utilisation of such a separation leads to what has become known, in nuclear 
physics, as a pionless EFT, in which the only fields are those of the asymptotically free 
particles with no exchange particles [9-1 1, 35, 36]. 

We will consider a system of two non-relativistic identical particle scattering via an 
eff'ective Lagrangian. For weakly interacting systems the terms in the expansion can be 
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organised according to naive dimensional analysis, a term proportional to (Q/AqY being 
counted as of order d. This is known as Weinberg power-counting and is familiar from 
chiral perturbation theory (ChPT) [8]. However, for strongly interacting systems there can 
be new low-energy scales which are generated by non-perturbative dynamics e.g. the very 
large s-wave scattering length in nucleon-nucleon scattering. In such cases we need to 
resum certain terms to all orders and arrive at a new power-counting scheme, KSW power- 
counting [9-11]. The origin and use of these schemes becomes quite transparent using the 
renormalisation group method. The contents of this chapter largely follows the work of 
Birse et al [14,37] 

2.2 The RG Equation 

In a Hamiltonian formulation, the effective Lagrangian is written as a potential consisting 
of contact interactions, which in coordinate space takes the form of 5-functions and their 
derivatives and in momentum space takes the form, 

V(k', k, = Coo + C2oo^" + Co2o^" + Cnok.k' + Qoa/ + • • • , (2. 1) 

where p = '^ME is the on-shell momentum corresponding to the total energy E in the 
centre of mass frame and M is twice the reduced mass. We shall consider s-wave scattering 
only at this stage, then since V must be independent of k.k', Cno = 0. 

Scattering variables are found by solving a Lippmann-Schwinger (LS) equation. In 
particular we shall work with the reactance matrix, K, which is related to the phaseshift, 5 
by, 

1 Mp 

^ cot6, (2.2) 



K(p,p,p) An 
and solves the LS equation, 

T^(T 1 \ ^ui' 1 ^ 2 . V(k',q,p)K(q,k,p) 

K(k ,k, p) = V(k ,k, p) + —r f q dq , (2.3) 

In^ Jo - 

where the bar on the integral sign indicates a principal value prescription for the pole on the 
real axis. Unitarity of the 5 -matrix is ensured if the ^-matrix is hermitian. The on-shell T 
and K matrices are related by 

+ —r^- (2.4) 



K(p,p,p) T(p,p,p) An 



r\ r\ Til. _ T-» x-' _ 
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For short range potentials, it is well known that the inverse of the /^-matrix may be 
written in terms of the effective range expansion, [16] 

47r 1 1 1 9 

-TIT} -=pcoi6 = — + -r,p' + ... (2.5) 

M K{p, p,p) a 2 

where a is called the scattering length and the effective range. The effective range, as its 
name suggests, is related in an indirect manner to the range of the potential [18] (provided 
it is "smooth"). 

For contact interactions the integral in Eqn. (12.31) is divergent. From the EFT point of 
view this is to be expected, the divergence occurs because we have allowed high momentum 
modes to probe the EFT beyond its range of validity. The solution is to remove the high mo- 
mentum modes and incorporate their effects into the effective potential, V. Mathematically, 
we apply a cut-off. A,' to the integral then derive a differential equation for V{k' ,k, p, A) 
based upon the constraint that K, the observable, is independent of that cut-off. Rather than 
acting as a UV regulator the cut-off acts as a separation scale and 'floats' between the low- 
energy scales, Q, and the higher energy scales, Aq. This leads us the concept of a Wilsonian 
renormalisation group (RG) [15]. The interesting limit is A ^ 0, in which all high energy 
physics has been integrated out. 

The regularised LS equation can be written in operator form as, 

K{p) = V(p, A) + V(p, A)G^(p, A)K(p), (2.6) 

where G^ip, A) indicates the free Green's function with standing wave boundary conditions 
and sharp cut-off A. The first step towards an RG equation for V is differentiating Eqn. (12.61) 
with respect to A and then eliminating K to obtain, 

dV(p,A) dGl{p,A) 

^ = -V{p,A) ' v{p,A). (2.7) 

dA dA 

Taking matrix elements and expanding the Green's function this equation takes the form, 

dV(k',k,p,A) M , A2 

= :r-jV(k',A,p,A)— -V(A,k,p,A). (2.8) 



dA 2n^ ^ ' 'A^-p 



^Any type of regularisation and renormalisation should produce the same results, here we just 
take the simple case of a sharp cut-off. 
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To obtain the RG equation for V we must now rescale all low energy scales with A. 
We define the rescaled potential V by, 

V(k',k,p,A) = ^V(Ak',AkAp,A). (2.9) 

where k = k/A, k' = k' /A and p = pi A. The resulting RG equation is, 

dV dV ^dV ^,dV . . 1 . . 

= + ^Tf + + ^ + hp,A)—-:^V(l,k,p,A). (2.10) 
OA dp dk dk' 1 - P^ 

The reason for this rescaling is that it allows distinct treatments of the low-energy physics 
and the parameterisation of high-energy physics. To explain, let us consider solutions that 
are independent of A. These solutions, known as fixed points, are scale free since they are 
independent of the only scale in the problem. Therefore they must be related to the rescaled 
low energy physics and be independent of the unsealed high energy physics. 

Fixed-point solutions provide systems in which only the low-energy physics has been 
included and it is this that provides the key to constructing power-counting schemes. Solv- 
ing the RG equation in the region of a fixed point by perturbing about it in powers of A 
introduces parameters that must scale in inverse powers of Aq, the scale of the high en- 
ergy physics. The perturbations around a fixed point organise the efi'ective couplings into a 
power-counting scheme. 

Since a fixed point solution is independent of A, all solutions of the RG equation tend to 
one as A — > 0. We may classify fixed points by examining the stability of the perturbations 
about them as A ^ 0. If all perturbations tend to zero as A ^ 0, i.e they scale with positive 
powers of A, the fixed point is stable. If there are perturbations that scale with negative 
powers of A then the fixed point is unstable. A perturbation that does not scale with A 
is known as marginal, and as we shall see in later chapters is associated with logarithmic 
behaviour in the cut-off. 

In this problem, we shall find two fixed points that are of particular interest. The first 
of these, which we shall refer to as the trivial fixed point, is simply the obvious solution 
V = 0. This solution corresponds to the scale free system in which no scattering occurs. We 
shall see that the power-counting associated with the fixed point is the Weinberg scheme. 
A second fixed point, which we shall refer to as the non-trivial fixed point, gives a scale 
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free system in which there is a bound state at exactly zero energy. The power-counting 
associated with this fixed point is the KSW scheme. 

Eqn. (12.11) and the assumption of s-wave scattering implies that V should be an analytic 
function of p^,k^, and k'^. These constraints constitute a boundary condition upon the RG 
equation, physically they follow because the energy is below all thresholds for production 
of other particles and that the effective potential should describe short-ranged interactions. 



2.3 The Trivial Fixed Point 

The Trivial fixed point solution is the A-independent solution V = 0. To examine the 
perturbations around it we write the potential as 

Vik, k', p. A) = CA^cpik', k, p), (2.11) 

insert it into the RG equation, linearise, and obtain the eigenvalue equation 

p^+k^+k'^ + (P = v(p. (2.12) 
dp dk dk' 

The solutions of this equation are easily found to be 

(f)(k',k,p) = k'^'B'"p^", (2.13) 

with RG eigenvalues v = 2(1 + m + n) + \. The analyticity boundary conditions imply that 

/, m and n must be non-negative integers, so that v takes the values 1, 3, 5, The solution 

to the RG equation in the region of the trivial fixed point is 

V(k',k,p,A)= y qJ^-] k'H'"'p'", (2.14) 

where the coefficients C/„„, have been made dimensionless by extracting the high-energy 
scale Aq'' and for a Hermitian potential we must take Ci^n = Cmin- Since all the RG eigen- 
values are positive this fixed point is stable as A ^ 0. 

The RG eigenvalues, y, provide a systematic way to classify the terms in the potential. 
In unsealed units we can see how this provides us with a power-counting scheme and straight 
away identify the couplings in the naively dimensioned potential suggested in chapter 1 : 

2n^ A . k'^'k^'-p^" 
V(k',k,p,A) = > Cto„ ^ , ■ (2.15) 

MAn ^ A2(/+m+n) ^ ^ 

^ l,m,n=0 ^^0 
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The power-counting scheme is the Weinberg scheme if we assign an order J = v - 1 to each 
term in the potential. That this power-counting scheme is useful for weakly interacting sys- 
tems now comes as no surprise considering its association with the trivial, or zero-scattering, 
fixed point. 

The corresponding ^-matrix is simply given by the first Bom approximation, as higher- 
order terms from the LS equation are cancelled by higher-order terms in the potential from 
the full, nonlinear RG equation [14]. It is important to notice that terms of the same order 
in the energy, p^, or momenta, or k'^, occur at the same order in the power-counting. 
It is therefore possible to swap between energy or momentum dependence in the potential 
without aflFecting physical (on-shell) observables. 



2.4 The Non-Trivial Fixed Point 
2.4.1 Momentum Independent Solutions 

Let us consider a fixed point solution, Vq{p), that depends on the energy, p, but not upon the 
momentum. It satisfies the RG equation 

p-r^ + yo + -T—Ti = 0- (2-16) 

dp I - p^ 

A convenient way to solve this equation, as well as other momentum-independent RG equa- 
tions, is to divide through by and obtain a linear equation for 1 /Vq, 



This equation is satisfied by the basic loop integral. 



= 0. (2.17) 



-f 

Jo 



l-^lni±^ 
2 1-p 



(2.18) 



p^-q^ 

Since J(p) is analytic in p^ as p^ —> it is a valid solution to the RG equation. Hence we 
take Vo = 1 /J- We shall refer to this as the non-trivial fixed point. When Vq is inserted into 
the LS equation we obtain 

K(p,p,pr'=0, (2.19) 
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which corresponds to a system with a bound state at zero energy. Since the energy of the 
bound state is exactly zero, the system has no scale associated with it, which is why it is 
described by a fixed point of the RG. 

To find the power-counting associated with this fixed point we look for perturbations. 
This task is considerably easy if we consider momentum independent solutions. We con- 
sider perturbations of the form 

1 1 



giving the eigenvalue equation. 



+ CAXp), (2.20) 



dd> 

P^-cp = v(l>. (2.21) 
op 

Notice that since the momentum independent RG equation for 1/ V is linear, no approxima- 
tion has been made to obtain eqn. (l2.21l) . The potential obtained using these perturbations 
will be an exact solution to the RG equation. The eigenvalue equation is easy to solve: 

m = P^\ (2.22) 

with RG eigenvalues v = In - \. The analyticity boundary conditions demand that n is 

a non-negative integer, so that the eigenvalues are v = -1,1,3, The first eigenvalue 

is negative, which means that the fixed point is unstable as A ^ 0. The full, momentum 
independent solution is given by 

+ y,C2.-i 4^ p'". (2.23) 



V{p,K) Voip) \A. 



^0 



The power-counting about this fixed point is the KSW scheme. As in the Weinberg case we 
assign each term an order J = v - 1, so that the KSW scheme is characterised by the values 
d = -2,0,2,.... 

When this potential is inserted into the LS equation, we obtain an expression for the 
K-matrix, which upon using eqn. (l2.2l) gives the phaseshift, 

pcotS = -^-^yC,,J{-T. (2.24) 

We can see that the terms in the expansion around the non-trivial fixed point are in one-to- 
one correspondence with the terms in the effective range expansion (12.51) . In particular, we 
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b„(A) 



Figure 2.1: The RG flow of the general momentum independent solutions, V(p,A) = 
2^o^2«(A)p^". The fixed-points are shown as dots. The eigenvector flows are shown in 
bold with the arrows indicating flow as A ^ 0. More general flows are shown as dashed. 

find 

C-1 = -— , Ci = f^. (2.25) 

If the theory is natural then we expect the coefficients €2,1-1 ~ 1- In this case the scattering 
length and the effective range are of order a ~ ~ I/Aq. 

To understand which of the two power-counting schemes is appropriate in a given 
situation, it is fruitful to examine the RG flow of purely energy dependent RG solutions 
illustrated in Fig l2.1l This figure shows the RG flow as A ^ in the plane (bo{A), ^2(A)) 
where 

V(P, A) = ^o(A) + b2(A)p^" + .... (2.26) 

The trivial fixed point lies at (0, 0). All flows near to it flow into it as A ^ 0, illustrating its 
already noted stability. The non-trivial fixed point lies at (-1, -1), only flows lying on the 
critical line bo = -I flow into it. All other flows in the region of this critical line flow initially 
towards the non-trivial fixed point before the unstable perturbation becomes apparent and 
takes the flow into the domain of the trivial fixed point (possibly 'via' infinity). 
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The values of A for which the flow is controlled by the non-trivial fixed point is depen- 
dent on the value of the coefficient C_i . In particular if 

A>^= AoC_i (2.27) 
2a 

then the unstable perturbation is suppressed and the RG flow is controlled by the expansion 
around the non-trivial fixed point. As A becomes less than n/(2a) then the unstable pertur- 
bation becomes important and the flow moves away from the non-trivial fixed point and into 
the trivial fixed point. 

Li a natural theory we have 1 /a ~ Aq so that as A floats between the scale of the high 
energy physics, Aq, and 0, the condition (12.271 ) is not met. This means that the RG flow is 
in the domain of the trivial fixed point and that the appropriate power-counting is the one 
obtained by perturbing about that point, the Weinberg scheme. 

If there is some fine tuning in the theory so that C-i is small, then the unstable per- 
turbation is suppressed. This fine tuning, according to eqn. (I2.25I) . will reveal itself in a 
surprisingly large scattering length, a » 1 / Aq. As A varies from Aq to 1 /a the condition 
(12.271) is true and the power-counting associated with the non-trivial fixed point provides a 
suitable expansion. 

For example, in s-wave nucleon-nucleon scattering the scattering lengths are large^ 
suggesting a finely tuned effective potential with a small coefficient C-i, in this case the RG 
flow is dominated by the non-trivial fixed point and the power-counting associated with it, 
the KSW scheme, is the appropriate one to use in organising the terms [9-1 1, 35, 36]. 

2.4.2 Momentum Dependent Solutions. 

In the analysis above the Weinberg and KSW counting schemes occur quite naturally. Study 
of the RG flow leads to conclusions about the usefulness of each of the schemes in different 
systems. However, the analysis is not complete as we have not considered the momentum 
dependent solution about the non-trivial fixed point. 

Around the trivial fixed point, it was found that the momentum dependent perturbations 

provided no new physics since they occurred at the same order in the power counting as the 

^Naturally the scattering length for NN scattering should be of the order of the pion mass, a ~ 
l/m„, however this is not what is empirically observed. 
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energy dependent perturbations. This meant that momentum and energy perturbations could 
not be distinguished within on-shell observables. 

Around the non-trivial fixed point things are considerably more complicated, as we 
shall see the momentum and energy dependent perturbations occur at different orders in the 
power counting. This result means that it is not a simple matter to interchange momentum 
and energy dependence of physical observables near the non-trivial fixed point. 

To find the momentum dependent perturbations we write the solution in the vicinity of 
the fixed point solution as, 

V(k, k', p. A) = Voip) + CAXi k\ p). (2.28) 
After linearising, the eigenvalue equation for ^ is 

# ^ .50 ^ ^ Vo(^ r^^^ ^ ^ (f>(l,k',p)\ = (V - 1)^. (2.29) 
op dk dk' ^-P^ ^ 

Useful solutions to this equation [37] are 



^n{k,p) = 



Voip), (2.30) 



with eigenvalues, v = 2n = 2, 4, 6, . . .. These solutions lead to two forms of perturbation that 
are symmetric in k and k' and so may be used within an hermitian momentum dependent 
potential: 

(f>(k,k',p) = p^ l<P„(k,p) + <Pn(k',p)} , (2.31) 
with eigenvalues v = 2m + 2n = 2,4,6,... and 

<P(k, k', p) = p^'" {cpnil p)<Pn4',p)] , (2.32) 

with eigenvalues v = 2m + 2n + 2n' + I = 5,7,9, .. .. These perturbations plus the the 
purely energy dependent perturbations form a complete set that can be used to expand any 
perturbation about the non-trivial fixed point that is well-behaved as k^,k'^,p^ 0. The 
momentum dependent perturbations vanish on-shell and the terms in the on-shell potential 
are still in one-to-one correspondence with the terms in the effective range expansion [37]. 
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The most general solution in the vicinity of the non-trivial fixed point may be written as 



V(k, k', p. A) = Vo(p) + g C2n-i |— j P^"Vo(pf+ 

2 D2(„+„) — p^"[cpUk,p) + <Pm(k',p)] + 

n=0,m=l ^ ^' 

2 £2(n+m+m')+l It-) p^" (Pm{k, P)^m'{k' , p) . (2.33) 

n=0;m,m'=l ^ ^' 

When the on-shell ^-matrix is calculated using this general form of expansion about the 
non-trivial fixed point, the only terms that contribute are precisely those that appear in the 
expansion given in the previous section, namely the coefficients of the purely energy de- 
pendent eigenfunctions, Cy. The momentum dependent eigenfunctions do not contribute 
to on-shell scattering, i.e. physical observables are independent of and Ey at this or- 
der [14,37]. 

The solution given above results from linearisation of the RG equation close to the non- 
trivial fixed point. It is not clear that if a full solution to the non-linear RG equation was 
constructed that the one-to-one correspondence betweens terms in the potential and terms 
in the effective range expansion persists. In order to clarify the situation, Birse et al [14] 
calculated corrections to the potential around the non-trivial fixed point due to the neglected 
non-linear terms in the RG equation up to order and found the additional term, 

— C_iZ)2 + k'+Ap^ + — Voip)\ Voip), (2.34) 

where A is a constant of integration. In general this second order piece will contribute to the 
eff'ective range term along with the Ci term. Since the constant A is unfixed, we are free to 
set it as we wish to ensure a nice correspondence between terms in the potential and terms in 
the effective range expansion. Setting A = -2 maintains the direct correspondence between 
the effective range expansion and the energy dependent perturbations around the non-trivial 
fixed point. In general it is assumed that at each order in the coefficients further degrees of 
freedom will become available to ensure this correspondence is continued to all orders. 

The degrees of freedom that become available as the full non-linear solution is explored 
allows the possibility of generating the effective range terms from momentum dependent, 
rather than energy dependent, perturbations. For example, it is possible to generate the 
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effective range, r^, solely from the momentum dependence in the correction given in eqn. 
(12 .341) . By setting Ci = and A = the effective range would be given by, 

re = (2.35) 
Ala 

However, since the momentum dependent perturbations occur at different RG eigenvalues, 
V = 2, 4, . . ., than the energy dependent perturbations they are replacing, the coefficient Dy 
is forced to take an unnaturally large value to compensate for the additional factor of Aoa. 

Unlike the trivial fixed point the possibility of using the equations of motion to move 
between momentum and energy dependence is not manifest at any particular order in the 
power-counting. This trade has to be between terms resulting from non-linear corrections 
in the solution to the RG equation in the vicinity of the fixed point, making it difficult to see 
and implement. Importantly, the momentum dependent perturbations offer no new degrees 
of freedom to on-shell observables than the easily elucidated energy dependent ones. 



2.5 Summary 

In this chapter we have introduced the RG for short range forces and used it to derive the 
Weinberg and KSW power-counting schemes. The key to solving the RG equation are fixed 
point solutions. Perturbing around the fixed point solutions gives us a simple recipe for 
constructing power-counting schemes. We have identified two fixed point solutions, the 
trivial fixed point that leads to the Weinberg scheme and a non-trivial fixed point that leads 
to the KSW scheme. 

The terms in the perturbations around the trivial fixed point are in one to one correspon- 
dence with the terms in the DW Bom expansion. Those about the non-trivial fixed point are 
in one-to-one correspondence with the terms in the effective range expansion. 

By studying the RG flow with the cut-off A we may make conclusions about the use- 
fulness of each of the fixed points. The non-trivial fixed point is unstable and requires the 
fine tuning of a parameter in the expansion around it. This fine tuning results in large scat- 
tering lengths. Because of the large scattering lengths observed in nuclear systems the KSW 
counting is the appropriate scheme. 



Chapter 3 



The Distorted Wave Renormalisation 
Group 

3.1 Introduction 

In this chapter we shall introduce the distorted wave renormalisation group (DWRG) [17]. 
This is an extension of the RG which allows the summation of some physics to all orders. 

The ability to construct an EFT with non-trivial low-energy physics is important. Con- 
sider proton-proton scattering at low energies. The Coulomb interaction between the pro- 
tons is a very long-ranged interaction with a characteristic momentum scale, k = ^ = 
3.42MeV. If we are interested in scattering between protons of typical energy p ~ k,k must 
be treated as a low-energy scale, so we must sum all photon exchange diagrams since they 
occur at LO in the EFT [21]. 

Proton-proton scattering is not the only example that is important to nuclear physics. 
Also of interest is nucleon-nucleon scattering at energies comparable to the pion mass [37- 
39, 45]. In the previous chapter, we constructed a EFT in which all exchange particles 
are absorbed into the effective couplings. In nucleon-nucleon scattering the scale at which 
this breaks down is the pion mass, m„, which acts as a high energy scale in the EFT. If 
we wish to examine nucleon-nucleon scattering at energies comparable to the pion mass, 
pion fields must be introduced into the EFT Lagrangian. Dependent upon our identification 
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of low-energy scales, one approach may be to sum all one pion exchange diagrams [45]. 
There are many issues surrounding how pions are to be including in an EFT for nucleons 
[8,40,45,49,50], we shall consider these towards the end of this chapter. 

We will introduce the DWRG equation by showing how the short and long range 
physics can be neatly separated. Then we will apply the equation to a number of exam- 
ples. Our first example, in which the long-range physics will be modelled by the Coulomb 
potential, will demonstrate the methods that will be used in later examples. In the DWRG 
analysis of this system we will derive a known result for the distorted wave effective range 
expansion [19-21]. 

Our second example, will be the scale-free, repulsive inverse square potential. Our 
interest in this example is two-fold. Firstly, it will allow examination of short range forces 
for higher partial waves in the EFT for short range forces and secondly, the scale-free nature 
of the potential makes it very similar to three body systems [25]. 

As a final example we shall consider a more general class of long-range forces, namely 
non-singular potentials that facilitate the definition of the Jost function [18]. This general 
analysis is of interest as it includes the Yukawa potential and because the methodology used 
in this section acts as inspiration for the DWRG analysis in later chapters. 

3.2 Separating Short- and Long-Range Physics 

We shall assume that all non-perturbative EFT diagrams can be summed to give a long range 
potential. To this end, we consider a system of two-particles of mass M interacting through 
the potential. 



where Vs is an effective short-range force that consists of contact interactions only. The 
terms in V5 can be related to the couplings in the effective Lagrangian. 

The full T-matrix that describes scattering from both the long and short range potentials 
is given by a LS equation. 



V = Vl + Vs, 



(3.1) 



Tip) = (Vl + Vs) + (Vl + Vs)G^,(p)T(p). 



(3.2) 
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We wish to find an RG equation that allows us to determine the power-counting for the 
short range force. Suppose that we were to attempt to find this by applying the cut-off to 
the free Green's function in eqn. (I3.2I) . The resulting differential equation for Vs, found by 
demanding that T is independent of A, is 

dVs ~ dGt 

-1 = -(Vl + Vs)^(Vl + Vs). (3.3) 
dA dA 

After taking matrix elements we obtain an equation that in general contains complicated A- 
dependence in terms that are quadratic, linear and independent of Vs ■ In essence this com- 
plicated A-dependence is because of the cut-off, which applied to the free Green's function 
in the LS equation, not only regularises the divergent loop integrals between contact inter- 
actions but also cuts-off elements of the long-range potential. Those parts of the long-range 
potential "removed" by the cut-off then have to be absorbed into the short-range potential 
resulting in the equation above, it then becomes difficult to justify any boundary condition 
on the potential, Vs . 

To circumvent this complication and obtain an RG equation containing the short-range 
potential alone, it is useful to work in terms of the distorted waves (DWs), |i/^p), of the long- 
range potential. The DWs are simply solutions to the Schrodinger equation containing the 
long-range potential alone: 

/^o + Vi-^j|<Ap) = 0. (3.4) 
The T-matrix, Tl, describing scattering by Vl, is simply given by the LS equation: 

Tdp) = Vl + Vi^Gl{p)Tdp). (3.5) 

The corresponding full Green's function is 

Glip) = Gtip) + Gt{p)T:^{p)Glip). (3.6) 

To isolate the effects of the short-range potential and enable a renormalisation group 
analysis we write the full T -matrix as [18] 

T{p) = Tdp) + (1 + Tdp)Gl{p))ts{p){\ + Gl{p)Tdp)). 0.1) 

The operator (1 -l- G'^Tl) is the MoUer wave operator that converts a plane wave into a DW 
ofVL, 

= (1 + G^^(p)Tdp))\v), (3.8) 
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so that the matrix elements of the full T-matrix are given by, 

{k\T(p)\k') = <k|rz.(7?)|k') + {r^\fs{p)K)- (3.9) 

An equation for the operator fsip) can be found by substituting eqn. (13.71 ) into the full 
LS equation (I3.2I) . After identifying and cancelling the terms in the pure long-range LS 
equation and cancelling the MoUer wave operator on the right hand side we are left with, 

[1 + Tl(p)G^o(pWs(p) = Vs+ VsG^(p)[l + Tdp)G^o(pWs(p) 

+VlG;{p)[1 + Tdp)Gl{pWs{p). (3.10) 

The third term on the right-hand side can be seen to cancel with the second on the left-hand 
side after identifying the LS equation for Ti. The second term on the right hand side can be 
simplified by identifying the form for the full long-range Green's function, eqn. (13.61) . What 
remains is a distorted wave Lippmann-Schwinger (DWLS) equation for the ts, 

Ts(p) = Vs + VsGl(p)fs(p). (3.11) 

This equation is now a far more promising starting point for the RG analysis since the effects 
of the long-range potential have been dealt with separately. The operator describes the 
interaction between the short-range potential and the DWs of the long-range potential 

We shall assume that Vl is a central potential and work in the partial wave basis. In 
particular we shall concentrate on s-wave scattering. In that case the on-shell T-matrix may 
be expressed in terms of the phaseshift,^ 

(pIWIp) = . (3.12) 

A similar relation holds between and 6t, the phaseshift due to the long range potential 
alone. We may write the full phaseshift, 6, as 6i plus a correction, 6s, due to the effect of 
the short-range potential: 

5 = 5i^ + ~5s. (3.13) 

^Note that the normalisation used in this chapter is slightly different from that implicitly assumed 
in the previous chapter, hence the relation between the J-matrix and the phaseshift differs by a factor 
ofp2. 
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Substituting this form into eqn. (13.121) and subsequently into eqn. (I3.7I) we obtain a simple 
relationship between the on- shell matrix elements of and the corrective phaseshift ^5 , 

4-TTn «2!(Ss(p) _ 1 

3.3 The DWRG equation 

In order to derive the DWRG equation for the short-range potential it is once again more 
convenient to work with a reactance matrix, Ks , which satisfies the DWLS equation, 

Ksip) = Vsip) + Vsip)Giip)Ksip), (3.15) 

where in this equation indicates the Green's function with standing wave boundary con- 
ditions. The relationship between the on- shell matrix elements of f and K (c.f. eqn. (12.41 ) 
is, 

1 e2,-<5z.(p) 

+ -—, (3.16) 



(i(,p\Ks(.pMp) (rplfsipMp 4;rp 
where similarly li^/p) without superscript indicates standing wave boundary conditions on 
the DW. Hence the on-shell ^5 -matrix is given by 

1 M 
= = -— cot^s. (3.17) 

i^fplKsipMp) 4np 

To regulate eqn. (13.151) we expand using the completeness relation for the DWs, and 
apply a cut-off to the continuum states, 

<^L(i^' A) = — + dq— J + — > . ^ . (3.18) 

Jo p -q ^^"^ p + Pi 

Demanding that Ks is A-independent we can obtain, as in the previous chapter, a differential 
equation for Vs by differentiating eqn. ( 13.151 ) with respect to A and eliminating Ks to obtain, 

— f (;?,A) = -Vs{p,A)—^{p,A)Vs{p,h). (3.19) 
oA oA 

In the remainder of this chapter, to simplify the analysis, we shall assume that the 
matrix elements of V5 depend only on the energy, p and not on the momentum. As shown 
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in chapter |2l the off-shell momentum dependent solutions to the RG equation have more 
complicated forms but are not needed to describe on-shell scattering [14, 37]. 

In the previous chapter we simply used a contact interaction proportional to a delta 
function, however, such a choice cannot be used in combination with some of the long-range 
potentials of interest. Long-range potentials that are sufficiently singular that their DWs, 
il/p{r), either vanish or diverge as r ^ will result in a poorly defined DWRG equation if a 
delta function contact term is used. 

To construct a contact interaction we chose a spherically symmetric potential with a 
short but nonzero range. By choosing the range, R, of this potential to be much smaller 
than 1 /A, we ensure that any additional energy or momentum dependence associated with 
it is no larger than that of the physics which has been integrated out, and hence the power- 
counting is not altered by it. The precise value of this scale is arbitrary and so observables 
should not depend on it, we may consider the results to be equivalent to those that would 
be obtained in the limit R ^ 0. R should be thought of nothing more than a tool to avoid a 
null DWRG equation. A simple and convenient choice for the form of the potential is the 
"5-sheU" potential, 

Vs(r) = VsS(r-R). (3.20) 
With this choice, the eqn. (13.191 ) for Vs becomes 

OA = -i;^^^^^^^''^'^^- ^'-^'^ 

The final step in obtaining the DWRG equation is to rescale each of the low energy 
scales, expressing them in terms of the cut-off" A. Dimensionless momentum variables are 
defined by ^ = p/A etc., along with a rescaled potential. Vs. The exact nature of the 
rescaling required for the potential is dependent on the form of i^a(R)- 

The solution to the DWRG equation should satisfy the analyticity boundary conditions 
in that follow from the same arguments given in the previous chapter. Vs must also be 
analytic in any other low energy scale, k, associated with the long-range potential. The exact 
analyticity condition for each scale k is case specific. In most cases we need to demand only 
that the eff'ective potential is analytic in k. An example is the inverse Bohr radius which is 
the low-energy scale associated with the Coulomb potential, and which is proportional to 
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the fine structure constant, a. Since the short-distance physics which has been incorporated 
in the effective potential should be analytic in a it should be analytic in k. An important 
exception is the pion mass. This is proportional to the square root of the strength of the 
chiral symmetry breaking (the current quark mass in the underlying theory, QCD) and so, 
as in ChPT, the effective potential should be analytic in m^. Under these restrictions, we see 
that Vs should have an expansion in non-negative, integer powers of and k (or k^). 



3.4 The Coulomb Potential 

Having discussed the general route to the DWRG equation we shall now provide a con- 
crete and physically interesting example to see how the whole procedure fits together. The 
Coulomb potential, Vtir) = ar~\ is highly interesting physically for obvious reasons. The 
wavefunction for the Coulomb potential may be written in terms of the confluent hypergeo- 
metric function, 0(a, b, z), 

^(;)(r) = pre-5'^T(l + /77)0(1 + ij], 2, -2ipr)e'P\ (3.22) 

where rj = k/p, k = aM/2 [18]. The definition of the phaseshift for the Coulomb potential 
is unusual because the long tail results in logarithmic asymptotic contributions. The 
Coulomb phaseshift is defined by writing the asymptotic behaviour of the wavefunction as, 

i/rW(r) ^ e'^' s'mipr -r]\n2pr + 6^), (3.23) 

which gives the phaseshift 6c = Ar^{r(l -I- ir])}. Incidentally, this expression also defines the 
normalisation of the DWs. A simple calculation shows that for R «; A"' we have, 

IMRf - A^^^.-'t(i + ,-^)r(l - ,|) = = R^A^CMA). (3.24) 

where C(ri) is known as the Sommerfeld factor. Hence, in the Coulomb case eqn. (13.211) 
becomes 

dA = -^^^'l^^^-pT^^ysiP,'^^)- (3.25) 
To obtain the DWRG equation we must rescale the variables p = pA and k = kA. Vs must 
be rescaled to factor out the scales M and R. Hence we define, 

MAR^ 

Vsip,k,A) = -^Vs{Ap,Ak,A\ (3.26) 
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resulting in the DWRG equation for the Coulomb potential: 

dVs dVs dVs ^ C(k) 

= p—1 + + Vs + -^^Vc. (3.27) 
OA ^ dp dk ^ \- ^ ^ 

The boundary conditions that follow from the discussion above are analyticity about p^,k ^ 

0. We shall assume that k is positive, i.e. that the Coulomb potential is repulsive, the mod- 
ifications to this discussion for an attractive Coulomb potential are examined in Appendix 

m 

Our study of the Coulomb DWRG equation will mirror the RG analysis of the previous 
chapter. In particular we shall concern ourselves with fixed-point solutions and their corre- 
sponding power-counting schemes. A trivial fixed point, Vs = 0^ is readily identified, but 
the existence of a non-trivial fixed point is not obvious. We will see that no true fixed point 
other than the trivial one exists. However, we will find a logarithmically evolving 'fixed 
point' that for the purpose of constructing a power-counting scheme may be regarded as a 
true fixed point. Furthermore, we will find a marginal perturbation in the RG flow about this 
fixed point that is associated with its logarithmic evolution. 

3.4.1 The Trivial Fixed Point 

An obvious fixed point solution of the DWRG equation for the Coulomb potential is the triv- 
ial fixed point solution, Vs = 0. This solution corresponds to a vanishing short-range force, 

1. e. a system in which only the Coulomb force is active. The power-counting associated 
with the point can be obtained by solving the eigenvalue equation for a small perturbation 
about it. Putting 

Vs{p,k,A) = CA''((>ip,k), (3.28) 
we obtain an eigenvalue equation that is very similar to that in the previous chapter, 

+ K— = (v - (3.29) 

op OK 

The solutions that satisfy the boundary conditions are = p^"k^, where n and m are positive 

integers. The RG eigenvalues are v = 2n + m + 1 = 1, 2, 3, All the RG eigenvalues, 

in common with the pure short range case of the previous chapter, are positive and the 
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fixed-point is stable. The potential in the vicinity of the fixed-point is given by, 

oo / . \2n+m+l 

ys=y]C2nJi^] p'"}r. (3.30) 

If we assign d = v - I then the power-counting scheme is the Weinberg scheme with ad- 
ditional /c-dependent terms. Upon substitution into the Lippmann Schwinger Equation, all 
higher order terms vanish and what remains is a distorted wave Born Expansion, 

— = -^^(^)L^--b^ ' (3.31) 

n,m=i ^ U ' 

which, it should be noted, is independent of the delta-shell range, R. This power-counting 
scheme, in parallel to the trivial fixed point for purely short-range potentials, is appropriate 
for systems with weakly interacting short-range potentials that provide only minor correc- 
tions to the Coulomb potential. 



3.4.2 The Non-Trivial Fixed Point 

The starting point for finding another fixed point solution is to divide the DWRG equation 
through by to obtain a linear PDE for 

dV,^ dV,^ 5t>-i . , C(k) 

Substitution shows that this equation has a fixed point solution given by the basic loop 
integral, 

J(p,K)= f q^dq^P^. (3.33) 

Jo p - q 

However, we may not simply take Vs = as the solution to the RG equation because it 
does not satisfy the necessary analyticity boundary conditions. A source of non-analyticity 
in are the poles atq = +p. As p —> these poles move to the endpoint of the integral, q = 
0, causing singular behaviour in /. A second source of singular behaviour is the essential 
singularity in C(ri) at 77 = 0. This singularity means that the integral J only converges 
for Re{^} > resulting in singular behaviour about ^ = 0. In order to isolate these non- 
analyticities we write / as: 

J(p, k) = - dqCiklq) - dq^^-^ + M(p, k), (3.34) 

Jo Ji p^ - r 
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where 

^2 



Jo 



M(p,k) = p''jr dq^P%. (3.35) 
p - q 



The non-analyticity in p^ caused by the q = endpoint now appears in M., whilst the non- 
analyticity resuking from the essential singularity in C is contained in both M. and the first 
integral in eqn. (I3.34I) . The second integral in eqn. (13.341) avoids the troublesome endpoint, 
4 = 0, and is analytic in both k and p. The first integral in eqn. (13.341) may be written as (see 
Appendix lAb. 

-1 

dqC(k/q) = -1 - tt^ In ^ + Analytic terms in k, (3.36) 



Jo 



So that overall we have, 

J(p, k) = Mip, k) - nklnk + Terms Analytic in p^, k. (3.37) 

To construct a solution to the DWRG fixed point equation that satisfies the analyticity 
boundary conditions we must remove the non-analytic terms M and kink. The removal 
of the non-analytic term M is simple. It is not difficult to show that it satisfies the homoge- 
neous form of eqn. (13. 321) . 

pf^.f.M^O. (3.38) 

op OK 

SO that J - M satisfies eqn. (13.321) . However, the removal of the logarithmic term in 
eqn. (13.361) is not as simple and cannot be done within the confines of the fixed point equa- 
tion. It is because of this term that we are forced to introduce logarithmic A dependence 
into the fixed point solution. The term, 

£ik,A) = -nkln—, (3.39) 

where yu is some arbitrary scale, satisfies the homogeneous form of the eqn. (13.321) . 

^^-^^+^ = 0' (3.40) 
dA dK 

and so may be used to remove the logarithmic term from 7. Bringing all the terms together 
we take 

Vf\p, k. A) = (j\p, k) - M{p, k) - l{k, A))~' (3.41) 
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as an analytic solution to the DWRG equation with logarithmic dependence on A. We will 
refer to this as the non-trivial fixed point ^. 

The perturbations around the non-trivial fixed point can be found in the same way as 
before. Writing, 

^ ^ +h^(P{p,k), (3.42) 



we obtain a linear equation for (p which is readily solved, 

(Pip,k) = P^"ir, (3.43) 

where n and m are positive integers and the RG eigenvalues are v = 2n + m- I. The leading 
order perturbation around this fixed point has a negative RG eigenvalue and so, like the non- 
trivial fixed point seen in the previous chapter, is unstable. Indeed, the power-counting in 
the energy dependent terms around this fixed point is simply the KSW scheme. In contrast 
to the power-counting observed there, the existence of a zero RG eigenvalue means that 
there is a marginal perturbation that does not scale with a power of A. The full solution in 
the vicinity of the non-trivial fixed point is, 

1 1 °° / * \2n+m-l 

1 1 . A 



oo I , \zn+m—i 

+ ^C2„,™|— J p-"r\ (3.44) 



The marginal perturbation is the key to understanding the need for logarithmic dependence 
on A in the fixed point solution and the arbitrary scale fi. Since the marginal perturbation is 
independent of A it cannot be separated unambiguously from the fixed point solution . 
The degree of freedom associated with the scale p is interchangeable with the coefficient 
Coi • Indeed the coefficient, Cqi can be chosen to depend on p in such a way that the term, 

A . 

nkln - + CQi{p)k, (3.45) 
)" 

and hence the full solution is independent of p. 

This system is very similar to the case of considered in ChapterEl The power-counting 
in the energy dependent terms around both fixed-points is the same as observed in the system 
with just short range forces. As a consequence, the RG flow and the conclusions drawn from 



^Although not strictly a fixed point, it shall be referred to as such for ease of nomenclature. 
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its examination are very similar. Writing, 

CO 

Vsip,k,A) = b2n,Mp''"r, (3.46) 

n,m=0 

the flow in the (Z?o,o(A), Z?2,o(A))-plane is as given in Fig. 12. IP . In that diagram the logarith- 
mic behaviour associated with the non-trivial fixed point is not apparent but is illustrated in 
Fig l3.ll which shows the flow in the plane (Z7o,o(A.), Z?o,i(A)). As before, the fixed point is 
shown as a dot, the flows along the RG eigenvectors as bold lines with the arrows showing 
the flow as A ^ and the dashed lines showing more general flow lines. The flow associ- 
ated with the marginal perturbation carries the potential up the vertical flow line Zpq o = -1 at 
a logarithmic rate. Although the coefficient appears to be tending to infinity along this flow 
line, it will eventually become so large that the expansion (13.461) breaks down. A general 
potential, for which Zjq.o is not exactly -1, will flow into the trivial fixed point. However, 
provided the coefficient of the unstable perturbation, Co,o is small, it is still possible to ex- 
pand around the logarithmic flow line. 

^The exact position of the non-trivial fixed point may be different but this does not affect the 
conclusions. 
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Despite the logarithmic flow we are led to a similar conclusion arrived at in the previous 
chapter. The organisation of terms using the expansion around the non-trivial fixed point is 
only suitable if there is some fine-tuning of the parameters, which give an unnaturally small 
value for Co,o, otherwise the suitable expansion is that about the trivial fixed point and the 
distorted wave Bom expansion. 



3.4.3 The Distorted Wave Effective Range Expansion 

In the pure short range case, the terms in the expansion around the non-trivial fixed point 
are in one to one correspondence with the terms in the ERE. The requirement of small C-i 
for the suitability of this expansion could be reformulated as the need for a large scattering 
length, 1/a Aq, as found in nucleon-nucleon systems. 

The generalisation of the ERE is the DW efl'ective range expansion (DWERE) [19, 
20,43,45]. We will see that the energy dependent perturbations around the DWRG non- 
trivial fixed point are in one-to-one correspondence to the terms in the DWERE. To show 
this we insert the solution Vs into the DWLS equation for Ks to obtain an expression for 
the corrective phaseshift, 6s ■ The DWLS equation can be solved by expanding the Green's 
function in terms of a complete set of DWs and iterating to get a geometric series. This 
series can then be summed to obtain 



Notice that the integral in the above expression is equal to AR^J{p, k). Inverting the equa- 
tion and substituting the expression for the ^-matrix element in terms of the phaseshift, 
eqn. (I3.17I ). and the expression for fori? p'^, eqn. (13.241) . gives, 

C{7i)p cot~6s = — Up, k) - \ ] . (3.48) 

Substituting in the expression for Vs and writing everything in terms of the physical vari- 
ables gives the final expression, 

2A. / D^"/^" \ 

Cir])pC0t~6s-Mip,K) = -2/cln -lJ^C2n,rni^^\, (3.49) 

^ ^ n.m ^0 ' 
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where [41], 

dq Ak 



M(p,K) = —M(p,K)= t 

^ Jo 



= 2/fRe 



ZIT] 



(3.50) 



and ifj is the logarithmic derivative of the F-function. This expansion is equivalent to the 
distorted wave effective range expansion (DWERE) first derived by Bethe [19] expanding 
on earlier work by Landau and Smorodinski [47] and more recently examined from an 
EFT viewpoint by Kong and Ravndal [21]. In the expansion all non-analytic behaviour 
has been isolated in the functions C and M allowing an expansion in p and k. The use 
of renormalisation group methods in deriving this equation is a new result and provides an 
interesting insight into the use of not only this expansion but also the distorted wave Bom 
approximation. Beyond that it also provides the full power-counting for this system; more 
than was shown in the work of Kong and Ravndal [21]. We may write the distorted wave 
effective range expansion as, 

C{ri)p cot ~6s - M(p, K) = -^ + l-fcp^ + ..., (3.51) 

ac 2 

allowing definition of a Coulomb-modified scattering length and effective-range. 



1 _ 2AoCo,o 
ac 



rc 



2£m , ,M(ln ^ + C„,a,)) + 2^ g C„,„ (^) , (3.52) 
= .^f^cj^]', (3.53) 



These expansions in aM correspond to nucleon loops in photon exchange diagrams in the 
EFT. 



3.4.4 Proton-Proton Scattering. 

The DWERE was first derived specifically to model low-energy proton-proton scattering 
and proved very successful in doing so [19-21]. The RG analysis tells us that the DWERE 
is only likely to provide a systematic expansion if the parameter Co,o is small allowing organ- 
isation of the terms around the non-trivial fixed point. In the case of neutron-neutron scat- 
tering the same criterion was met because of the large scattering length, «„„ = -18.4fm. In 
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proton-proton scattering the issue is clouded since the Coulomb-modified scattering length 
depends logarithmically upon the the fine structure constant a. If we assume an isospin sym- 
metry for the strong force, then as or ^ 0, the DWERE for proton-proton scattering should 
become the ERE for neutron-neutron (or proton-neutron) scattering length, suggesting that 
eqn. (13.521) can be written as, 

_ = _ + In — + Co,i(ji) ] + .... (3.54) 

Gpp a„n \ 2/u I 

Hence, the criteria for the use of the DWERE to systematically describe proton-proton scat- 
tering is met because Co,o ~ 7r/(2m;ra„„) is small. 

This rather naive analysis must be taken with caution, the isospin symmetry of the 
strong interaction is only approximate. Furthermore, the large separation of scales in the 
nuclear system between Aq and 1 /a tends to amplify the isospin asymmetry of the strong 
interactions due to Coulomb interactions at short scales. For example, the naive argument 
above suggests a„p = a„„. However we find that a„„ = -18.4fm and a„p = -23.7fm in 
the spin-singlet channel, where the 25% difference in the scattering length is a result of 
a much smaller difl'erence in the unsealed eff'ective potential. That said, we can only cite 
studies that have successfully used the Coulomb modified ERE in modelling proton-proton 
scattering [19-21] and the corresponding EFT to model proton-proton fusion [42] 



3.5 Repulsive Inverse-Square Potential 

The next example in this chapter is the repulsive inverse-square potential, 

V,ir) = ^. (3.55) 

This potential is of interest because, firstly, the centrifugal barrier is of this form and the 
analysis will show how the power-counting is constructed for different angular momenta, 
and secondly, because of the potential's relevance to the three body problem [25]. 

The Schrodinger equation is easily solved in the case where /3 > -1/4. The DWs 
satisfying the boundary condition of vanishing amplitude at the origin are 



Mr)= J^Mpr), (3.56) 



Chapter 3. The Distorted Wave Renormalisation Group 



49 



where J viz) is Bessel's function and v = ^jp + 1/4. If < -1/4 then v becomes imaginary 
and the boundary condition at the origin becomes an issue for discussion, this case is con- 
sidered in the next chapter. The long-range phaseshift is given by 6l = n{l/4 - v/2). The 
DWRG equation for the short-range potential is determined by the value at the wavefunction 
close to the origin, 



-t) • 



giving from eqn. (13.211 ) the differential equation for V5, 



dVs M lARy^' V. 



2v+l ^72 
S 



(3.58) 



OA 4;z-r(l + v)2 \ 2 ) - 

The rescaling of V5 is markedly different to the previous examples because of the odd di- 
mension in the A dependence on the right-hand side. The equation is rescaled via the rela- 
tionships, 

Ap = p, Vs(p,A) = ^-^:^^^[-) A^^VsiAp,A), (3.59) 
resulting in the DWRG equation for % , 

dVs dVs ^ 

= + 2^^^ + T^- (3.60) 

oA op 1 - p^ 

The RG analysis of this equation is straightforward but interesting [17,48]. There are two 
fixed points, the trivial fixed point, Vs = 0> and a non-trivial fixed point. If v is an integer the 
latter of these has logarithmic A-dependence. The perturbations around these fixed points 
scale differently to the cases seen so far. In the two previous cases the trivial fixed point 
was stable, with the LO perturbation scaling with A, while the non-trivial fixed point was 
unstable, with the LO perturbation scaling with A ^ In this example, we shall see that 
the trivial and non-trivial fixed points are still stable and unstable respectively (for v 0), 
but that the nature of the stable and unstable perturbations are different, resulting in rather 
different power-counting schemes. 

Perturbations around the trivial fixed point can be determined, as before, by writing 
Vsi^, p) = CA"'^(p), substituting this into the RG equation and linearising. The resulting 
equation for 0, 

^=(cr-2v)<f>, (3.61) 
dp 
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has solutions <p(p) = p"' which upon application of the analyticity boundary condition 
yields the potential, 

V5(p) = gC2„|^j (3.62) 

As promised, for v ^ 0, the fixed point is stable. However, the LO perturbation, instead of 
scaling with A, now scales with A^". The case of v = is interesting, in this case, the LO 
perturbation is marginal, which will be associated with some logarithmic dependence on A. 
The term proportional to p^" is of order d = 2n+2v- 1 in the corresponding power-counting. 
The non-trivial fixed point is determined by solving the DWRG fixed point equation 

P^It^)- 1^ + 7^ = 0, (3.63) 
which is easily solved by the integral, 

J\p) = \ dq^^ = y ' + ^p'^Mip, V), (3.64) 
Jo p2-r 2^ n-v 2 

where, 

cot nv, V i 



M{p,v) = 



(3.65) 



21np, vgN. 

The prime on the sum here indicates that the term with n = v must be omitted when v is 
an integer. The non-analytic term in the solution J{p) can be subtracted off in the case of 
V ^ N as it satisfies the homogeneous fixed point DWRG equation. In the case of v e N the 
logarithmic dependence upon p must be removed in the manner outlined in the Coulomb 
example to give an analytic solution which may be expressed as. 



1 f/ -X ^ - 



= j\p)---p^^M{^,v\. (3.66) 
The perturbations around this fixed point are easily found, 

Vs(p) Vo(p) £^ \Ao/ 
This fixed point is unstable with the number of unstable eigenvectors being determined by 
V. If V lies between the integers - 1 and N then the first perturbations are unstable. If 
V = N then there is also a marginal eigenvector, p^^, which is associated with the logarithmic 
behaviour in the usual manner. 
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The power-counting around the non-trivial fixed point is = 2n - 2y- 1 for a term pro- 
portional to p^". Both of the power-counting schemes are quite different from the Weinberg 
or KSW schemes seen so far. Since the inverse- square potential is scale-free, its strength 
does not provide an expansion parameter in the low-energy EFT, instead it appears in the 
energy power-counting itself. 

Because of the difference in stability in the fixed points from the cases examined 
thus far, the RG flow is quite different. The RG flow is illustrated in the familiar way in 
Figs. (I3.2l3.3ll . Fig. (Q shows the flow in the (Z?o(A), ^2(A)) plane with v = 0.8. The 
trivial fixed point occurs at (0,0) with the non-trivial fixed point at (-l/(2v), 1/(2 - 2v)). 
The solutions flow towards the non-trivial fixed point close to the critical line as in the short 
range case (Fig. 12.11) . As the unstable perturbation becomes important the RG flow peels 
away from the critical line very quickly and into the trivial fixed point. As the strength of 
the inverse square potential increases the rate at which the flow peels away from the critical 
line increases until at v = 1 the flow along that critical line becomes marginal resulting in a 
flow like that illustrated in Fig. 13. 11 For v > 1 the flow on the critical line becomes unstable. 
Fig. l3.3l shows the RG flow for v = 1 .2. In this figure any flow in the region of the non-trivial 
fixed point is pushed away by the unstable perturbations and into the trivial fixed point. 

As the strength of the inverse square potential increases, the instability of the non- 
trivial fixed point 'increases', i.e. the number of unstable perturbations and the order v of 
the unstable perturbation increases. At the same time the stability of the trivial fixed point 
also 'increases'. In order to keep the RG flow in the vicinity of the non-trivial fixed point, 
all coefficients associated with unstable perturbations, Co up to C2N-1 niust be finely tuned. 
For large values of v this fine tuning becomes more and more contrived. 



The expansion around the non-trivial fixed point can be expressed in terms of the 
DWERE, 



However, in order for this expansion to provide a systematic organisation of the terms we 
require fine tuning of all effective couplings associated with unstable perturbations. 

In the case of scattering of a particle with angular momentum / by a short-range poten- 
tial, we have v = Z -I- 1 /2 and there is no non-analytic energy dependence in cot{6 + In/l). 




(3.68) 
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b„(A) 



Figure 3.2: The RG flow for y = 0.8 in the plane (Z7o(A), Z72(A)), where Vs(p,A) = 
Y,n=o ^2n(A)p^". The trivial fixed point is stable, the non-trivial fixed point is unstable with 
a stable NLO perturbation. 
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Figure 3.3: The RG flow for y = 1.2 in the plane (Z7o(A), Z72(A)), where V5(^,A) = 
2^=0 ^2n(A)p^". The trivial fixed point is stable, the non-trivial fixed point is unstable with 
an unstable NLO perturbation. 
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We can identify two expansions, one based on the trivial fixed point 

p-^'-4an(5 + | 
the other on the non-trivial fixed point, 

/'-icot(c5 + |^ 

If the parameters in the latter expansion are natural then it is equivalent to the expansion 
around the trivial fixed point as can be seen by inverting the left and right hand sides. How- 
ever, if Cm-, (n < N) are: finely tuned to a small amount then the expansion obtained by 
inverting the equation will have increasingly large coefficients and will not be a systematic. 
The fine tuning required for this expansion will show itself as shallow bound states or res- 
onances. In nucleon-nucleon scattering there are no shallow bound states or resonances in 
the higher partial waves and the suitable power-counting is that associated with the trivial 
fixed point. 




3.6 "Well-behaved" potentials 

Having looked at two examples in detail, we shall now turn our attention to a more general 
class of potentials, those for which the Jost function exists. That is all potentials for which, 

drr\VL(r)i < oo, [ drr^\VL(r)\ < oo. (3.71) 
Jo 

We shall assume that Vl depends on the low energy scales /f,. 



3.6.1 The Jost function 

We present here a quick overview of Jost's solutions to the Schrodinger equation and the Jost 
function, for a more thorough analysis see Newton [18]. We are interested in two types of 
solutions of the Schrodinger equation, the regular solution (p(p, ku r) and the Jost solutions 
f±{p, Ki, r). These satisfy the boundary conditions: 

<p{p, Ku 0) = 0, <p\p, Ki, 0) = 1 , (3.72) 
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where the prime indicates differentiation with respect to r, and 

f^(p, Ki, r) e-'P' as r ^ oo. (3.73) 

These solutions are interesting as their simple boundary conditions allow us to consider their 
properties as an analytic function of the complex variable p. By writing these solutions as 
power series and assuming the constraints (13.711) upon the potential Vl, we may show that 
ip(p, Ki, r) is an entire function of p^ and that /+(p, Ki, r) {f-{p, Ki, r)) is an analytic function 
of p in the upper (lower) half of the complex plane [18]. 

By analytically continuing f+{p, Ki, r) through the upper half of the complex p-plane 
we arrive at f+{-p, Ki, r), which satisfies the same boundary condition as f-{p, Ki, r). Hence, 
for p > we have, 

f.{p,Ki,r) = U{pe-''\Ki,r). (3.74) 

Since /+ and /_ are linked by analytic continuation we shall write /+(p, Ki, r) = f{p, Ki, r) 
and f-{p, Ki, r) = f(-p, Ki, r). From the reality of Vl and the boundary conditions, it follows 
that ip{p, Ki, r) is real for real p and that 

U{p\Ki,r)r =f{-p,Ki,r). (3.75) 

We introduce the Jost function, T{p, Ki), by writing the regular solution as a superpo- 
sition of the Jost solutions, 

(p{p, Ki, r) = ^(^(-P' Ki)f(P' Kh r) - Tip, Ki)f(-p, Ki, r)). (3.76) 
From the boundary conditions on ip(p, Ki, r) and flux conservation it follows that 

T{p,Ki)=f{p,Ki,Q). 0.11) 

The analytic properties of !F as a complex function of p follow from those of /. i.e T is 
analytic in the upper half of the complex plane. The Jost function also satisfies the conjugate 
relation ( 13.751 ). The Jost function is extremely useful as it provides all the information we 
require to obtain both the long range phaseshift, 6l(p), and the DWRG equation. 

Since the physical wavefunctions must vanish at the origin they must be proportional to 
the regular solution, (p{p, Ki, r). The DWs, ippir), are found by demanding the normalisation: 



^p{^) ~^ s'm(pr + 6{p)) as r ^ oo. 



(3.78) 
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Comparing the asymptotic form for i^p(r) to that for ip{p, ku r) obtained from eqns. (l3.73l3.76l) 
we obtain, 

i^pir) = p (3.79) 

yJr(p,Ki)ri-p,Ki) 

From the asymptotic form, eqn. (13.781) we can also obtain an expression for the S -matrix: 

^2mp) = ZtElJ^, (3.80) 

rip,Kt) 

In the case of the bound states we have the boundary condition of vanishing amplitude 
as r ^ oo. For p positive imaginary and for large r we have, from eqn. (13.761) . 

^(p,Ki,r) ^ l-{r(-p,Kde-^''^'- -np,Kde^''^'-), (3.81) 

which will vanish for large r if Tip, /c,) = 0. This shows that the bound states are given 
by the zeros of the lost function on the positive imaginary axis. The normalisation of the 
bound state solutions is quite tricky and its proof (see Newton [18]) offers no insight so we 
just state the result: 

f/fn(r) = (fjipn, Kj, r), (3.82) 

where the dot signifies differentiation with respect to p. 



3.6.2 The DWRG equation 

In the DWRG equation we require the magnitude of the DWs near to the origin, this follows 
easily from the boundary condition on (p{p,Ki,r). For small r, (p(p,Ki,r) ^ r so that for 
i? <sc we have, 

l^Rf = ■ (3.83) 

The differential equation (13.211) becomes 



dVs{p,K) MR^ 



Vjip^X). (3.84) 



dK In'^TiK Ki)T{-K Ki) p^ - A2 

To rescale this equation we note that the lost function T is dimensionless (this follows from 
the dimensionless boundary condition on the lost solution /), so that the rescaling of the 
potential is exactly the same as that for the Coulomb potential, 

MAR^ 

Vs(p, ki, A) = -^Vs(Ap, Aki, A), (3.85) 
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so that the rescaled DWRG equation is, 



where 



. dVs .dVs dVs ^ .y ^ cm .2 



Cikd = I (3.87) 



which is independent of A since it is dimensionless. 

Once again, we identify the trivial fixed point, Vs = 0- Given the similarity of the 
DWRG equation to those considered in the pure-short range case and also in the Coulomb 
case, the analysis of this fixed point is identical to those examples. It is stable with the LO 
perturbation scaling with A. The power-counting is the Weinberg scheme augmented with 
additional terms in /c, that are easily resolved. The correction to the phaseshift is given by 
the DW Bom expansion, eqn. (13.311) . 

Dividing the DWRG equation through by V| and writing it as a linear PDE in V^^ we 
obtain. 

In parallel to the solutions for the Coulomb and repulsive inverse square potentials our 
starting point for a non-trivial fixed-point solution to this equation is the basic loop integral, 

J(p,kd = f dq^:^C(kJq). (3.89) 

Jo - r 

To isolate the non-analytic behaviour in this integral we need to understand the analytic 
properties of C. Since 

Cik/q) = . (3.90) 

and Tip, Kj) and Ti-p, /c,) are only analytic in the upper and lower half of the complex 
/>-plane respectively, we cannot analytically continue C as function of q into the complex 
plane at all. To circumvent this problem we note that since, 

^'{p, Ki, 0) = lim ^iT{-p, Kdf'ip, Ki, r) - Tip, Ki)f'(-p, Ki, r)) = 1 , (3.91) 

we can write 

r 1 (f'iq,ki,r) f'(-q,ki,r)\ 

C{Kilq) = lim — . • (3.92) 

r-»o 2iq \ T{q, Ki) T{-q, Ki) / 
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This form for C is more promising as it is now written as the difference of two functions, 
one analytic in the upper half of the complex ^-plane and the other in the lower half. If 
the long-range potential does not have an r~' singularity as r ^ then we may define the 
limit\ 

M{Ki,q) = (3.93) 
T{q, Ki) 



SO that, 



C{kdq) = ^ (M{q, ki) - M{-q, kd) , (3.94) 



and the basic loop integral can be written as, 

J(p,kd = :r. dq^:j^M(q,kd. (3.95) 

The integrand in this integral is meromorphic in the upper half of the complex q plane. It has 
poles at the zeroes of the lost function, which correspond to bound states, and propagator 
poles atq = +p. We may ensure analytic properties of the integral by moving the contour of 
integration into the complex plane and not allowing the singularities to 'pinch' the contour. 
We define the contour of integration, C, to run from - 1 to 1 and follow a path in the upper 
half of the complex ^-plane that avoids the point q = and ensures analyticity in p. This 
means that it must also go outside all the bound state poles (see Fig. 13.41 ). to avoid getting 
'pinched' between two of them as ki 0. Hence we write the non-trivial fixed point 
solution as. 



— =-r 

ip,ki) 2/Jc 



dq-^Miq^ki). (3.96) 



Vf\p,ki) 2/Jc ^p^-q' 
It is an analytic function of all scales ki provided the potential does not violate the constraints 
(13.711) or develop an r"' -singularity as each goes to zero. If these conditions are violated 
as some scales go to zero then there will be a need to subtract logarithms using the method 
outlined in the Coulomb example. These logarithmic counterterms will be associated with 
marginal perturbations in the RG. 

The construction of the perturbations around Vf^ follows in much the same way as in 
the Coulomb example. 



1 1 A . /A 



oo / . \ 2«+mi+m2+...-l 



See later for dealing with the r singularity 
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Figure 3.4: The contour for the non-trivial fixed point solution of the DWRG equation. It 
must run from -1 to 1 going outside all bound state poles. 

The LO perturbation is unstable as A ^ 0, the NLO perturbations are marginal and cor- 
respond to the logarithmic counter terms^. All other perturbations are stable as A ^ 0. 
The power-counting scheme in the ^(-independent terms is precisely the KSW scheme. The 
RG flow is the same as described in the Coulomb example and the arguments about the 
usefulness of each power-counting scheme are the same. 

The form for the basic loop integral 13.961 is extremely useful and will form the basis 
for non-trivial solutions in later chapters. The important property it possesses is that it is 
always analytic in the energy p. 

3.6.3 The DWERE and Interpretation of the Fixed Points 

In all our examples, the physical interpretation of the trivial fixed point has, for want of a bet- 
ter word, been trivial. In the pure short range case the non-trivial fixed point corresponded 
to a system with a bound state at exactly threshold. Since we began our analysis of the 

^Not all low-energy scales will have this perturbation as they may satisfy an analyticity in k^j 
boundary condition, rather than simply in kj. 



y tfinr 11 1 1 
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DWRG we have used the non-trivial fixed point as a tool for constructing power-counting 
schemes and for understanding the RG flow but have failed to give an interpretation of it. 
To remedy this situation let us derive a DWERE by using the the non-trivial fixed point in 
this example. We substitute the potential into the equation for the T-matrix, 



M r' 

2^ Jo 



4:t ^ p2 + 



Pi 



The term in square brackets on the right hand side may be written as 



[...] = 



MAR^ 



p^ - Tiq, h) 



p^ -q^ + ie Tiq, ki) 



,p^ iPn 



(3.98) 



, (3.99) 



where 'R{f{z),z Zq} indicates the residue of f{z) at z = zq. This result follows from 
definitions (13.8313.9013.941) and from the result for the normalisation of the bound states, 
(13.821) . Inverting eqn. (13.981 ) and using the expression for the f -matrix in terms of the 
phaseshift correction we obtain after some algebra, 



P ~ I ^ ^ 

— 7J7(COt 6s -l) + M{p, Ki) = — > C2n,mumi, 



A -2«-mi-m2-...„2n,_, mi,. m2 



(3.100) 

where A1(Ap, A^,) = KM{p,ki), we have absorbed any logarithmic counterterms into the 
marginal couplings and we have used the result: 

q f(q,0) 1 



2A 

Mip, Ki) = — 

n 



2i L -q' + ie Tiq, k) vf\p, k) 



q f'{q,0) 
p^ - q^ T{q, ki) 



,p ^iPn 



(3.101) 



which follows from evaluating what is now a closed contour integral using Cauchy's the- 
orem, cancelling the bound state residues and evaluating the remaining residue dX q = +p 
(See Fig. 13.51) . The final result is expressed as. 



\T{p,Ki)rp{coi~5s-i) + 



f'(.P,'<i,0) 11. 2^ 
Tip, Ki) a 2 



(3.102) 



where the distorted wave scattering length and effective range may contain logarithmic de- 
pendence on some of the scales a-, . Relationship (13.941) shows that the imaginary parts on 
the LHS cancel to ensure a real phaseshift, 6s ■ 
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Figure 3.5: Evaluation of the removed non-analytic term in the general DWERE 

The DWERE (13.1021) has been derived before by van Haeringen and Kok [43], fol- 
lowing the work of many others (see their references), in a somewhat more mathematically 
inspired manner. The advantage of this derivation is that, because of the EFT philosophy 
behind the derivation, we gain a physical rather than mathematical interpretation of how and 
when it may be used, i.e. when the series is likely to converge well. 

To gain insight into the non-trivial fixed point let us examine what the DWERE looks 
like in the limit of all perturbations around the non-trivial fixed point going to zero. The 
DWERE when the solution is taken to be the non-trivial fixed point with no perturbations 
(i.e. a ^ oo, fg ^ etc) can be written as: 

cot §5 = --np, Ki)f'{-p, Ki, 0) + /. (3.103) 

P 

This may be re-written as an equation for e^'^* : 

2,g^ ^ T{-p,Kdf'{p,Ki,0)-2ip ^ rip,Kdf'i-p,Kt,0) 
^ T{-p,Kdnp,Ki,0) T{-p,Ki)f'{p,Ki,Oy 

where the second identity follows from eqn. (13.911) . When this is combined with the form 
for e^'^^ , eqn. (13.801) . we obtain an equation for the full S -matrix taking account of both long 
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and short range physics. 



f'(p,Ki,0) 



This equation leads to two interpretations of the non-trivial fixed point, one mathematical 
and the other physical. 

Mathematically, to understand what this equation says let us look at the resulting wave- 
functions, if/p(r). Given the form of the S -matrix we know that at large r the wavefunctions 
must go like 



Mr) ^ —if'(-P, ^y"' - f'ip, Ki, 0)e-'n (3.106) 



1 

2ip' 

and hence are given by (recalling that the short range force acts at zero radius, so that for 
r > the Schrodinger equation is unchanged from its long range counterpart), 

Mr) = ^{f'{-p,Ki,Qi)f{p,Ki,r)-f'{p,Ki,Qi)f{-p,Ki,r)). (3.107) 
2ip 

As r ^ these solutions satisfy the Neumann boundary condition, i}j'{p, /c,, r) = 0, rather 
than the usual Dirichlet boundary condition. This off"ers us a new interpretation of the 
trivial and non-trivial fixed points in terms of boundary conditions on the DWs. The trivial 
fixed point corresponds to the normal boundary condition of vanishing wavefunctions at 
r = 0, while the non-trivial fixed point corresponds to the boundary condition of vanishing 
derivative at r = 0. 

The two solutions to the Schrodinger equation with these two different boundary con- 
ditions are linearly independent. Hence, the non-trivial fixed point changes each DW into a 
solution linearly independent to it and, in eff'ect, changes it 'by the maximum amount pos- 
sible'. This interpretation of the non-trivial fixed point is also appropriate (as of course it 
should be) for the example in chapter 2, in which there was no long range potential, and had 
been implicitly noted by van Kolck [10]. 

Physically, we can understand the non-trivial fixed point in terms of a bound state at 
zero energy. As p ^ 0, f'(p,Ki,0) —> ip,^ so that e^^ —> -I and 6(p = 0) = (n + l/2)n, 
which is the condition for a zero energy bound state. This interpretation is also in parallel 
with the pure short-range case. 



^This is only true for potentials without an r ^ singularity and follows from a power series solu- 
tion in r of the Schrodinger equation. 



_19! 



Chapter 3. The Distorted Wave Renormalisation Group 



62 



3.6.4 Identification of Scales and the Yukawa Potential 

One particular moot point where the Jost function DWRG analysis is important is in the 
inclusion of pions in an EFT for nucleons. We will now look at the issues surrounding the 
Yukawa potential and pions in an EFT for nucleons. 

However, before we can discuss the physics of this problem, we must quickly tidy up 
one final mathematical issue. Our definition of M, eqn. (13.931) . relies upon the existence of 
the limit of the derivative of the Jost function as r ^ 0. This limit only exists for potentials 
that do not have a divergence as r ^ 0. For potentials that do diverge, such as the 
Yukawa potential, we must modify the definition of M. To understand the problem let us 
consider the expansion, in powers of r, of the Jost solution for the Yukawa potential for one 
pion exchange between nucleons of mass M. This potential is defined by, 

MVdr) = 1k„ , (3.108) 

r 

where, 

^- = %r7r' (3.109) 

the inverse 'pionic Bohr radius'^ and m;r, the pion mass are the low energy scales associated 
with the potential. The Jost solution for the Yukawa potential is given by 

K„, m„)r + 2K„r\nm„r + 0(r^)^ . (3.1 10) 

In this case the term ~ r In m^r means that the derivative of the Jost solution is not defined at 
r = 0. However, inserting eqn. (13.1101 ) into eqn. ( 13.911) shows that the identity (13.941 ) holds 
with yVtreg replacing M. The rest of the arguments for the non-trivial fixed point solution 
hold with the substitution M replaced by M"^^ which we may write as 



M"(q, m^) = lim 



— ZKjr In nij^r 



(3.111) 



Without going into the details this results in a non-trivial fixed-point, , defined by 
eqn. (13.961 ) with logarithmic dependence on rhj^, which can be removed in the usual way. 
See refs. [9,44,45] for more details. 



'^f„ = 93MeV is the pion decay constant, gA = 1-26 is the axial coupling of the nucleon. The 
treatment of this scale k,^ as a low energy scale in an EFT for nucleons with explicit pions is still 
open to debate, see below. 
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There has been some debate in the Uterature as to how to handle OPE in an EFT. One 
scheme (WvK) proposed by Weinberg [8] and further developed by van Kolck [40, 49], 
iterates it to all orders. Another, proposed by KSW [50], treats the force perturbatively. 

So far in this section it has been understood that the long range potential has only low 
energy scales associated with it. In this case the analysis is simple and results in a choice 
of two expansions. A Weinberg scheme based upon the trivial fixed point and with terms in 
one-to-one correspondence with the terms in the DW Bom expansion and a KSW scheme 
based upon the non-trivial fixed point and with terms in one-to-one correspondence with the 
terms in the DWERE. We shall see that it is the identification of the low-energy scales in 
the Yukawa potential that has led to such confusion. 

Since the DWRG analysis promotes all low-energy physics to a fixed-point, if it was 
applied to the Yukawa potential with the two identified low energy scales Kj^ and m„ it would 
be equivalent to the WvK scheme, summing the effects of one pion exchange to all orders 
and leaving a choice of power-counting schemes. For the strongly interacting nucleon sys- 
tem we would expect the non-trivial power-counting scheme and the DWERE to provide a 
parameterisation that converges up to the mass of the p-meson, nip = 770MeV [45]. This 
choice of low energy scales has resulted in some moderate success [40,49]. 

Unfortunately there is an issue of debate. In the power-counting around the non-trivial 
fixed point, there are terms like 



The first of these terms is proportional to ml" and occurs at order d = 2n -2m the power- 
counting. The second of these terms occurs at order d = n - linthe power-counting, yet 
because k„ «: ml it is also proportional to ml". Hence, terms of the same order in m„ occur 
at diff'erent orders in this power-counting and the direct link to ChPT has been lost [50, 51]. 

If we wish the nuclear EFT including OPE to be consistent with ChPT, then one solu- 
tion is the KSW scheme [50]. In this case we treat the scale. 




(3.112) 



(3.113) 



= 300MeV, 



(3.114) 
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as a high energy scale so that the Yukawa potential is given by 




(3.115) 



and occurs at the order in the EFT and so does not require summing to all orders. 

From the DWRG point of view, the effects of the Yukawa potential vanish as the cut-off 
A —> and so do not get promoted into the fixed point. Thus the fixed points are the simple 
ones of the pure short range case. In this case the rescaled DW Green's function Gl can 
be expanded in powers of A/A^n and treated as perturbations around the non-trivial fixed 
point. 

Unfortunately, despite consistency with ChPT, the resulting expansion turns out to be 
slowly convergent [37,38,45,46] because of the small separation of scales between the pion 
mass, = 140MeV and A^viv- 

Although the DWRG analysis gives the possible power-counting schemes that can be 
obtained from the non-trivial or trivial fixed point they can only ever be as good as the 
information put in. One must consider very carefully the scales in any particular problem. 



In this chapter we have introduced the DWRG equation and solved it for several examples. 
The analysis for the Coulomb and general "well-behaved" potentials yielded known results 
that have both recently been used in EFTs. Despite these results being well-known, the 
method of derivation is novel and in itself is interesting. The conclusions regarding the 
usefulness of the two possible power-counting schemes are very much the same as arrived 
out for the pure short-range potential. 

The method of solution of the DWRG fixed point equation in the final general example 
will prove extremely useful in later chapters and provides a robust way of solving these 
types of equation. 

The other example considered in this chapter was the repulsive inverse square potential. 
As far as we are aware the results derived here are new. The power-counting schemes that 
are derived are novel and depend upon the strength of the potential. The conclusions about 
the use of either power-counting scheme are also different. If the inverse square potential 



3.7 Summary 
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is 'strong', the use of the non-trivial fixed point scheme may require the tuning of several 

parameters to unnaturally small values, something that may seem a little contrived. It is 
likely that in such a system the trivial fixed point should provide the correct power-counting 
scheme. 



Chapter 4 

The DWRG for Singular Inverse Square 
Potential 

4.1 Introduction 

The inverse square potential, eqn. (13.551) is an interesting example. When /3 > -1/4 the po- 
tential has a well-defined set of DWs. However, if yS < -1/4 the potential becomes singular 
with no well-defined set of DWs. The reason for this is that any two linearly independent 
solutions of the Schrodinger equation cannot be resolved by a boundary condition at the ori- 
gin. Indeed no solution of the Schrodinger equation has a well-defined value at the origin. 
Furthermore, it may appear at first glance that there is a continuum of possible bound states. 
The solution, 

i/roc ^f?KU ^-MEr), (4.1) 

where v = sJ-l/A- - j3 and K^^iz) is the modified Bessel function of the third kind, satisfies 
vanishing boundary conditions at infinity and the Schrodinger equation, yet the lack of a 
boundary condition at the origin means these bound states are not resolved into a discrete 
spectrum. 

Mathematically, the problem with the potential is that the Hamiltonian is no longer self- 
adjoint [27,28]. The solution is to make it so by introducing an extra boundary condition 
to uniquely define a solution [28,29]. This extra boundary condition can be introduced in 
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several ways, either by requiring vanishing wavefunctions at some radius, Rq, or by defin- 
ing some bound state, po [29]. However the boundary condition is defined it necessarily 
introduces a new scale into the problem. Mathematically, this is equivalent to forming a 
self -adjoint extension of the Hamiltonian [27,28]. 

Once the self-adjoint extension is formed we have a complete set of DWs including a 
discrete spectrum of bound states, however there are still problems to be addressed. The 
resulting discrete spectrum consists of an infinite number of geometrically spaced bound 
states. These bound states accumulate at zero energy and have no ground state. The lack of 
a ground state shows that even after the self-adjoint extension is formed the full resolution 
of the short-range physics has not been achieved. Without a ground state the system would 
'implode' radiating an infinite amount of energy. 

The resolution of the attractive inverse square singularity in an EFT is critical to the 
understanding of the three-body KSW EFT [22, 24, 25, 30, 52, 65]. Several papers have 
attempted to resolve the problem by replacing the singularity with an 'effective' potential 
well at short range and then matching the logarithmic derivative of the wavefunctions in the 
two regions [30-32]. In each case the range, R, of the eff'ective potential has been assumed 
to be far less than all inverse momenta. The matching criteria acts, in effect, as the boundary 
condition required to form the self-adjoint extension and does not, as correctly pointed out 
by Bawin and Coon [30] but not by Camblong and Ordonez [31] give a ground state to the 
system. Further to this, using this method, the running of the 'strength' of the effective well 
is multi-valued [32]. 

The DWRG method provides the explanation to these problems. We find that the choice 
of a self-adjoint extension is equivalent to the LO EFT, which in turn corresponds to a 
marginal perturbation in the RG. Because the LO EFT is marginal it cannot distinguish 
between high and low energy states, which is why the self-adjoint extension cannot provide 
a ground state. 

We shall also see that because of the lack of a ground state in the DWs, the DW Green's 
function must be truncated in both the continuum and the bound states to obtain an unique 
solution to the DWRG. The multi-valuedness of the effective potential can then be attributed 
to different methods of truncating the bound states. 
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4.2 Defining the Distorted Waves 

The distorted waves of the inverse square potential satisfy the Schrodinger equation, 



^Mr) - + P^Mf) = 0' (4-2) 

of which the general solution can be written as a superposition of Bessel functions of order 
+ a/1/4 When /3 < -1/4 the orders of the Bessel functions are imaginary and the 
general solution is, 

^p(r) = ^^r{AMpr) +A2J-w{pr)), (4.3) 



where y = ^-p - 1 /4. It is impossible to find a unique solution by defining a boundary 
condition at the origin. To resolve the ambiguity and form a self-adjoint extension we begin 
by considering the lost solution /(/?; r) that satisfy the Schrodinger equation and has the 
unambiguous asymptotic boundary condition, /(p; r) e''"' as r ^ oo. / is easily found to 
be 

mr)= ^e-'^^'^-^^H^ipr), (4.4) 

where hII\z) is the Hankel function of the first kind. We would now like to define a lost 
function, T in terms of the value of the Jost solution, /, at the origin but this limit does not 
exist. Explicitly, in the limit of r «; \/p we have, 

^^/iP\ r) - sin (v log ^ - - /y ) , (4.5) 

where 6 = arg{r(l + iv)}. The solution is to define the Jost function in terms of the value of 
the Jost solution at the point r = Rq rather than at r = 0, resulting in, 

r\ —is. 

2e 4 / 7rv\ 

Tip) = . . , ^ ^ sin r]ip) + (4.6) 

where 

r](p) = -vln— + e = -vln—. (4.7) 

2 po 

Eqn. (14.71) defines the scale po, which as we shall see is simply related to the binding energies 
of the bound states. The definition of the Jost function introduced here is equivalent to 
applying a boundary condition of vanishing wavefunctions at the small distance r = Rq. 
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With these definitions, the free DWs are given precisely as in eqn. (I3.79I) . from which 
follows the usual relation between the S -matrix and the lost function, 

^mdP) ^ ^(-P) = . sHrj(p)-inv/2) 
Tip) sm(r]ip) + inv/2y 

With the lost function we may begin to investigate the bound states. The eigenvalues, -pi, 
of the Hamiltonian are given by the poles of the S -matrix or equivalently the zeros of the 
Jost function. Using the definition of Tip) and ri{p) we find that these zeroes are defined 

by, 

— j neZ. (4.9) 

Thus, these bound states, i^„{r), form an infinite tower of states with geometrically spaced 
energies. It is clear that a ground state does not exist and that the states get infinitesimally 
close to zero energy. Since the distorted waves now form a complete set, the familiar spectral 
decomposition of the Green's function follows. The Green's function with standing wave 
boundary conditions is given by, 

M Ur)Un , M ^ ^„(r)^„(rO 



27r2 Jo p'^ - An p^ + pi 

Explicitly the DWs, given by eqn. (13.791) . are 

M^) = J^ o-i w r L- /oJ ^'^^'^-^'vCpr) - e-'''(^V^,,(pr)), (4. 1 1) 
y 2 2i\ sm{T]{p) + inv/2)\ ' 



and the bound states are. 



V2r sinh (nv) 
—PnKi^Pnr). (4.12) 
nv 

The definition of the Jost function in terms of the Jost solution at r = i?o provides the 
additional boundary condition, which in turn forms a self-adjoint extension of the Hamil- 
tonian. However, in obtaining the Jost function in this way we have implicitly assumed 
Rq Ijk for all k, i.e. Rq is some infinitesimal distance, which makes it difficult to un- 
derstand the exact nature of the boundary condition. Fortunately, this method lends itself to 
different interpretations. Firstly, and perhaps most simply, Ro fixes the phase of the trig-log 
behaviour close to the origin and hence defines a regular and irregular solution. Secondly, 
Pa acts as a choice of a particular bound state, which forces all other bound states to fall into 
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the exponential tower defined by eqn. (14.91) . Importantly this means that the relationship be- 
tween po and choice of self-extension is not one-to-one. Any choice given by p'^ = poe""^'^ 
is equivalent to the choice po. 

4.3 The DWRG equation 

We now move on to the construction and solution of the DWRG equation for the short range 
force. Now that we have arrived at an equivalent of the lost function for the attractive in- 
verse square potential it may be hoped that the general analysis used in section 3.6 will yield 
a method of solving the DWRG equation. However, we are immediately faced with a prob- 
lem. In the previous analysis the contour of integration in eqn. (13.951) for the DWRG basic 
loop integral solution had to go outside all the bound state poles to avoid being 'pinched' 
as the scales associated with the long range potential went to zero, in this case it is quite 
impossible for the contour to go outside all of the bound state poles because there is no 
ground state. 

With a little thought, our inability to define a contour that goes outside the bound state 
poles is not a problem. The position of the poles is controlled by the self-adjoint extension 
defining scale, Pq. As long as we are not concerned with analyticity of the short-range force 
as — ^ then there is no worry with being pinched as pq 0. 

With this problem resolved a new one becomes apparent, namely non-unique solutions 
to the DWRG equation. To explain, if we are to define our basic loop integral solution in a 
method like that outlined in section 3.6 our contour of integration must cross the imaginary 
axis between two particular poles, with no apparent correct choice of which two poles. 
Without resolving this issue, any resulting physical observables will invariably depend upon 
the choice of contour. 

It is clear that this problem is a result of the lack of a ground state which in turn is a 
result of the singular behaviour of the long-range potential close to the origin. The philos- 
ophy behind the DWRG is to 'replace' all the interactions between high energy distorted 
waves with an effective short-range interaction, it therefore makes sense that the deeply 
bound states, which are a result of the inverse square singularity, should be cut-off. Thus, 
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we define the truncated Green's function to be, 

GL(p,A;r,r) = + d^ — — + y' /i — 2^ 2 ' (4.13) 

in which all distorted waves with energies outside the range (-A^/M, A^/M) are removed. 
As we shall see the truncation of bound states and subsequent renormalisation specifies a 
non -trivial fixed point solution and leads to well-defined results for physical observables. 
The method of truncating the bound states chosen here is far from unique, other truncation 
methods will specify diff"erent solutions, however because of the renormalisation procedure, 
physical observables will be independent of the truncation method and of the specific loop 
integral. 

The difl'erential equation, (13.211) . for V5 is now given by. 



dV^ ^ _M_ 



2 + r>^ 



(4.14) 



_ ^ n=—oo ^ ^ " 

where the cutting off of bound states has resulted in a series of discontinuities expressed in 

terms of a sum of delta functions. Since R is taken as very small, the distorted waves go to, 

^ AB(R)smhnv 2 pi 
cosh nv - cos 2?7(A) v 

where 

5(i?) = ^sin^(vln(^)-0). (4.16) 

The final step to obtain the DWRG equation is to rescale. In this problem, the only low 
energy scale is the on-shell momentum, p, which is rescaled to j5 = p/A. V5 is rescaled by 
the relation, 

MB{R) 

Vsip,A) = -^Vs(pA,A), (4.17) 



finally resulting in the DWRG equation. 



sinh nv 



^dA\Vsj ^dp\vj (cosh;rv-cos2vln(A/;?o))(l -^^) 

00 - 

+ - y T7— rA5(A-p«), (4.18) 
V ^ p^ + I 



n=~co 



where we have divided through by to obtain a linear equation in 1/ V5 . Unlike all the RG 
equations seen so far this one has A dependence in the inhomogeneous term on the RHS. 
This dependence is a result of the scale po, which may not be treated as a low energy scale. 
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4.4 Solving the DWRG equation 

The DWRG equation, as always, yields a trivial fixed point, % = 0- The perturbations, 
Vs = CA^(py(p), satisfy the equation 

<P',(p) = vMp)- (4.19) 

We find (f>v(p) = P^- The corresponding RG eigenvalues are v = 0, 2, 4, . . .. Hence, the RG 
solution in the region of the trivial fixed point is given by, 

,2« 



Vs(A,p)= yC2n\^\ r"- (4.20) 



The LO term is marginal, i.e. it does not scale with any power of A and is expected to 
be associated with some logarithmic behaviour in A. However, before we can give a full 
explanation of this we must obtain the full solution to the RG equation. 

It is clear, because of the A-dependence on the right-hand side of equation (14.181) . that 
no other fixed point solution can be found. However, since that dependence is only loga- 
rithmic we may hope to find a slowly evolving solution. (Such a solution will be no worse 
than that encountered when considering the Coulomb potential, in which logarithmic depen- 
dence on A was introduced to remove logarithmic dependence upon k.) A logarithmically 
evolving fixed point solution may still be used to construct a power-counting by perturbing 
about it in powers of A. 

Since (c.f eqn. (IX^ ). 



sinh nv 1 



cosh;7rv - cos 2vln A/po 2i 



, , A inv\ ( , A inv 
cot V In — - cot V In — -I- — — 



(4.21) 



Po 2 / \ pq 2 
we can immediately write down a solution to the continuous part of the DWRG equation in 
parallel to eqn. (13.951 ) considered in section 3.6: 

7ko)7TTT = T^ dq^T-^ cotivln^ , (4.22) 

y^"Xp,A) 2iJc P^-q^ \ Pq 2/ 

which is clearly analytic as p ^ 0. The contour C is shown in Fig. 14.11 it follows a path 
from q = -lio q = I. The position of bound state poles poles shown in Fig. 14. H are simply 
given by, 

ip^ = iEl = iEle'^^l^. (4.23) 
A A 

A A n_l • ^1 T-^n TT\ ^ ■ 
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A 



Figure 4. 1 : The contour in the complex ^-plane for the non-trivial fixed point solution of the 
DWRG equation. It runs from -1 to 1 and crosses the imaginary axis dX q = i. The bound 
state poles occur aXq = ip„ = ipoe""!^ IK. 



The position of these poles varies with A. This is in contrast to the examples considered so 
far where the rescaling of the scales in the long-range potential meant that the bound state 
poles in the complex ^-plane were independent of A (which of course is what one would 
hope if looking for a fixed point solution). The poles in this example move up the imaginary 
q axis as A ^ because the scale that initially controlled their positions, po, is treated as a 
high energy scale and not absorbed into the fixed-point solution. 

Fortunately this migration of the bound state poles with A is to our benefit as it enables 
us to create the necessary discontinuities in the DWRG solution by allowing the bound state 
poles to cross the contour of integration. Dividing eqn. (14.181) through by A and integrating 
it with respect to A from p„ - e to /»„ -I- 6 we get the discontinuity at A = p„, 

A=p„+6 



1 



n 1 



(4.24) 



If we chose the contour C so that it crosses the imaginary axis at the point q = / then from 
eqn. (14.231) it follows that the bound state pole dXq = ipn crosses the contour of integration 



A A n _ i_ 
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when A = Pf, and produces the correct discontinuity, namely 



1 



A=p„+e 



A=p„-6 



Hence, is a "fixed point" solution to the DWRG. From its definition it is readily ob- 



served that V: is invariant under the transformation A he™'^ so that V: has log- 



5 

periodic behaviour in A. Rather than use the rather stretched notation of a fixed point we 
shall, more appropriately refer to it as a limit cycle solution. 

The explicit demonstration of limit-cycle behaviour in the RG for this problem con- 
firms what has been shown at LO in several articles [22, 24, 32]. The first example of a RG 
limit cycle solution was given by Wilson and Glazek [53] in a toy model. Since then, be- 
cause of its importance in the three-body EFT, many authors have taken an interest in such 
solutions [30, 32, 54, 55, 65]. 

Perturbations around this limit cycle are obtained in the usual way. In fact the eigen- 
value equation for the perturbation is exactly the same as that for the trivial fixed point, 
(14.191 ) with the same result. Because the DWRG equation (14.1811 for is linear, these 
perturbations give an exact solution to the RG equation, 

1 1 - / A 



The LO term is marginal and is associated with the logarithmic flow in the limit cycle solu- 
tion. Since there is a marginal (or A-independent) perturbation, it cannot be unambiguously 
removed from any limit cycle solution. This means that there is, in fact, a family of limit 
cycle solutions parameterised by the value of the marginal perturbation Cq. 

Previous examples of logarithmic behaviour in the scale A resulted from a need to 
remove logarithmic behaviour in some low-energy scale and necessitated the introduction 
of an arbitrary scale, /i. In this example however, no new scale needs to be introduced since 
the scale po fills its purpose. 

Where an arbitrary scale was introduced previously it could be traded off against the 
marginal perturbation leaving one remaining degree of freedom. It is not clear that po and 
Co represent the same degree of freedom, resolution of this issue must wait until we consider 
the effects of the long and short range potentials in tandem as po also effects the long-range 
physics. 



A A n_l ■__ ^1 T-^n TT\ ^ 
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l.OxlO'' i.oxio' i.oxio" i.oxio' i.oxio'* 



Figure 4.2: The running of 1 /Vs{0, A) with A for v = 0.3. The dotted lines show disconti- 
nuities. 

The RG flow is illustrated in Figs. 14.21 and 14.31 These diagrams show the flow of 
solutions Vs with A. Fig. 14.21 shows the running of Vf\p 0, A)"'. The bold line 
shows the solution 1/V^ , the dotted lines show the discontinuities. The solution follows 
a logarithmically evolving cycle with discontinuities as each bound state is cut-off". Any 
member of the family of limit cycle solutions can be obtained by adding a marginal term, 
Cq. In Fig. 14.21 this would simply appear as a constant shift up or down without change 
in form or period. From this diagram alone it is clear that Co and po cannot represent the 
same degree of freedom in the limit-cycle solution because variation of po will produce a 
horizontal rather than vertical shift in the form of l/V^. 

The limit-cycle behaviour is illustrated in fig. 14.31 which shows the RG flow in the 
familiar slice through the space (Z?o(A), Z?2(A), . . .). Again the solution is shown in bold. 
The dashed lines show general solutions generated by perturbing with C2. The solution 
acts as a limit-cycle solution that loops at a logarithmic rate in A with a discontinuity as each 
bound state is cut-oflr . All more general solutions move in a cycle (with discontinuities) that 
tends to vf like A^ as A ^ 0. 

Having obtained a general solution to the DWRG we may now understand the nature of 



A A n _ 1 ' 
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Figure 4.3: The flow of VsiP-, A) = 2 bin{ts)p'^^ with A of (bold line) and stable pertur- 
bations about that (dashed lines). The dotted lines show discontinuities. 

the trivial fixed point. As Co ^ oo in eqn. (I4.26I) . the general solution 0, i.e the limit 

cycle solution converges on the trivial fixed point. This means that the trivial fixed point is 
a special case of the more general family of limit-cycle solutions, in which the whole limit- 
cycle has been compacted into a point. If the RG flow is perturbed away from the trivial 
fixed point with a marginal term it moves into a nearby limit cycle. 

The nature of the renormalisation group flow means that no particular tuning of the 
parameter Co is required. Its exact value merely specifies which limit cycle the RG flow 
will converge upon. As we shall see each of these diff"erent limit cycles corresponds to a 
difi'erent 'shift' in the self-adjoint phase r] and that there is a one-to-one relationship between 
all possible limit-cycles and all possible self-adjoint extensions. 

4.5 Full S -matrix and bound states 

Using very much the same procedure that yields the distorted wave eff'ective range expan- 
sion in section 3.6 we may write an expression for the correction to the phaseshift 6$ as, 
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2 CO I %2n 

sinh nv cot 6s = sin 2t]{p) (cosh nv - cos 2t](p)) } C^n — ■ (4.27) 

Equation (14.271) still leaves the relation between the marginal perturbation and the scale 
Po clouded. Since p^ also appears in the equation for the long-range component of the phase- 
shift, 5l, the relationship between p^ and Co may become more transparent in consideration 
of the full S -matrix. Study of the full S -matrix may also highlight the significance of the 
RG solution and make examination of the bound states possible. The full 5 -matrix is 
obtained from equations (14.81) and (14. 271) . 

2iS .2 ^P) 



e = I- 



Zip) 



Z(;?) = cos (;/(;?) + ^)--sin(77(;?) + ^)gC2„|^J . (4.28) 



In the limit of all perturbations, Cin, going to zero (or equivalently the RG solution reducing 
to Vf^) we have, 

/ ITTV \ 

Zip) COS T]ip) + — ]. (4.29) 



2 

Comparing this to equation (14.81 ) for the pure long-range force, we see that there has been a 
maximal change in the arbitrary phase rj. The effect of the solution is to move the choice 
of self-adjoint extension, made in constructing the distorted waves, by a factor n/2. The 
significance of the is now resolved, it is the solution to the RG equation that is 'furthest 
away', in the sense of maximum possible change of physical observables, from the trivial 
solution Vs = 0, which of course would leave the full 5 -matrix unchanged from (14.81) . This 
interpretation is very much in parallel to that arrived at in the previous chapter for the non- 
trivial fixed point for well-behaved potentials. In general we may say that the non-trivial 
fixed point solution changes the boundary condition that defines the regular solution to its 
linearly independent boundary condition. 

Physically one may interpret the limit cycle in terms of its action on the bound 
states, which are shifted to the geometric mean of the original set of states. 

The existence of the marginal eigenvalue, Co, and the scale po leads to an embarrass- 
ment of parameters that may be resolved if it is apparent that they do in fact represent the 
same degree of freedom in the effective theory. This freedom is not manifest in equation 
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(I4.28I) but with the correct rewriting of the parameters, C2„, 

Co = --cotcr, _■> C2J— = -— ^, (4.30) 

2 ^tt sin^ + coscrZ-iq„(^) 



2n' 



in terms of the new parameters, cr and C'^^ , n > 2 we are left with a far more transparent 
parameterisation, 

,2n 



= ,|^, Z(rt = sin + . . f ) + COS . . g Ci, (^] 

(431) 

Recalling that ri{p) = -vln(p/po), eqn. ( I4.31I) for the 5 -matrix clearly shows the overlap- 
ping roles of Co and po- The scale po acted originally to define a self-adjoint extension, 
through the marginal term in the short range force we may change it to any other choice of 
self-adjoint extension. 

The space of all limit-cycle solutions is mapped out as Co varies from -00 to 00 or, from 
eqn. (14.301) . as cr varies from to n. In turn the variation of cr from to tt with po fixed maps 
out all possible self-adjoint extensions. Hence, the relationship between possible limit-cycle 
solutions of the DWRG equation and self-adjoint extensions is one-to-one. This is in stark 
contrast to the relationship between choices of po and self-adjoint extensions where many 
equivalent po's existed. In practice, we may ignore the phase cr and allow the self-adjoint 
extension to be defined by our original choice, po. 

The EFT bound states are given by the zeros of Z{p). The position of the shallower 
bound states are virtually unmodified from the original theory and are controlled by the scale 
Po with minimal dependence upon C2„. There is of course still an accumulation of bound 
states at zero energy, though this is expected as these states exist because of the inverse 
square tail of the potential, which is unaffected by the short-range forces. The deepest 
states within the validity of the EFT, p < Aq, are strongly affected by the higher order 
perturbations. Beyond the range of validity, p > Aq, we can say nothing about the bound 
states except that some other physics must come into play to ensure that a ground state 
exists. 
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4.6 Truncation of Bound states and Uniqueness of Solu- 



Questions still remain about the the truncation of the bound states. It is obvious that some 
sort of regularisation of these states is necessary, yet the method of truncation of the bound 
states chosen above, though convenient and intuitive, is rather arbitrary. However, with a 
little thought it is quite straightforward to show that no matter how the bound state truncation 
is performed, a solution to the RG equation similar to that found above can be found and 
that more importantly the physical results will be identical. 

For example, suppose we choose to truncate the bound states with a second cut-off yu. 
Then the renormalisation prescription states that physical variables should be independent 
of n as well as A. The DWRG equation will be the same as (14.181) without the (5-function 
terms. It can be solved in a similar way to the discontinuous DWRG equation. 



where the contour C(A,/i) still runs from 1 to -1 but now follows a A-dependent path that 
ensures that as A varies no bound state poles cross it. The exact path of the contour (i.e. 
which two bound state poles it passes between) is then uniquely specified by the renormal- 
isation condition on ji. The solution to this DWRG equation is of course quite arbitrary 
without the second renormalisation condition. 

The family of solutions that result are illustrated in Fig l4.41 which shows the previous 
discontinuous solution, 1/V^ , in bold with the discontinuities dotted and the various con- 
tinuous solutions as dashed lines. Which solution is correct depends on jj.. Notice, that the 
correspondence between /i and the continuous RG solutions is not bijective; any two values 
of n between |p„| and |p„+i| will map to the same continuous solution. More generally any 
prescription for truncating bound states will give a unique RG solution, jumping to the next 
branch of the solutions as a bound state is cut-off". 



tions. 




(4.32) 



_1 T T 



_ J? O _ 1 
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Figure 4.4: The discontinuous solution l/V^ and the family of continuous solutions that 
result from an alternative method of truncation of the bound states (y = 0.3). 



4.7 Summary 



In this section we have looked at the attractive inverse square potential and the use of the 
DWRG equation in resolving its singular behaviour. The first step was to introduce a self- 
adjoint extension to the Hamiltonian by defining a lost function. This allows definition of 
the DWs but does not resolve the singular behaviour. In particular, even after definition of a 
self-adjoint extension there is no ground state. 

In order to produce a unique result for the EFT 5 -matrix it was necessary to truncate 
the deeply bound states that are a result of the singular behaviour of the long range potential. 
The basic loop integral solutions to the DWRG take the form of limit cycles towards which 
all more general solutions tend as A ^ 0. 

The limit cycle solutions are associated with a marginal perturbation in the DWRG. 
The degree of freedom embodied in this perturbation is the same as that in the choice of 
self-adjoint extension. Despite the complications in obtaining it, the EFT 5 -matrix, (I4.31I) . 
takes a satisfyingly simple form. 

It is difficult to link the limit cycle behaviour seen in our approach to the those seen by 
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Bawin and Coon [30] and by Beane et al [32]. Their methods rely upon constructing a self- 
adjoint extension and short-range force in one step rather than the two step process that our 
method requires. So whereas their LO force defines the self-adjoint extension ours serves 
to change the original choice. What is particularly pleasing is the common approach of 
this chapter and section 3.6. The limit-cycle solutions do not require any particular special 
treatment apart from the truncation of bound states to specify the 'branch' of the multi- 
valued solution. 

The truncation of bound states in this problem is forced upon us in seeking unique 
physical results. It resolves many of the issues raised about multi-valued solutions to the 
RG in refs. [30, 32]. The need to truncate bound states in this problem occurs because their 
position is controlled by a scale, pq, that cannot be rescaled and absorbed into the limit- 
cycle. 



Chapter 5 

The DWRG for Three-Body Forces 



In this chapter we shall consider the DWRG equation for three body forces (3BDWRG). In 
particular we shall derive the 3BDWRG equation for a general system for which the two- 
body forces are known but the three-body force is to be constructed using the 3BDWRG. 
The equation is applied to a system in which the pairwise forces are "well-behaved" in the 
sense of section 3.6 and then, more interestingly, used to examine the power-counting for 
the three body force in the KSW EFT for short range forces. 



5.1 The Faddeev Equations and the Three-Body Force. 

For a system with three particles the LS equation is inadequate. Although the three-body 
T-matrix satisfies the LS equation, it is not possible to use this equation to find it because 
its kernel is not compact. This is because a single LS equation is not enough to specify all 
the boundary conditions for a three-body DW. The solution is to use the Faddeev equations 
(see Appendix O- 

Consider the case of three identical bosons interacting via a pairwise interaction V2b- 
The wavefunction is written as 

|T^) = |(Ao) + (l+/')r), (5.1) 
where Ii/tq) is the 'in'-state and solves the free Schrodinger equation and P is a permutation 
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operator, which permutes particle indices and has the matrix elements 

<xi2,y3|P|x;2,y'3) = 6\xn - x!,,)6\y, - y\) + 6\xn - x^,)6\y, - y^) 
= 6' |xi2 + ^x;2 + y^ j 6' |y3 - ^x[^ + ^y^ j 

+6' |xi2 + ix;^ - y^j 6' |y3 + ^x',, + iy^j. (5.3) 

The Faddeev equation for the wavefunction is written as, 

\r) = G^,(pr(pMo) + G^o(py(p)p\r), (5.4) 

where is the T -matrix given by the equation, f^ip) = Vie + V2sGq(/7)?"^(/7). The corre- 
sponding equations for the full three-body Green's function and T-matrix now follow from 
similar decompositions. If we write the full T-matrix, T, in terms of components, 

r\p) = {\^P)T\p), (5.5) 

then the Faddeev equation for the T -matrix component is, 

T\p) = t\p) + t\p)Gl{p)PT\p). (5.6) 

It is possible to incorporate three-body forces directly into the Faddeev equations, how- 
ever this is not a suitable method for obtaining the DWRG for the three-body force as we 
wish to separate it from the two-body interactions in order to treat it distinctly. Fortunately 
it is quite possible to use the method outlined in chapter 3 to separate the T-matrix. Writing 

T\p) = T^p) + (1 + nBiP)Gt{p))t;j,{p){l + Gl{p)Tii,{p)), (5.7) 

where 7'2g(p) is the T-matrix in the presence of pair-wise forces only and must be obtained 
using Faddeev-like equations, we obtain the equation for T:^g{p), 

TUp) = V3B + V^B^MnBiP)- (5.8) 

^ We will work in the coordinate system (x,y, y^), given by 

Xij^Vi-rj, yk^rk-^iVi + Vj), (5.2) 

where [i, j, k} is an even permutation of { 1 , 2, 3} and ri , r2, r3 are the absolute coordinates of the three 
particles. 
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Notice that this equation is connected since all three particles must interact by the very 
nature of the three-body force. 

Eqn. (15.81) leads, after introducing standing wave boundary conditions^, to the distorted 
wave renormalisation group equation for the three-body force, 

= -y^B^V.B. (5.9) 
dK OA 

Solving this equation now offers a recipe for constructing the power-counting for the three- 
body force. 



5.2 Full Green's Function and Projection Operator. 

We shall consider the simple case of zero total angular momentum. We may identify three 
different types of wave function characterised by their incoming boundary condition. The 
wavefunctions with three free incident particles, 1^^^), ^ are given by eqns. (15.115.41) . If there 
is a two-body bound state with binding momentum y„, there are wavefunctions, l^py,,) with 
the incoming boundary condition of the bound state and one free particle, the equations that 
define these are found by taking the residues of the equations for 1^^^) atk = ijn, 

I^;^„) = (i + p)|.a;,„), (5.10) 

l'A;,r„) = kp,r„> + G,;r^/'|^;,,), (5.11) 

where the incoming wavefunction \Xp,y„) satisfies the equation, {Ho + V2B + P^/M)\Xp,y„) = 
and has a bound pair. Finally there are the three-particle bound states, |¥„); these satisfy the 
homogeneous version of eqn. (15.1 II) with the boundary condition of vanishing amplitude at 

^Standing wave boundary conditions on the full three-body Green's function are quite tricky 
to visualise. We use the expression to refer to a Green's function which uses a principal value 
prescription on the propagator and avoids any complications in the DWRG arising from using a 
complex valued Green's function. 

^The wavefunctions are labelled by their total centre of mass energy, p^jM, (M is the mass of 
one particle) and the relative momentum of one pair, k. Of course it doesn't matter which pair since 
all three particles are identical, and the wavefunction is symmetric with respect to interchange of any 
of them. 
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infinity. We normalise these states with the conditions, 

(*I'«l*Pm) — Sfjjn, (5.12) 

(^/.rJ^/ArJ = "^Sip - p')6nm, (5.13) 

2 

{'^p,k\'Vp',k') = j6ip-p')6(k-k'). (5.14) 

We may now define the full Green's function, Q2b{p), using the three-body complete- 
ness relation. 



2 2q p^ - + ie 



+ -J dqj dk , I (5.15) 



p^ - q^ + 16 



An ^ p^ + pI 2n^ j_y2 2q p^ - q^ + ie' 



•,2 



where P(q) is a projection operator that projects out all states with energy q'^/M and is given 

by, 

nq) = X KyJ<%J if - To < < 0, (5.17) 

yn>\q\ 

nq) = 2 Ky,)<Ky,} + " f ^^I^^Xn^l if > (5-18) 

The Green's function with standing wave boundary conditions is obtained by applying a 
principal value prescription to the pole at p = q. The cut-off to obtain the DWRG may be 
applied to eqn. (15.161) by changing the upper limit to A^, 

g2B(p;A) = — > 2^ 2 2- (^-i^) 

4;/r ^ p-^ + pf^ 2n^ j_yi 2q p^ - q^ 

We have assumed that the bound states do not need truncating (which as we shall see is not 
the case in the KSW EFT). The resulting equation for V^b from eqn. (15.91) is 

= TT^35-7^V3B. (5.20) 
OA 2n^ - p^ 

In the 3BDWRG equation the projection operator replaces the distorted waves. 
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The three-body effective interaction, which acts at some small hyperradius R and some 
arbitrary hyperangle or is given by 

CVpAViBip, Ki, A)|¥p,,^> = R V^Bip, Ki, k, k'; A) ^;,iR, a)^p,AR^ a), (5.23) 

where the scales /c, are low-energy scales associated with the two-body long-range potential. 
We shall assume, as in the previous two chapters, that the off-shell matrix elements of 
the potential are not required to describe on-shell scattering. However, due to the extra 
degrees of freedom in the three body system, the potential must in general depend upon the 
unconstrained momentum k and k'. To obtain the matrix elements of the potential between 
states with incoming and/or outgoing two body bound states we simply set k and/or k' to the 
appropriate iy„. The boundary condition for Vi,b is the analyticity condition, that is it has an 
expansion in even powers of p,k and k' and an expansion in positive (possibly even) powers 

of Ki. 

Substituting the form for the potential into eqns. (15. 18115.20b we obtain, 

d 

— Vwip,Ki;k,k;K) = 
aA 



MR 1 



^ ViBip, Ki', k, iy„; N)Vw{p, Ki; iy,,, k'; A)|*PA,^„(i?, a)\^ 

' . (5.24) 



2;r2 A2 - 

+- r dk'y^B{p,Ki\k,k'';K)V,B{p,Ki;k'\k';K)\^>!,riR,af 

Jo 

The last step to obtaining the 3BDWRG equation is to rescale all quantities in terms of A, a 



process that of course depends upon the explicit form of the distorted waves. 



"^The s-wave wavefunctions depend only on Xij = |x,j| and yt = \yk\- The hyperpolar coordinates 
are defined by, 

^ = ^4 + ^yl «^ " ^''''^''''(■^] ■ ^^'^^^ 

R is independent of the labelling of particles, while a is not. Since these are s-wave DWs we have 
factored out the term Xijyic so that they are normalised as, 

RdR da'i'* ^(/?, a)*P^' k'(R, a) = —5{q - q)6{k - k'), (5.22) 
Jo Jo 4 



consistent with eqn. (l5.14t . 
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5.3 Properties of the Projection Operator. 

The 3BDWRG equation, which comes from rescaling eqn. (15.241) is obviously going to be 
very complicated. This complication arises because the equation describes the coupling 
of the three-body force to each of the different distorted waves. Let us briefly summarise 
our strategy for solving it. The study of the trivial fixed point that invariably results in 
these equations remains simple since in linearising the equation about the point we are able 
to ignore the complicated coupling to all the different channels. However, the basic loop 
integral solution which may lead to a non-trivial fixed point or limit cycle solution and to a 
general solution to the DWRG equation will depend on these details. 

Our strategy in finding a basic loop integral solution will be to generalise the approach 
used in section 3.6 and subsequently in chapter 4. If we look for solutions to the 3BDWRG 
equation that are independent of momentum, then these will satisfy a considerably simpler 
equation because the three-body force decouples from each of the DWs allowing us to divide 
through by V^^ to obtain a linear PDE for V^^ . In unsealed coordinates, eqn. (15.241) becomes, 

dV^B{p,Ki-K) MR . ^^ {R,a\P(A)\R,a) 

^ = ^V,,(p,.r,A) . (5.25) 

In order to use the methods of section 3.6 to solve the fixed point equation that results 
from this simplification we need to know about the analytic properties of the projection 
operator, just as we knew the analytic properties of the lost functions. 

Consider the expression for the full Green's function as a function of energy, E, where 
£■ is a general complex number, 

4nZ-j E + E„ 271-2 E-E' 

where, 

PiE)= J^nVME). (5.27) 
V 4-E 

Eqn. (15.261 ) defines the function Qib^E) as an analytic function of E except when E e 
[-yl/M, oo) and when E = -E„. It is a simple matter to see from its definition that ^2b(E) 
also has simple poles at E = -En with residues, 

nQ2B(E), -E„} = "' ■ (5.28) 

An 
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We can also see directly from its definition that Q2b{E) has a branch cut running from -yl/M 
to infinity along the real axis. Using the result, 

, f(x) , fix) . 

ax ;- = T ax inj{xo) (5.29) 

Xo — X + le J_oo Xq — X 

and its complex conjugate we may use eqn. (15.261) to show that, 

P{E) = i7T(02B(E + ie) - g2B(E - ie)). (5.30) 

Hence, for E e [-y^/M, oo), the projection operator P(E) is defined by the discontinuity 
across this branch cut. From eqn. (15.181) that defines the projection operator for real E in 
terms of the distorted waves it follows that the discontinuity across the real E axis must in 
fact be the result of several 'overlapping' branch cuts that run from each two-body threshold 
E = -yl/M to £■ = oo and from the three-body threshold E = Oto E = co. 

Substituting eqns. (15.2815.301 ) into eqn. (15.261) we can see that the spectral decomposi- 
tion of the Green's function is equivalent to a dispersion relation for Q2b{E), 

^ /I7\ , \^ 'f^^&2B{E),-En] 1 C ,„,@2b{E') /c-iix 

where the contour C is illustrated in Fig. 15. 11 This relation is consistent with the conclusion 
that the Green's function is an analytic function of E with no other singularities than those 
observed^. 

When we consider the Green's function as an analytic function of p = ^ME the 
singularity structure found as a function of energy is mapped into the upper-half of the 
complex p-plane. The singularity structure in the lower-half of the complex p-plane is 
difficult to determine, fortunately we do not need to know any analytic properties in that 
region. 

Now we know the analytic properties of ^2b(p), expression (15.301) can be used to obtain 
an analytic continuation of the projection operator, Pip). Using eqns. (15.2715.301) we can 
write 

np) = ^{02b{p) - 02b{-p)). (5.32) 

^ There is the possibility of other symmetric singularities whose contributions cancel. If these 
exist they do not effect the validity of the results, which only rely on the dispersion relation. 



X 
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Figure 5.1: The singularity structure of an example of the full Green's function for complex 
E with three three-body bound states and three two-body bound states. Also shown is the 
contour C used in the dispersion relation. 

Due to the form of the equation we cannot readily continue the projection operator into the 
complex p-plane without assumptions on the analytic properties of Q2b{p) in the lower half 
of the complex plane. 

Eqn. (15.321) is the generalisation (in unsealed coordinates) of eqn. (13.941 ) and provides 
us with the means to solve the 3BDWRG quite generally ^. 



5.4 The Solution of the 3BDWRG for "Well-behaved" Pair- 
wise Forces. 



If the pairwise forces are "well-behaved" in the sense of section 3.6 then the distorted waves 

must vanish as as R ^ 0. We can immediately write down a form for the full Green's 

^The extension of the analysis so far in this chapter to general N body forces is straightforward, 
Because of eqn. ( I5.32t . which is entirely independent of the number of bodies interacting, the method 
used in finding a non-trivial solution to the DWRG equation also generalises. 



C A 



I O _ 1 



15! T-1 • ■ 
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function and the distorted waves at the point at which the three-body force acts R, a, 

\^A,k{R, a)P = A'R'\<fi(afmk/K Ki/A), (5.33) 
|^A,,„(^, at = A'R'\^iatDliKi/A), (5.34) 
{R, m2B(A)\R, a) = MA'R\{afDG{KilA), (5.35) 

where (p{a) is an unimportant functions that describes the hyperangular behaviour of the 
DWs close to the origin. Inserting these forms into eqn. (I5.24I) . the correct rescaling for the 
three-body force is, 

MR^A"^ 

V,b(p, Ki, k, k'-A) = —^V^b{Ap, Aku Ak, A, k'. A). (5.36) 
Resulting in the 3BDWRG equation, 

dVsB .d%B , \^.d%B fd%B , f,dV3B , 
A = p h > Ki h k h k h 4V3B 

dA ^ dp ^ • dki dk dt 

1 



2 Visip, Kf, k, i%; A)V3b(p, Kf, i%, k'\ A)£)'^(^,) 



1 -p^ 

2 ^ ^ 

+- dk''%B{p,h\kk''\A)%B{p,h\k'\k'\A)D3{k'\ki) 

Jo 



(5.37) 



This equation is quite complicated as it has to describe the coupling of the three-body 
force to each of the distorted waves. The trivial fixed point solution, V^b = 0, corresponds 
to a zero three-body force. The perturbations around it, 

li,m,n,n'=0 ^ ^' 

have RG eigenvalues v = 2m + 2n + 2n' -I- 2 + 4 = 4, (5), 6 All the eigenvalues are 

positive so the fixed point is stable. This fixed point is simply Weinberg counting, in which 
the first three-body force term occurs at (Q/Aq)'^. 

To find a non-trivial fixed point solution, we generalise the arguments of section 3.6. 
A non-trivial fixed point that is independent of the momentum variables may be found by 
solving the equation. 



P- 



dV-^ ^dV;l 2m 



+^,-^-41/3-^-— ^(2)g(^,)-£)g(-^,)) = 0, (5.39) 

op OKi I- p^^ ' 
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which comes from using the form for the projection operator (I5.32I) . This is solved by the 
contour integration, 



where C is a contour going from -1 to 1 via the upper-half of the complex ^-plane outside 
all three-body bound state poles and two-body elastic thresholds on the imaginary axis. 
This solution is analytic in p but may have simple logarithmic dependence on some of the 
scales that can be removed using counterterms with logarithmic dependence on A. The 
momentum-independent perturbations around this fixed point solution take the form 



with RG eigenvalues v = 2m + Ij - A = -4, (-3), -2, This fixed point is unstable 

with 2 energy-dependent unstable perturbations and possibly 2 more /Crdependent unstable 
perturbations. Like the power-counting schemes in the two-body higher partial waves it 
seems implausible that this fixed point is ever likely to provide a useful expansion in physical 
systems because of the fine-tuning required in the unstable perturbations. 



5.5 Three Body Force in the KSW EFT for Short Range 



In this section we will return to the EFT considered in the chapter 2 and its extension to 
systems of three particles. When we construct an EFT Lagrangian, we have to include all 
terms that do not violate the observed symmetries of the system. In a system in which all 
energies are far lower than the inverse range of the interactions this leads to a Lagrangian 
containing not only all the four-point vertex interactions (including derivative couplings) 
that correspond to the energy-dependent contact potentials but also to six-point and higher 
vertex interactions needed to model interaction channels that only become available in the 
presence of more than two particles. Of course the six point interactions are only required 
in a system with at least three particles, while the eight point terms require four particles 
and so on. The two-body couplings are entirely determined by two-body observables and 




(5.40) 




(5.41) 



Forces 
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should be taken as given when considering a three-body system. The three-body couplings 
are fixed by three-body observables. 

In reactions between a bound pair and a third particle the pair- wise interactions result in 
a long range force due to the exchange of one particle. The range of this force is comparable 
to the scattering length of the two-body interaction. 

If the two-body system is weakly interacting, the correct power-counting in the two- 
body system is the Weinberg system. The two-body scattering length a2 ~ 1 /Ao, and hence 
provides a short-range exchange interaction. This system does not provide any particular 
difficulty in the three-body system since we can use naive dimensional analysis to deter- 
mine at which order each diagram occurs. In this expansion the three-body force occurs 
at {Ql^of as in the previous example. (There is, of course, the alternative finely-tuned 
three-body force outlined in the previous example.) 

In a KSW system the fine tuning of the LO perturbation around the non-trivial fixed 
point of the RG results in a large scattering length, 1 /a2 ~ <2 1 / Aq, and hence a long- 
range exchange interaction. This causes a significant complication as all the diagrams in- 
volving the leading four-point vertex need to be resummed [22,66]. Since naive dimensional 
analysis is no longer applicable the power-counting for the three-body forces must be deter- 
mined by renormalisation of the diagrams to which they contribute. 

The 3BDWRG allows us to resum the effects of these diagrams and promote them to a 
fixed point. In the KSW EFT we need only resum the effects of the scattering length terms 
since it must be treated as a low-energy scale, it is not necessary to resum the effects of the 
effective range terms since these are perturbative 1 Irj ~ Aq 

The study of this problem has a long history. Skomiakov and Ter-Martirosian (STM) 
[56] were the first to derive equations for a system of a bound pair and particle with two- 
body contact interactions defined solely by a scattering length. There work was followed up 
by Danilov [57] who realised that in some systems the STM equations did not have unique 
solutions and required the input of a single piece of three-body data, namely the binding 
energy of the three-body state. The work of both of these predates that of Faddeev [58] 
in resolving the disconnected nature of the three-body equations. The STM equations rely 
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upon the zero range nature of the two-body interaction to reduce the problem to a one- 
dimensional equation. 

An alternative approach to the problem was given by Efimov [25]. Efimov's equations 
show the singular nature of the system far more transparently than those of Skomiakov and 
Ter-Martirosian. We shall look at these equations below: in a nutshell Efimov's approach 
allows us to express the system in terms of a hyperradial inverse square potential. The 
strength of the potential depends upon the statistics of the system. In the case of three 
Bosons and three nucleons in the i channel the inverse square potential is attractive and 
hence singular, requiring a choice of self-adjoint extension. 

Since the conception of EFTs this system and the related singular inverse square poten- 
tial have received much attention [22-24, 30-34]. In particular, Bedaque et al [22, 24] have 
used extended versions of the STM equations to create a EFT that has proven highly effec- 
tive in describing three-body data. In nuclear systems, the so-called pionless EFT proves 
effective in the s-wave 7 = 3/2 system [66] and in higher angular momentum systems [59] 
without three-body forces. These systems are repulsive at short distances because of the 
centrifugal barrier or the Pauli exclusion principle and are therefore insensitive to short dis- 
tance physics and three-body interactions. In Efimov's approach these systems correspond 
to repulsive hyperangular inverse square potentials and so, from the DWRG analysis in sec- 
tion 3.5, the power-countings correspond to a weak system with LO three-body force scaling 
with ((2/Ao)^'' or a finely-tuned system with a LO three-body force scaling with (Q/Aq)"^'' 
where v is the strength of the hyperangular potential. Given the apparent success achieved 
without three-body forces the first counting is appropriate. These systems are now well 
understood and we shall not concern ourselves with them. 

The systems that correspond to singular hyperangular potentials are still being stud- 
ied. Bedaque et al have incorporated three-body forces into the STM equations and used 
them to resolve the singular behaviour of these system [22, 24]. The effective three-body 
forces contain free parameters that must be fixed by fitting to three-body data. The power- 
counting proposed by Bedaque et al has three-body forces occurring at orders (Q/Aq)^" 
where n = 0, 1, 2, . . .. The LO force is marginal. Each of the three-body force terms have 
cyclic behaviour in the cut-off, A. Phillips and Afnan [33] have also looked at the problem. 
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their solution is to incorporate the 2+1 scattering length into the equations in a manner that 
'fixes' the singularity and ensures that they give the correct 2+1 scattering length. The in- 
put of a single piece of three-body data is surely equivalent to the LO three-body force of 
Bedaque et al and both are equivalent to a particular choice of self-adjoint extension of the 
Hamiltonian. The difference in their stances is in attributing it to a three-body force and the 
need for higher order three-body forces. 

In this section we shall introduce Efimov's equations for three Bosons and explicitly 
demonstrate the connection between the inverse square potential and the KSW EFT three- 
body problem. We shall discuss the case of infinite two-body scattering length, where Efi- 
mov's equations may be solved analytically, and show that in this case the power-counting 
for the three-body force is the same as proposed by Bedaque et al [24]. We shall then use 
the 3BDWRG to discuss the more general case of systems with finite two-body scattering 
length. We will show that the introduction of a finite scattering length does not change the 
power-counting for the three-body force and that the 3BDWRG solution is governed by a 
limit-cycle solution. 



5.5.1 Efimov's Equations 



The two-body interaction may be written as a boundary condition on the wavefunction at 
zero separation^. The full three-body wavefunction, is defined by the wave equation 
and the three boundary conditions. 



(//o - Em = 0, 



dX: 



Xiji^ij, ym 



(5.42) 
(5.43) 



From here we shall drop the particle number indices by writing the wavefunction in terms 
of the permutation operator, P, 

|T) = (1+/')|«A). (5.44) 
Substituting this form into eqn. ( 15.431 ) we obtain. 





d 




— x<x, y|(A) + < 
ox 



x=0 



^Compare this to the observations about Neumann and Dirichlet boundary conditions and their 
connections to the fixed points of the two-body DWRG. 
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We consider the case of zero angular momentum and write the wavefunction as il/{x, y) 
xy{x, y|i/r), so that eqns. (l5.42l5.45l) become, 

^2 o a2 



dx 



+ -Ay^ y = «2V(0'J)• 



(5.46) 
(5.47) 



.v=o y \ - I 

In addition we have the boundary condition which follows from the definition of if/: ip{x, 0) = 
0. Eqns. (15 .4615.471) are the wave equation with non-trivial boundary condition and are, in 
general, extremely difficult to solve. To examine them more closely we rewrite them in 
terms of the hyperpolar coordinates, x = R cos a,y = y/3/2R sin a. 



■7r^ + + —— + ME\ MR, a) = 0, 

dR^ RdR R^da^ IV'v , ; 



di//{R, a) 
da 



8 



n 



n 



\a=n/2 V3 

These provide us with the information we require to study the DWRG in this problem 



(5.48) 
(5.49) 



5.5.2 The Infinite Scattering Length Limit. 

Efimov's equations, (15.4815.491 ) may be solved analytically in the limit of infinite scattering 
length. We shall turn our attention to this case as it is considerably simpler than the more 
general case of finite scattering length. In the limit a2 ^ oo the boundary condition becomes 
separable. Writing, 

i/fiR, a) = F,(R) sin sa, (5.50) 



we obtain. 



sn S . sn ^ 

-5 cos h — - sm — = 0, (5.52) 

2 V3 6 



where p = ^ME. These equations are now a simple matter to solve. The transcen- 
dental equation (15.521) must be solved numerically. It has two imaginary solutions sq = 
006237 . . . and an infinite number of real solutions, the smallest of which is si = 4. 
The imaginary roots of Eqn. (l5.52l) result in an attractive inverse square potential for the 
hyperradial equation ( 15.511) while the real roots give repulsive potentials. 
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The hyperradial part of the wavefunction is given by, 



FsiR;p) = \ 




'"'^P^J-iApR)) if ^ ^ R, 

(5.53) 



where s = \ms. In the imaginary hyperangular momentum channel, since this is an attractive 
channel, we have, as was necessary, defined a self-adjoint extension of the Hamiltonian by 
choosing a phase, ri{p). The bound states of the system only occur in this channel and are 
given, in parallel to the bound states of the attractive inverse square potential, by 



These states, like the attractive inverse square potential, form an infinite tower of bound 
states. The lack of a ground state in three body system acting under such forces was first 
noted by Thomas [26] and is known as the Thomas effect; of course such a state of affairs 
is a result of an inexact treatment of short-range physics and will be resolved by the EFT 
precisely as it was in the chapterlH The accumulation of bound states at zero energy is due 
to the tail of the inverse square potential, and hence a result of this special case of infinite 
scattering length. These states are known as Efimov states [25] and are not a result of any 
mistreatment of physics but are genuine phenomena ^. 

Since the Hamiltonian conserves hyperangular momentum we can expand any three 
body scattering amplitude in terms of the hyperangular eigenfunctions and corresponding 
hyper-phaseshifts. This makes the DWRG analysis particularly simple since the equations 
become one-dimensional. The effect of the three-body force may be analysed in each hy- 
perangular channel, in other words dependence upon the unconstrained external momentum 
terms k and k' , favoured earlier in the chapter, can be dropped in favour of dependence upon 
the hyperangular momentum s. The differential equation for the three body force in each 

^It is hoped that Efimov states may be observed experimentally by tuning a diatomic bound state 
in a Bose Einstein Condensate. [68] 




(5.54) 
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channel is then, 



aV>,A)_M^|TU^.«lV«(,,^). if.eR. (5.55) 



OA 2n^ 



if^^R. (5.56) 



In the repulsive channels (real s) the resulting DWRG equation is precisely that con- 
sidered in section 3.5 and the analysis there can be followed to obtain the power-counting 
for the three-body force in those channels. The effects of the force can be expressed as a 
distorted wave Bom or effective range expansion for the correction to the hyper-phaseshift. 
Following the discussion in section 3.5 the fine-tuning required for a non-trivial fixed point 
power-counting for these forces would be rather contrived. For example, the least repulsive 
channel has the eigenvalue si = 4 and would require three unnaturally small parameters to 
ensure the RG flow remained in the region of the non-trivial fixed point. Hence, the appro- 
priate power-counting for the three-body force in these channels is based upon the trivial 
fixed point and provides only very small (at least of order {Ql Aof) contributions to the 
scattering observables. 

In the attractive channels we must turn to the analysis considered of chapter HJ The 
DWRG analysis is precisely as found there and is very different from that in the repul- 
sive channels and results in a three-body force that contributes at a surprisingly low order, 
(2/Ao)°, to the scattering observables. In this role it is more important than effective range 
corrections which occur at an order, {QIAqY. This power-counting for the three-body force 
agrees with Bedaque et al [24]. 

The implications, for both EFT and potential modelling of the three-body problem, of 
such a prominent three-body force will be examined later. Before then we shall show how 
the DWRG analysis for the three body force needs to be modified after the introduction of 
a finite two-body scattering length. 

The method outlined in this section is useful because everything may be dealt with 
analytically. However, in order to fully grasp the general problem with finite scattering 
length, our general analysis must be based upon the the 3BDWRG. This not only allows us 
to deal with the problem of the different channels but also allows us to presents results for 
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observables in more practical format. 

5.5.3 Finite Scattering Length 

The Distorted Waves and the DWRG Equation 

In reaching eqns. dS.S 115.521) we took the two-body scattering length to infinity. This limit is 
of course also reached if we take the hyperradius to zero, meaning that at small hyperradii 
the general distorted waves have the same functional form (in R) as their infinite scattering 
length counterparts. In particular we must have for any DW, 

WiR, a) = Difia) sin(5o In p,R), (5.57) 

where we must chose the scale in order to define a self-adjoint extension. The magnitude 
of the DWs close to the origin will in general depend upon the two-body scattering length. 
As in the attractive inverse square case, the relationship between self-adjoint extensions and 
Pt is not one-to-one. All observables are invariant under the transformation, 

p, ^ (5.58) 

where sq = 1.00623 ... is the magnitude of the imaginary solution of eqn. (15.521) . 

The Green's function, since it solves the Schrodinger equation, also observes the be- 
haviour (15.571) . We write the Green's function matrix elements evaluated at small hyperra- 
dius as 

(R, a\02B(p)\R, a) = MDg {yip. In pi p.) \ipm^ sin^ (5,, In p^R) , (5.59) 

where 7=1 /a2 is the binding momentum of the two-body bound state. Because of the 
relationship ( 15. 581 ). DG{x,y) is periodic in y with period Inlso- AH the analytic properties 
of Qib{p) as a function of /?, described earlier, are conferred on the function Do- Namely, 
there is the elastic two-body scattering cut at p = iy that runs down the imaginary p axis 
then along the real axis and there are poles at each of the bound states. 

The singular behaviour at small hyperradii that leads to the trig-log behaviour in the 
DWs also means that there is no ground state. At very high momentum the two-body scatter- 
ing length becomes insignificant and the DWs take the form they would in the limit y ^ 0. 
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This means that the deeply bound state take on the geometrically spaced form of the attrac- 
tive inverse square potential's states. At low energies the value of y is important and disrupts 
the exponential ladder of bound states; the Efimov effect no longer occurs. 

The lack of a ground state means that, as in chapter 4, we must truncate the full Green's 
function in both the continuum and bound states. Notice there is no need to continue this 
truncation into the elastic threshold since this is controlled by the low-energy scale y. The 
result of truncating the bound states is slightly more complicated than in the attractive in- 
verse square example because the position of each bound state is not only controlled by 
the high-energy scale but also by the low-energy scale y. The effect this has on the 
renormalisation group flow is important and requires careful consideration. 

For the individual distorted waves we write down the forms for small R as 

= Dj-,^,\n-^]\<p(atsm^solnp,R), (5.60) 
\P P P*l ^ ' 

|^p,,(^,a)|2 = p2)2(-,ln-^)|^(a)psin2(Solnp*^), (5.61) 

\P P* I 

|¥„(^,a)P = pn{y,p.fVB{y,Pn{y,pM^{af^^ri^[so\^P*R), (5.62) 

where we have explicitly shown the dependence of the binding momenta, p„, on y and p^. 
From eqn. (15.281) it follows that ©^(y, Pn) can simply be found by taking the residue of the 
Green's function at the appropriate bound state pole to give, 

J^Biy, Pn) = —nDciylp, In pip.), p ^ ipn). (5.63) 

Pn 

The projection operator matrix element for small R and p^ > follow from its definition 

a\np)\R, a) = pC U, In ) \ip(a)f sin^ (^o In pS) , (5.64) 
\P P* I 



where 



c{lM^\ = D^UMi-\^l r dkD,{K^-M^\. (5.65) 

" P*l \P P*l ^Jq \P P P*, 



The function C may also be written using eqn. (15.321) as. 



C|-,ln — 1 = Ini 



DG\-M-^\-DG\—M-^+in 

P P* \ P P* 



(5.66) 



P*, 

We are now in a position to rescale the differential equation for V-ig to obtain the DWRG 
equation. The scales p, y, k and k' rescale in the usual way. The three-body force is rescaled 
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by the relation, 

Vssip, r; k k';A) = ^\<P(mf sin' (^o In p.R) VU^p, Ay; A^, A^'; A). (5.67) 

Using eqn. (l5.24l ) with the addition of the truncated bound states the resulting DWRG for 
the three-body force is, 

d%B ^dV-iB ^dVsB fdVsB f^dVjs 
1 



1 -p^ 



VwiP, r; h if, A)V3b{p, r; iy, k'; A)2)2(r, in A/p,) 



+- [ dk"V3B(p,7\k:k''-,A)V3Bip,f k'\k';A)D3{k'\y,\nA/p,) 
^ Jo 

^S{A- p„{Ay , p*))A, 
1 + p^ 

(5.68) 



^ V f/ r -/a am"/ -/a 7' AN^(^-Ai(A7,p*))A , , . 

-T > , V3b(p, r; k, i/3; A)V3b(p, r; i/3, k ; A) Vsiy/Pn, Inpjp,). 



where we have suppressed the obvious functional dependence of Vsb in the terms on the 
first line to save space. 

The perturbations around the trivial fixed point solution, V^b = 0, are given by, 

(. \I+2m+2n+2n' 
^ • (5.69) 

The leading order perturbation is marginal, this is consistent with the trivial fixed point 
solution in the attractive inverse square (zero y) solution. Like the attractive inverse square 
DWRG this solution means very little by itself and in order to determine the nature of the 
marginal perturbation and the general flow of the DWRG we must examine a more general 
solution. Like the attractive inverse square DWRG the marginal perturbation occurs because 
the basic loop integral solution is an example of a limit-cycle solution. 

The Limit-Cycle DWRG solution. 

We look now for a solution to the 3BDWRG equation that depends on A logarithmically 
and that is independent of the momentum terms. This satisfies the 3BDWRG equation, 
.5 1 I d \ C(y,lnA/p,) n^6(A-p„iAy,p,))A 

(5.70) 
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Figure 5.2: The singularity structure of the integrand in the limit cycle solution. The bound 
state poles occur at § = ip„, the two-body elastic scattering threshold branch point occurs at 
q = iy. The propagator poles occur at ^ = +p. The contour C runs in a semi-circle from -1 
to 1 through the upper half of the complex plane and crosses the imaginary axis at § = /. 



The solution to this equation with logarithmic dependence upon A may be constructed in 
much the same way as the solution to the attractive inverse square potential DWRG. We 
write the solution as the integral, 



where the contour C is illustrated in Fig 15. 21 it runs from -1 to 1 through the upper half 
of the complex plane and crosses the imaginary axis at ^ = i. Substitution and use of 
eqn. (15.661) shows that it satisfies the continuous part of the 3BDWRG equation. The integral 
is clearly analytic about p ^ since the propagator poles to do not approach any part of the 
contour C in this limit. The integral is also analytic as y ^ since it matches continuously 
onto the attractive inverse square solution in that limit. The discontinuities caused by the 
truncation of bound states result from allowing the bound state poles to cross the contour. 
That the residues of these poles produce the correct discontinuities is readily confirmed 
using eqn. (15.631 ). 

The migration of the bound state poles up the imaginary axis is more complicated in 




(5.71) 



c c 
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this example than in the attractive inverse square case. The position of each bound state 
in the unsealed theory is determined by the scales y and p^,. In the rescaled theory they 
are determined by y and the ratio Ph./A. This has the interesting consequence that as A 
varies only that element of the rescaled poles position controlled by the high energy scale 
Pt changes. The element controlled by y, as low energy physics, remains fixed. 

We know that the transformation p^e"^'''' leaves the bound state spectrum un- 

changed. This transformation is equivalent, in the rescaled spectrum, to 

A ^ Ae-''"^ (5.72) 

which must therefore leave the rescaled spectrum unchanged. 

As A decreases the poles move up the imaginary q axis. To ensure the rescaled bound 
states poles are invariant under transformation (15.721) . new poles must 'condense' out of the 
elastic threshold cut when necessary. 

The solution (15.711) describes a logarithmically evolving limit cycle solution. This can 
be seen by observing that since DG(x,y) = DG(x,y + 2n/so), 

Vf^{p, y. A) = Vf^{p, y, Ae^"""'). (5.73) 

This cyclic behaviour means that there must be an infinite number of geometrically spaced 
discontinuities as A goes to either zero or infinity. At first it seems strange that we should 
be able to produce an infinite number of discontinuities as A ^ because for any particular 
A there are only a finite number of bound state poles inside the contour C in the complex 
^-plane. However, bearing in mind the discussion above we can see that as A ^ more 
rescaled bound states appear and, because of eqn. (15.721) . these will cross the contours in a 
strictly geometrically spaced pattern. 

The limit cycle behaviour is only exact in the rescaled theory. The unsealed potential 
will approximate the limit cycle when A » y but as A becomes less than y it will deviate 
significantly from it. 

The momentum independent perturbations around the limit-cycle solution take the 
form observed in chapter 4, 

1 1 °° / A \2n+m 

j-(p, y, A) = j^Xp, y,A)+y C^n^n t- ^"V"- (5.74) 
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The marginal perturbation Co,() must, as in the inverse square case, embody the same degree 
of freedom as the self- adjoint extension defining scale, p*- The DWRG flow is very much 
as described in section 3.3 and as illustrated in Figs. 14.214.31 

This momentum independent solution may be substituted into the three-body distorted 
wave Lippmann-Schwinger equation (15.81) to obtain results for the rss-matrix. The analysis 
follows much the same form of that in section 2.6. In this case the truncated Green's function 
also includes the integral around the branch cut along the imaginary axis, which corresponds 
to elastic bound pair and particle scattering. The result is 



£>3 -,-,lni- ©3 -,^,lni- 

P P Pj \P P P. 



1/2 



\p pJ 

~ 2-i ^2n,"' ^2n+m (5.75) 

m,n=0 



for a general 3^3 scattering. For elastic scattering of a bound pair and a single particle we 
may write the Ti,b matrix element in terms of a correction to the phaseshift in accordance 
with eqn. (ITT?ll . The DWRG limit-cycle potential then yields the DWERE, 

Dj^,ln^ ycot~6,s-i)+4nDoll,ln^\ = -l J] C.,.^. (5.76) 



If we have <0 then by eqns. (15.6515.661) we have, 

4nlmDGi-,\n^] = D2i-,\n^], (5.77) 
\P P*l \P Pj 

so that the imaginary parts on the LHS of eqn. (15.761) cancel and the correction to the phase- 
shift, 63B is real. If > then eqn. (15 .771) no longer holds because of the 3^3 states. 

The results above are for the momentum independent solutions of the DWRG only. 
More general solutions to the DWRG equations can be obtained by perturbing about the 
limit cycle solution with momentum dependent perturbations. These perturbations allow 
the three-body potential to difi'erentiate between difi'erent types of DW. It is clear that they 
will be important in amplitudes like (15.751) . However, in amplitudes such as (15.761) they will 
not provide any extra degrees of freedom in the power-counting since the external momenta 
k and k' are set to iy in this amplitude. For completeness the momentum perturbations are 
found in Appendix iDl 
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In any particular system that the EFT is describing the two-body bound state binding 
momentum, y, is not adjustable. It is therefore impossible to determine the three-body force 
terms that correspond to these perturbations around the limit-cycle solution. For practical 
purposes these terms can be absorbed into the energy dependent perturbations to obtain for 
elastic scattering, 

DjZ,in^ ](cot~S,B-i) + 4nDG{-,\n^] = -- V (|-) ■ (5.78) 

The power-counting in the energy dependent terms for the three-body force about the 
limit-cycle solution is precisely that found by Bedaque et al. The leading three-body force 
term is marginal and occurs at order (2/Ao)°. This leading order force fixes the phase of 
the three-body DWs close to the origin, or equivalently forms a self-adjoint extension. The 
use of the DWRG method in deriving the power-counting for the three-body force in the 
EFT is a new result. It provides support for the work of Bedaque et al [22], which relies 
on introducing a dimer^ field into the EFT Lagrangian and then practical use of the three- 
body force to renormalise the theory order by order. Despite the power of their approach 
Bedaque et al cannot provide a simple algebraic statement of the power-counting nor a 
simple formula like (15.781) . 

When looked at from a different angle, our analysis means that the introduction of any 
three-body data to fix a self-adjoint extension such as was done by Phillips and Afnan [33] 
is equivalent to a LO three-body force. It is not necessary to explicitly calculate the three- 
body force. However, if we were to do so after fixing a self-adjoint extension we would find 
that the degree of freedom associated with that three-body force would be equivalent to that 
associated with the self-adjoint extension choice. 

The RG methods used here are very different from those used by Hammer and Mehen 
[23], which start with the STM equation and supplement the approach of Bedaque et al 
[22,24]. The STM equation can only describe three-body bound states and 2-1-1 scattering 
and so does not include all of the three-body physics. Furthermore, since the equations of 
Hammer and Mehen start from a truncated STM equation, part of the role of their three- 
body force is to restore the truncated part of the STM equation, a complication the DWRG 
field with the quantum numbers of two particles, essentially used as a bound state. 
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method was designed to avoid. For similar reasons the form of our three-body force can not 
be directly compared to that given by Bedaque et al. 

A limit cycle form of the three-body force in this EFT has been suggested by many 
previous papers [30, 32] that have essentially tackled the infinite scattering length case of 
the attractive inverse square problem. Bedaque et al also found a limit cycle solution but 
because of their approach could only find the force order by order [24]. 

Our result is also applicable to the three nucleon force in the pionless KSW EFT [22,24, 
65]. Braaten and Hammer have discussed the tuning of the quark masses in QCD required 
for an infra-red limit cycle in the three nucleon problem [65]. They assume that to obtain 
the scale invariant relationships that define a limit-cycle, such as (15. 721) . the system must be 
tuned to the infinite scattering length limit to avoid the introduction of scales that break the 
symmetry. The analysis here, which uses a Wilsonian RG approach, readily produces limit 
cycles in the rescaled potential since the low-energy quantities do not affect the symmetry. 
The discussion of Braaten and Hammer [65] applies to our unsealed potential, which only 
exhibits infra-red limit cycle behaviour in the infinite scattering length limit. 

The complete 3BDWRG analysis in this system provides a concrete grounding for what 
has been intuitively understood. Beyond that, because of the common method of solving the 
DWRG equations in in this chapter and the previous two, the connection between all these 
results can be seen. 

5.6 Summary 

In this chapter we have examined the DWRG for three body forces. The approach used can 
be generalised quite simply to A^-body forces. The 3BDWRG equation is complicated by 
this multi-channelled problem because the 3BDWRG has to describe the coupling of the 
three-body force to each of the different distorted waves. 

The application of the 3BDWRG to the case of "well-behaved" potentials was included 
as an example. The two fixed points found correspond to Weinberg counting and to the three 
body equivalent of the KSW scheme. The instability of the non-trivial fixed point leads us 
to conclude that in almost all systems of this kind Weinberg counting and naive dimensional 
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analysis is appropriate. 

The 3BDWRG was also used to study three body forces in the KSW EFT. The solutions 
to the equation take the form of limit cycles and have much in common with the solutions 
found in chapter 4. The power counting found is marginal at leading order and matches that 
found by Bedaque et al [24]. 

Although our analysis has only considered the system of three Bosons, it is also appro- 
priate for three nucleons in the i channel as we shall see in the next chapter. 



Chapter 6 

The KSW Effective Field Theory of 
Three Bodies 



A pionless KSW EFT for three nucleons has been successfully described by Bedaque et 
al [24,59,66]. Their work also covers the related system of three Bosons. We have seen 
in the previous chapter how the 3BDWRG analysis of the three-body force in this prob- 
lem supports their conclusions. We shall now look to complement this work by deriving 
equations for the three body DWs. 

All the predictive literature on the three body KSW EFT is concerned with scattering 
observables [24, 33]. However, some applications of the EFT require knowledge of the 
wavef unctions. For example, if we want to model scattering of an electron from the triton 
we need to know its charge density, which can be ascertained from the wavefunction. 

The explicit solution of the equations for the DWs will also serve to illustrate the 
'phase-fixing' role of the LO three-body force. We have seen in the 3BDWRG the equiv- 
alence of the choice of self-adjoint extension and the LO marginal term in the short range 
force. Because of this equivalence, the introduction of the three-body force at LO is a trivial 
matter of choosing a boundary condition for the poorly defined DWs. 

As an introduction we shall consider the example of three Bosons. We shall then move 
on to the physically more interesting example of three nucleons. 



1 r\r^ 
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6.1 The EFT for Three Bosons 

The distorted waves for elastic scattering of a bound pair and a particle are not particu- 
larly easy to find using Efimov's boundary conditions, although Fedorov and Jensen have 
attempted to develop a method of doing so [60]. We will derive equations for the DWs using 
a useful property of systems with zero range interactions. Namely, that the value of the DW 
may be simply determined from the value of the DW when the relative coordinate of a pair 
of particles is zero. Because the potential is zero at any point other than where two particles 
are 'touching', the value of the DW anywhere else may be found by integrating the wave 
equation. Skorniakov and Ter-Martirosian (STM) [56] were the first to take advantage of 
this property in their derivation of equations for the 2+1 scattering amplitude for zero range 
forces. 

Our equations for the DWs will also use this property. We will derive equations that 
may be used to find all DWs below the three-body threshold. At one point in the derivation 
we will be able to link our method to the STM equations. Unlike Skorniakov and Ter- 
Martirosian, we will start our derivation with the Faddeev equation (15.1 IL which we rewrite 
here for convenience, 

\W^) = {1+P)\ilfl), (6.1) 
\K> = i^ko) + Gl(p)np)P\,lft,)^ (6.2) 

where we have now labelled the wavefunctions by ko, the vector of the third particle in the 
centre of mass frame which has magnitude, |ko| = ■yj4{p^ + y^)/3, and = ME < is the 
centre of mass energy, relative to the three-body threshold. Initially we will derive equations 
for the LO EFT, hence the two-body T-matrix is given simply by the scattering length term 
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in the effective range expansion^ 

<k, l\t(p)\k', V) = (2nf6\k - k'Mp, k), (6.5) 



T(p,k) = 



M 



-y + i\lP^--M ■ (6-6) 



To make the importance of the zero separation value of the DWs apparent we will work in a 
mixed co-ordinate system (k, x), i.e. in terms of the distance between two particles and the 
momentum of the third particle in the centre of mass frame. Two matrix elements that will 
be needed are: 

<k, x|Gj(/7)?^(/7)|k', x') = T{p, k')Li(p, k, k',x)6\x'), (6.7) 

M 



Li(/?,k,k',x) = -(2;r)V(k-k')— e'"^^^, (6.8) 

4nx 

<k, x\PG^(p)r(p)\k'x') = T(p, k')U{p, k, k', x)6\x'), (6.9) 

cos[x.(k' + ik)] 

The correctly normalised two-body bound state wavefunction is simply given by 

X{k,x) = (2;r)353(k - k,)^e~y' = -^L, (p,k,ko,x) . (6.11) 

4nx M 



6.1.1 Momentum Space Equations 

As a starting point, we may use the Faddeev equations to write an integral equation for 
i/r(k, x). Inserting the values for the matrix elements given above we get, 

JTv r d%' 

^(k,x) = — y-Li(p,k,ko,x) + J —T{p,k')U{p,k,k',x)^(k'), (6.12) 



and y are defined, as before, by 

Xij^Ti-rj, y;t = r^-^(r,+r^), (6.3) 

where {/, j, k} is an even permutation of {1, 2, 3} and ri, r2, rs are the absolute coordinates of the three 
particles. The momentum conjugate to these we label as 1,^ and and are given by 

hj = q< - q;, kk = (ik- ^(q« + qy), (6.4) 

where {/, j, k} is an even permutation of {1, 2, 3} and qi, q2, q3 are the absolute momenta of the three 
particles 
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where 

^(k) = ¥(k, X = 0) - iA(k, X = 0). (6.13) 

Already, the importance of the DWs at x = is apparent. ¥(k, x) is found by operating on 
(16.21) with 1 + P. This gives, 



mx) = -4^L(;?,k,ko,x)+ r ^T{p,k')L{p,kX,x)^{k'), (6.14) 



where L = L\ + L2. These equations may be used to find a one-dimensional equation for ^. 
Subtracting (16.121 ) from (16.141) and setting x = we get, 

m = -^L2(p,kM,0)+ [ ^^Tip,k')L2ip,k,k\0)^ik'). (6.15) 

Equations (16.151) and (16.141) may now be solved to find Our aim is to derive equations for 
the DWs in the coordinate space representation but these equations are not in a convenient 
form for converting into coordinate space. To find more suitable equations, let us define the 
projection of the wavefunction as 

0(k) = limxWiKx) = 2!^{2nf6\k - ko) - ^T(p,fc)^(k), (6.16) 

where the second identity follows from eqns. (16.816.1016.141) . Substituting this into equa- 
tions (16.141) and (16.151 ) we get. 

An C d^k' 

mx) = -— — — L(p,k,k',x)0(k'), (6.17) 
M J (2ny 

T(ky'm) = J ^L2(p,k,k',0mk'). (6.18) 

These equations were first derived by Skorniakov and Ter-Martirosian (STM) [56] in a 
somewhat different manner. The projection of the wavefunction, O, is the value of *F when 
two of the particles are at zero separation. It can, with caution, be thought of as a wavefunc- 
tion that represents the penetration of the third particle into the bound pair. Mathematically, 
O is simply a tool that allows the full DW, to be found using a one-dimensional equation. 
An equation for the elastic scattering amplitude can be obtained by assuming an asymptotic 
form for O, 

y //o-nScs... 1. n , ^tt a(ko,k) 




The resulting equation has been studied by Danilov [57] and has more recently been derived 
by Bedaque el al [22, 24] from the KSW EFT Lagrangian. 
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6.1.2 Coordinate Space Equations 

The STM equations may be Fourier transformed to obtain equations that may be solved to 
find the distorted waves in the coordinate representation (see Appendix El). We define, 

,((k.y-k'.y') 



V3^2 ^ ^^(^Q^) k^^kQ^) KM-) \ ... 

+ — tt; — + — -7^^ — r ' (6.20) 



Vo(y,y') = - 



r d% r d^k 

3mJ (2^ J W 



16;r r d^k C d^k' [ /3~ I 



L2(p,k,k',0)e'*-^-'''^') 



p3^ — + 4|r--j— |, (6.21) 



2V37r: 



where, 



<3o = J + (y - y')^ <3± = y ^-^^ + / + + y-y' ± ^^.y', 



Kn{z) is the modified Bessel function of the third kind. In the case of positive energy, the ie 
prescription gives k = -Hp I Vs so that each of the modified Bessel functions go to 

KM^^'if^'n'^n^W. (6.22) 

Using the transforms we may find equations for *F(x, y) and 0(y). A simple Fourier trans- 
form of eqn. ( 16.171) yields. 



¥(x,y) = JdVW^(y,y',x)0(y'). (6.23) 

In order to Fourier transform eqn. (16.181) we first multiple both sides by 

16;r. 



3M 

and then transform to obtain. 



(r + V3^V4 - p2), 



(V^ + k"^ 0(y) + J dV V()(y, y')O(y') = 0. (6.24) 

This looks very much like the Schrodinger equation with a non-local potential. This form is 
very suggestive of the interpretation of O mentioned above, namely a wavefunction describ- 
ing the penetration of the third particle into the bound pair. However since the potential is 
not symmetric, and so is not hermitian, this interpretation should be used with caution. 
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The s-wave wavefunctions are simply found by integrating over all angles, 

'¥(x,y) = xy J ^^^C^'y)' W=jJ ^«^(y)- (6-25) 

By integrating eqns. (16.2316.241) over all angles we can obtain equations for the s-wave 
wavefunctions: 

J/-»oo 

dy'w(y,y,x)(^(y'), (6.26) 


^ + fc^ j 0(3;) + dy'voiy, y')<!>(y') = 0, (6.27) 



where, 







27T ' 


\ Qo 



w(y,y',x) = 



y[?>ny' 

^ ,^ l^K^y (K2(kR.) K2(kR^) 



Qo 

+^^iKo(KQl) + Ko(kQZ) - Ko(kQ1) - Ko(kQ-)\ (6.28) 



and 



R^ = ^/y^+y^±yy, (6.30) 



2o = ^jlx^+y^+y^±2yy, Qt= ^J^x^+y^+y^±yy±^xy. (6.31) 
These s-wave results follow from the indefinite integral result [41], 



Kn(z) Kn-liz) 

z" ^ z" ' 



These equations can be used to find all s-wave DWs to LO in the EFT with total centre 
of mass energy less than zero, i.e. all bound states and all bound pair and particle interactions 
below the three-body threshold. If the centre of mass energy becomes greater than zero, k 
becomes imaginary and the modified Bessel functions in the 'potential' become Hankel 
functions. Since Vo(j, j') will then have oscillatory behaviour for large y, solution of the 
equations becomes extremely difficult. 
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6.1.3 Analytic Solution 

In general eqns. (I6.26I6.27I) have to be solved numerically. Before we look at their general 
solution we shall look at their analytic solution in certain limits. Eqns. (16.2616.271) may be 
solved analytically in the limit a:, 7 ^ as they take on the far simpler forms: 



The solutions to these equations correspond to the short range behaviour of the DWs. To 
solve them we assume an ansatz, (D(y) = y'\ which we justify by observing that since no 
other scales exist in the equation, the wavefunction must satisfy a power law. This argument 
may be put into a more rigorous form by considering Mellin transforms. Eqn. (16.331) gives 
the possible values of s which we find, by evaluating the integral, must satisfy the equation. 



This equation, as it should be, is precisely that obtained from Efimov's approach. It is in 
fact possible to show that, in accordance with the analytic results of Efimov's approach, in 
the limit of 7 ^ 0, but non-zero k eqn. (16.271 ) has the solutions K^iKy) where 5 is a root of 
eqn. (lOl) . 

As noted before, eqn. (16.341) has two imaginary solutions, s = +iso, with s = 4 the 
smallest real solution. Therefore, the solutions of eqn. (16.271) will be dominated by the so 
solution at small y. 

If we take 0(3;) = y\ the integral to find ^'(x,y) is far more involved. After some 
work (see Appendix 10) we find we may write the result in terms of the hyperradius and the 
hyperangle: 




(6.33) 




(6.34) 




when > |. 



when < |, 



(6.35) 



(6.36) 
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This result confirms that there is no good boundary condition for the DWs as i? ^ 0. The 
angular behaviour is continuous at a = 7r/6 but has a glitch. Since a = n/6 is where the third 
particle meets either of the particles in the pair, this glitch corresponds to their interactions. 

The behaviour of the solutions, O, of eqn. (16.271) for large y depend simply on the 
relative values of k and y. If > and k > then the system is above the threshold 
for elastic scattering of a bound pair and particle and below the three-body threshold. The 
projection, (!)(y) sm(kQy + 6), where 6 is the phaseshift. The 'potential', VQ(y,y'), is real in 
this limit so the phaseshift is also real. Using eqn. (16.261) we may find ¥(jc, y) in this limit. 
We find, 

"Vix, y) sinikoy + 5)e-^\ (631) 

which describes a non-interacting free particle far from a bound pair. 

If < then the system is below the two-body threshold and the only solutions are 
the bound state solutions (^(y) e'^'^^-. 



6.1.4 Numerical solution 

Solving eqn. (16.271) with some fixed phase in the trig-log behaviour close to the origin is 
equivalent to the three boson EFT to order (2/Ao)". The LO three-boson force is taken to 
be the mechanism by which the phase is fixed in accordance with the analysis in chapters 4 
and 5. 

In order to solve eqn. (16.271 ) numerically we must find some way to deal with the 
singular behaviour close to the origin and to fix the phase of the solution there. We introduce 
some scale Q Aq and define the solution to be. 



0(3;) = 



sinisoMyy) + T]), y < Q. 

(6.38) 



The three-body force may now be chosen by simply choosing the phase rj of the solution. 
Inserting this into eqn. (16.271) gives a homogeneous equation for O^j, 



dy^ 



dy'(!>n(y')vo(y, y') = - dy' sm(so \n{yy') + T])vo(y, y'), (6.39) 
a Jo 
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with the boundary conditions 



On(Q) = sin (^o In yO^ + T]), 



So 

= — cos ( In yQ. + rj) . (6.40) 

y=n 



dy 

Q. is simply a numerical tool and, up to numerical noise, results should not depend upon it. 

This equation may be solved using linear algebra methods, in practice the upper limit 
in the integral must be made finite, however, because of the exponential suppression of the 
kernel vo(y,y') for large y the effect upon the result is minimal. The first step in solving 
numerically is discretisation of the variable y. Since the solutions of the eqn. (16.271) depend 
logarithmically on y at small y and linearly at large y, the discretisation of the variable must 
reflect this. We define the points y„ by the relation, 

u„ = soln(yy„ + e'°y"'">-l), (6.41) 

where m„ varies linearly onn. At small y„ we have, 

u„ ~ So In yyn, (6.42) 

so that the y„'s are logarithmically spaced allowing them to describe the trig-log behaviour 
of the solutions in that region. At large y„ we have, 

Un ~ koyn, (6.43) 

so that the y,, are linearly spaced allowing them to describe the ordinary trigonometric or 
exponential behaviour there. In order to set the two boundary conditions at the points yo and 
yi we set y2 = fias the first point at which the value of the solution is unknown. 

Using this discretisation of the variable y, eqn. ( 16.391) can be written as an inhomoge- 
neous matrix equation. We define O,, = <l>(j„+i) so that the derivative term is found using 
the difference, 

"^'"^^"^ ^ ^^0(j„,i) - \ 0(3;,,) + -J— 0(3;„„i)) (6.44) 



dy^ 2(A„ + A„_i) \A„ (A„ + A„_i) 2A, 



71-1 



= £)„,„-2<I>„_2 + £>,v,-lO„_i + Dn,n^n, (6.45) 

where A„ = y,j+i - y^. The last line defines the derivative matrix D, y, 1 < i, j < N in 
which all other elements are set to zero. Notice that our definition also defines three other 
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elements, 2)i,o,2)i -i and £)2,o that do not appear in the matrix; these terms are required to 
set the boundary conditions at yo and yi . 

The integral term can simply be written as 

N 

dy'vo(yn,y')^(y') = J] (6.46) 

i=l 

where the coefficients, "Vij form the matrix "V and can most easily be found using Simpson's 
rule, in which case A'^ must be chosen to be odd. 

Altogether the eqn. (16.271) is written in matrix form as, 

N 

-^u((^i,o + ^o) sin(5o In 73^1 + n) + sin(5o Inyjo + ??)) 

-^;,2^2,o sin(5o In yyi + rj), (6.47) 

where, 

In = -\ dy'vo(y„+i,y')sm(so\nyy' + 7]). (6.48) 
Jo 

The second and third lines of eqn. (16.471) define the boundary conditions by ensuring that 
the solution matches on to the form chosen for y < yi- 

For scattering solutions, > 0, we solve the equation numerically for two values of rj, 
T] = and T] = n/2 for instance. All solutions with different values of rj can be constructed 
by superposition. 

Bound state solutions can be found for any energy by simply solving the equation for 
the two different t/'s then finding the superposition that vanishes for large y. This gives a 
constraint on the interior phase of the bound state solution at that energy. Only the bound 
states with the correct interior phase - chosen by the LO three-body force are the physical 
ones. 

Two sample solutions of eqn. (16.271) are shown in Figs. 16. 1116.21 Fig. 16. 11 shows a scat- 
tering solution to the equation with incoming momenta ko = y/2. The figure shows the 
solution on both log and linear plots so that the interior and exterior behaviour of the wave- 
function is visible. In the interior the trig-log behaviour is evident, externally the solution 
goes to sm{koy + 6) where 6 is the phaseshift for the elastic scattering amplitude. The full 
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*{p) 




*(p) 




Figure 6.1: Numerical 2+1 scattering solution to eqn. (I6.27I) with relative momenta ko 
7/2 and total centre of mass energy E = -13y^/(16M) 



*(p) 





Figure 6.2: Numerical solution to eqn. (16.271) . Example of a shallow bound state solution 
with centre of mass energy E = -lOly^/ClOOM). 
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Figure 6.3: Three-body distorted wave, numerical solution to eqns. ( I6.26l6.27t with relative 




Figure 6.4: The small R detail of fig. 16.31 matching the anticipated result, eqn. ( 16.351 ). 
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DW, *F(x, y) which results from this solution and was obtained using eqn. (16.261) is shown in 
Figs. 16.316.41 

The full DW solution is plotted on a hyperpolar scale (i?, or). The hyperradius R can 
be considered the mean distance between the particles and the hyperangle a describes the 
configuration of the three particles. Evident in the plots is the symmetry of the wavefunction 
with respect to the interchange of particles. In the fig. 16.31 the hyperradial scale is linear. 
There are two outgoing states at a = n/2, corresponding to the 'chosen'^ bound pair, and 
at q; = n/6 corresponding to the other two possible pairs. The discontinuity across the 
derivative at a = 7t/6, anticipated in eqn. (16.351) is also evident. Fig. 16.41 shows the small R 
detail of the DWs, which exactly matches the analytic result, (16.351) . 

Fig. 16.21 shows a bound state solution to eqn. (16.271) . again showing detail on both log 
and linear plots. The interior phase of this solution was determined by the choice of binding 
energy. From an EFT perspective the three-body force determines the interior phase which 
in turn determines the binding energy. 

Fig. l6.5l shows physical observables determined from solutions to the LO EFT equation 
(16.271) . Fig. l6.5r A) shows the possible different phaseshift curves for the elastic amplitude 
that result from different choices for the LO three-body force. Each of these curves cor- 
responds to a different value for the scattering length (Fig. l6.5r B)). The range of possible 
scattering lengths that can be obtained from different three-body forces is -oo < ^3 < oo. 
Notice that it is only because the three-body force occurs at such a low order, (2/Ao)° that 
it is able to affect the scattering length of these amplitudes so completely. Fig. l6.5r C) shows 
the shallowest bound states that result from our choice of three-body force. Figs. l6.5t B) 
and l6.5r C') together imply that there is a single parameter relationship between the elastic 
scattering length and the binding energy of the shallowest bound state, which is illustrated 
in Fig. l6.5r D). This curve cannot be seen experimentally as the physical world corresponds 
to just one point on it. However, the equivalent curve in the three nucleon problem, known 
as the Phillips line [22,24,61] has very interesting implications. 
^Chosen in the sense of our definition of a. 
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Figure 6.5: Numerical Results of the LO EFT for three-bosons. (A) All possible phaseshift 
curves that result from different values of the leading three body force. Each value for the 
LO three body force corresponds to a different interior phase rj and results in a different 
scattering length for 2+1 elastic scattering (B) and the shallowest bound state (C). Together 
these imply a single parameter relationship between the scattering length and binding energy 
of the bound states (D). 
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6.1.5 NLO Equations 



At NLO, (Q/AqY, in accordance with eqn. (I5.76I) there is no extra three-body force term' but 
there is an effective range correction in the two-body force. The equations for the effective 
range corrections can be found using eqn. (16.141) and the form for the two-body T-matrix 
with effective range corrections, 

An 1 



rip,k) = 



M 
An 







2 


7 + 




-1- . . . 









(6.49) 



This modified form of the effective range expansion ensures that there is a bound state pole 
at iy. Substituted into eqn. (16.141) and keeping just terms of order (Q/Aq) (i.e. neglecting 
terms like p^) gives, 

\6n( 



(ki-e)m) 



Stt 



3M 



7+ \l^k^-P^ 



V 



7+ \l^k^-P^) 



r d'k' 
J (2^ 



:L2(p,k,k',0)O(k') = 0. (6.50) 



Fourier transforming this equation as before we obtain, 

{y^ + kl)(l>{y) + J dV(yo(y,y')+P.Vi(y,y'))0(y') = 0, (6.51) 
where (see Appendix |E|), 

d'k' 



Stt r d'k r d%' 



7+ \l^k^-P^ 



L2(;?,k,k',0)e'*''-'''^') 



In^ 



8 '^-\yy--^ 



2y( 2\ K2m\ 



(6.52) 



'Recall that the general expansion of the three-body force contained a term at this order but 
proportional to y. This will in fact be seen as a small correction to our LO 'phase-fixing' three body 
force. 
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The s-wave equivalent follows quite simply, 



^ + klj(!>(y) + j^^ dy'ivoiy,y')+PeVi(y,y')<!>(y') = 0, (6.53) 



where 



Ri Rl 



We can solve eqn. (I6.53I) up to order pe using first order perturbation theory. If we define 
^(y) = 0(°^(j) + Pe^'-^Ky) and <l)^")(j) satisfies the LO equation (16.271) then up to order p^, 
O^'^Cy) satisfies 

(^+-tgJo(i)(3;) + d/vo(j,/)0(V) = -J^ vi(3;,/)0(V)- (6-55) 

This equation may be solved numerically by adapting the method used for the LO equation. 
Since no new parameters are entered at this order the results of solving the equations to this 
order give a refinement to the one parameter relationship between the 2+1 scattering length 
and the three particle binding energies. 



6.2 The Pionless EFT for Three Nucleons. 

One area where the two-body EFT for short range forces has proven very successful is 
nuclear physics. As noted in the earlier chapters, nuclear physics provides an excellent and 
physically important example of the KSW power-counting scheme in action. In this chapter 
we will extend the results of the previous section to look at the pionless KSW EFT for three 
nucleons. 

The spin and isospin degrees of freedom available to the nucleons as well as Fermi 
statistics makes the problem slightly more complicated then the analysis for Bosons. To 
avoid the complication of electromagnetic interactions we shall concentrate on the system 
of two neutrons and a proton"*. 

"^This system is clearly more interesting than the other electromagnetically neutral possibility of 
three neutrons as it contains a well-known two nucleon bound state, the deuteron, and three nucleon 
bound state, the triton. 
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If the total spin of the three nucleons is 3/2 then the three spins are aligned and the 
Pauli exclusion principle acts as a repulsive force between the two neutrons. This system 
is insensitive to short range physics and has been successfully described by an EFT without 
three-body forces [66]. However, if the total spin is 1/2 then, because it is possible for the 
three nucleons, now non-identical, to approach each other, the system is sensitive to short 
range physics, it is upon this system that we shall concentrate in this chapter. 



6.2.1 Wavef unction Symmetry 



We assume that the pairwise nuclear force conserves both isospin and spin so that these 
may be taken as constants, and we put the total spin and isospin and isospin component to 
S = T = T^ = j.lfwe identify one particular pair of nucleons, 1 and 2, then we can identify 
two linearly independent states of three particles with spin 1 /2 that may be used as a basis. 
Firstly, the singlet state in which the pair of nucleons are in the spin singlet state, s = 0, 
which we denote by and secondly, the triplet state in which the pair are in a spin triplet 
state, 5=1, which we denote by a>l^. The singlet state is antisymmetric and the triplet state 
is symmetric under interchange of 1 and 2. A spin basis may also be constructed using an 
alternative pair, 2 and 3 say, as the starting point for the basis. The transformations between 
the different possible bases are given by. 



CO 



23 



OJ 



23 



1 

"2 

V3 
2 



^ "1 

2 

_1 
2 J 



OJ 



12 



CO 



12 



31 



CO 



, ,31 



_1 V3 ^ 
2 2 

V3 _1 
V 2 1 ) 



CO 



(6.56) 



Similarly for isospin we have the singlet and triplet states, ■Q^^ and which satisfy the 
same transformations and symmetries. 

Since the pairwise potential is symmetric under interchange of the two interacting par- 
ticles the continuous part of the wavefunction must be symmetric under interchange of the 
particles in the pair. Thus a general state which is antisymmetric in 1 and 2 can be written 
in the form, 

= c^f ) + (6.57) 
The full wavefunction, which must be antisymmetric in the interchange of any pair can now 
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be written as, 

m) = \<p'^)) + \<p^')) + \ip'')) = (1 + P)\<p)), (6.58) 
where P is the permutation operator. To see that is totally antisymmetric notice that 

\<p'')) ^ 1^^')) = -1^'')), (6.59) 

under the interchange of 1 and 2. Since we have now introduced the permutation operator, 
P, we shall drop particle labels. Using the decomposition of we can write down the 
Faddeev equation for |^)), 

(Ho + V-E)\cp)) = -VP\cp)), (6.60) 

where the hat on the potential signifies that it is an operator in spin-isospin space. 

We shall assume that the potential V is diagonal in spin and isospin. (This assumption 
is approximate but is true to the orders at which we shall be working [22].) In the centre of 
mass frame of the two interacting particles the Pauli exclusion principle restricts V to just 
two non-zero matrix elements: 

(CO A, Vco,i9,) = y„ Vcor&,) = Vt. (6.61) 

The matrix elements, in spin-isospin space, of VP now follow from the transformation rela- 
tionships (I6.56I) . 

{ioA,VPcoA) = ^VsP, 
(oj,&,,VPoJt&s) = lv,P, 

(oj,&„VPajA) = ^V,P (6.62) 

Taking the inner product of eqn. (16.601) with iOs&t and cot&s yields coupled Faddeev 
equations for the and 

(Ho + V, - £)I«A.v) = -V.i'l^liA.v) + , 

(Ho + V,- E)\i/f,) = -VtPl^lilf,,) + , (6.63) 
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which can be written in the form, 



\ilf,) = Gl{E)t,{E)P\\\ili,) + ^-\ilf,)^, 
m = \Xt) + Gl{E)t,{E)P\^\ilf,) + 



(6.64) 



where \xt) is the asymptotic state with a bound pair (the deuteron) and free particle. There 
is no bound pair in the singlet channel so there is no asymptotic state. 

Using eqns. (16.5616.581) we can now write the full wavef unction in its final form, 

P\. . 3P. 



I P\ iP 

|Y,) = (l + -)|,A.) + — l-A.), 

/ P\ 3P 
|¥,) = (l + -)|^,) + — 1^,), 



(6.65) 



where, (xi2,y3l^ = (x23,yi| - (X3i,y2|. From these equations we may now parallel the 
analysis for three Bosons to obtain the equations for the full wavefunction which are, 

yl3n 



T(kx)--i^ 

^«(k,x) = - 



L,{p, k, k', x)0,(k') + L,{p, k, k', X) ^ io,(k') + ^0,(k') 



M 



f ^^3(;?, k, k', x)[o,(k') + 0,(k')], 



Taip,k)-'^a(k)= f -^L2(;?,k,k',0) 
where a,b £ {s, t} with a b, and L3 is given by, 
L3(;?,k,k',x) = -2iM 



^cD„(k') + ^(D,(k') 



(6.66) 



sin[x.(k' + ik)] 



F + k'^ + k.k' - - ie 



(6.67) 



6.2.2 The Coordinate Space Equations. 

These equations can as before be readily converted into coordinate space. The equations for 
T(x, y) do not depend upon the two-body T-matrices and are given by. 



Y„(x,y) = J 



Wi(y,y',x)(D„(y') + W2(y, y', x) | io,(y') + ^0,(y')] 



„(x,y) = J JVW^3(y,y',x)|^0,(y') + 70,(y')|, 



4 



,(6.68) 
(6.69) 
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where 



Wi(y,y',x) 



V3^^ K2(kQo) 



W2(y,y',x) = 
W3(y,y',x) = 



4n^ { Ql Ql 

3/^2 (K2{kQ^) K2(kQ.) 



(6.70) 

(6.71) 
(6.72) 



167r2 [ Ql Q^ 

In order to get the equations for (b(y) we must define the two-body r-matrices. The T-matrix 
in the triplet channel is given by, 

4n 1 



T,(p,k) = - 



M 
An 



+ 



(6.73) 



This has a pole at p = ijt, which corresponds to the deuteron, hence jt is given by the 
deuteron binding momenta y, = 45.70MeV. The effective range is p, = 1.764fm. In the 
singlet channel there is no bound state so the T-matrix is given by, 

471 1 



Ts(p,k) = - 



An 



( r; A 



1 



My2_ 3^2+^2 



+-r. 



ys+ ^J^k^-p' 



y] 



^\k^-p' 

■ \ + 



'ys-^IW-p\ 



(6.74) 



The difference in the expansion makes the NLO term slightly more complicated. Here we 
have = = -23.714fm and r, = 2.73fm. The Fourier transforms now follow in much 
the same way as in the Boson case. We get. 



v' + 37a-^')^«(y) + 



JdV'(vo"(y'y') 



+ra(n(y, y') + 6asVHy, y') 1 1 ^0,(y') + ^(D,(y' | = 0, (6.75) 
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where V'^iy, y') is given by eqn. (16.211) with the substitution y ja, V"(y, y') is given by 
eqn. (16.521) with the substitution y ^ ja and 



3M 



r d^k r d^k' 



L2(p,k,k',0y*-^-'''^'\ (6.76) 



We were unable to evaluate this final integral analytically. It can be reduced to a one dimen- 
sional equation as in appendix El and then evaluated numerically. The s-wave equivalents of 
these follow easily. The term in eqn. (16.761) is neglected by Bedaque et al [24], who simply 
use the modified form of the ERE in the singlet channel and take = r^. The difference is 
comparable to NNLO corrections. 



6.2.3 Analytic Results. 

It is useful to rewrite the equations for and O, in terms of 0+ = + O^. The LO s-wave 
equations for these are then given by, 

($ ^ ^ ^ ^' ) " ^) ^^^^^ ^ I " ^' ) ^^^^ 

+ dy' |vo%,/)0,(3;') - ^v^(j,/)0_(/)j = 0, (6.77) 

where, 

Vo%,/) = \(vl(y,y') ± v',(y,y')). (6.78) 
The equations for 0+ are useful because at small y they decouple to give 

^^^-'^^^rd/a..(/) r^?y:^: =o. ,6.79) 

3;r Jo 



yl3n Jo "^-^ >4 + 3,2^2 + ^4)2 
where A+ = I and = -1/2. Substituting 0+ = as before we obtain the equation for s, 

cos f?:^) sinful = 0. (6.80) 



2 1 V3. V 6 

This equation for s when /I = 1 is the same as the equation that describes the Bosonic 
wavefunction at short distances and hence has the imaginary solutions s = +sq [22,24,25]. 
The equation for /I = -1/2 only has real solutions, the smallest of which is = 2. In solving 
eqn. (16.771) generally we may identify solutions by their behaviour at small y. There are 
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solutions in which the interior behaviour is determined by a boundary condition on 0+(>'), 
for example, 

0+(3;)^sin(^o In 7^^ + 77) + ..., (6.81) 
C>.{y) Ji(y, - jt) ^m{\ny,y + 77 + tr) + . . . , (6.82) 

where the constants J[ and a can be determined by solving eqn. (16.771) order-by-order in 
y. There are also solutions in which the interior behaviour is determined by a boundary 
condition on 0_();), e.g. 

0.(3;)^/ + ..., (6.83) 
0^(3;) ^ S(r, - yr)y' + .... (6.84) 

These solutions must be uniquely combined in a manner determined by the boundary con- 
ditions at 3; ^ 00. All the physical distorted waves will be dominated at short distances by 
the trig-log behaviour of the 0+ solution. 

The three-body force will prove to be very important in this system. The 3BDWRG 
analysis for the three Boson EFT may be extended to this example and its conclusions are 
very much the same. Namely the LO three-body force will be marginal and will correspond 
to a self-adjoint extension of the Hamiltonian. It will in effect fix the phase of the distorted 
waves close to the origin. The next three-body force term that can be determined^ occurs at 
the order (Q/Aq)^. The eqns. (16.771 ) plus a phase fixing condition correspond to an EFT to 
order (Q/A^)". 

At large y, below the three-body threshold but above the deuteron binding energy, the 
singlet solution should vanish. This constitutes an additional boundary condition on the 
equation that determines the combination of solutions described above. After applying this 
boundary condition, the triplet wavefunction has the long-distance form, 

^t(y) ^ ^Hkoy + S), (6.85) 

where = 4/3yj + is the momentum of the free neutron in the centre of mass frame and 
6 is the phaseshift for neutron-deuteron scattering 

^Recall that the terms that scale with y cannot be measured since we cannot change the value of 
y, they are therefore absorbed into the energy-dependent force couplings. 
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In solving the NLO equations no new three-body force term is required. The 'yt,ys 
dependent terms at this order in the three-body force will appear as refinements in the LO 
'phase-fixing' three-body force. At this order there are two equations to solve, correspond- 
ing to the perturbations, and Pf Writing, 

^±iy) = of (y) + (r, + A)OL^^^ + (r, - p,Wt' + • • • , (6.86) 
where 0+ ^ satisfy eqn. (16.771) . using first order perturbation theory we obtain the equations. 



Jo 



v|(j,/)Of (/) - ^-v\{y,y'W^\y')\ ,(6.87) 



4 (r^ r?) or (,) + ^ (r^ r?) or (,) 



+ dy' [vffy, y')Q>V(y') - \vliy, y'W''^iy') | = 



where. 



(vr(j,/)Of (/) - \v\{y,y'W^\y')^ ,(6.88) 



V7 = {v\ + v\ + v\)IA. (6.89) 



Taking k and y to zero in these equations gives equations for the behaviour of the correc- 
tions for small y. Since these equations become scale free at small y it means that the two 
corrections 0+ "\y) ~ as y becomes small. Hence, these 'corrections' radically change 
the form of the wavefunctions at small 3^ but provide only minor corrections at large y. 



6.2.4 Numerical Results 

We take the mass of each nucleon to the be isospin averaged mass, M = 938.9MeV. 
Eqn. (16.771) may be solved numerically using a method adapted from that used in the three 
Boson example. To set the interior phase defined by the LO three-body force, we apply a 
cut-off, Q., that allows us to 'feed' the desired phase into the equation. Below the cut-off, 
y < Q, we use the form given in eqns. ( 16.8116.821 ). for y > Q. the solution is to be determined 
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from the resulting homogeneous equation. In order to set the exterior boundary condition 
on the singlet solution we must also solve eqn. (16.771) to obtain a solution that has the small 
y asymptotic form given by eqn. (16.831) . 

Numerically, the solutions are again most easily found using linear algebra. The dis- 
cretisation of the y variable is done as before, however, since we are solving coupled equa- 
tions, the matrix equation to be solved is now of order 2n. 

At NLO, the numerical solution of eqns. (16.8716.881) is achieved similarly but is more 
complicated. Since these equations are still differential equations they require boundary 
conditions. The correct boundary conditions are those that remove all solutions to the ho- 
mogeneous equations, which are identical to the LO equations. That is, the solutions of the 
NLO equations should be the particular integral. To find these boundary conditions we must 
solve eqns. (16.8716.881) order-by-order in y. The results give, 

O^i^) ^ Jlly-^ sin(5o \nyy + j] + <), (6.90) 

0^i+) ^ Jl^_y-^ sin(5o Inyy + t] + crl), (6.91) 

OL^-' ^ JlZiy, - jt) sin(5o \nyy + rj + (t~), (6.92) 

^ ^;(r. - Jt) sin(5o Inyr + 77 + (r~\ (6.93) 

where the exact values of the constants ^ and cr are not important for this discussion. These 
boundary conditions are now fed into the equation by using the cut-off Q.. 

Fig. 16.61 shows the projection of the wavefunction, 0(y), of an elastic scattering wave- 
function at LO and NLO. The solid lines show the triplet wavefunction and dashed lines 
show the singlet wavefunction. The LO result is shown in red and the NLO result in black. 
Since there is no outgoing bound state in the singlet state, the singlet wavefunction vanishes 
as J ^ oo. The wave in the triplet channel corresponds to an outgoing/incoming deuteron 
and neutron. At small y the trig-log behaviour is apparent in the LO solution. The NLO 
solution also displays this behaviour with an extra factor oiy~^ . The LO and NLO solutions 
become distinct at around 2fm ~ rs,Pf 

To obtain physical variables, the LO three-body force must be fixed. This can be done 
by matching to the deuteron neutron scattering length, Und = 0.65 + 0.04fm [63]. With this 
single input the EFT becomes predictive for all other scattering observables including the 
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y [fm] 



Figure 6.6: Scattering solution to eqn. (16.751) with = Sl.OMeV {k = 13.6MeV) at LO 
(red) and NLO (black). The solid lines are the triplet wavefunction 0,(y), the dashed lines 
are the singlet wavefunction (^s(y) 
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triton binding energy. The three-body force at NLO will differ from the three-body force at 
LO by a term of order p^y. 

The relationship between the choice of three-body force, the neutron-deuteron scatter- 
ing length and the triton binding energy is shown in figs. l6.7l01 and l6!9l Fig. 16.71 shows the 
values of the neutron deuteron scattering length that correspond to different interior phases 
of the DWs and hence to different three-body forces. Fig l6. 81 shows how the binding en- 
ergy of the Triton (the shallowest bound state) depends upon our choice of three body force. 
Many of these shallowest bound states actually lie outside the validity of the EFT defined 
by the pion mass. 

Bringing the results of figs. 16.71 and 16.81 together gives a one variable relationship 
between the neutron-deuteron scattering length and the triton binding energy, fig. 16.91 Al- 
though this curve cannot be seen experimentally as the physical values correspond to a single 
point in this space, shown by a cross, it does have interesting implications for the study of 
three nucleon problems with potential models. The dots in fig. 16.91 correspond to different 
predictions for the triton binding energy and the neutron-deuteron scattering length obtained 
from different pair wise nucleon potentials with the same two body scattering lengths and 
effective ranges as inputs [52]. Instead of clustering randomly around the physical values, 
these form a curve through the space.. This relationship was first noted by Phillips [61] in 
looking at the predictions for three-body variables by nuclear pairwise potential models. In 
fig- 16.91 the NLO result (dashed) is a clear improvement upon the LO result (solid). The 
connection between the Phillips line and the three-body force in the three-body EFT has 
already been illustrated by Bedaque et al [24], our NLO curve differs slightly from theirs 
because of the additional term 16.761 

The interior phase of the DWs and the prediction for the Triton binding energy are 
shown in table 16. 1[ Fig. 16.101 shows the projection of the triton wavefunction at LO and 
NLO. The full LO Triton wavefunction obtained from eqn. (16.681 ) are shown in fig. 16.1 li on 
hyperpolar plots. The configuration of the three particles, parameterised in a, depends upon 
our choice of pair in defining a. In the triplet wavefunction, the 'chosen' pair have isospin 
and so correspond to a neutron and proton, in the singlet wavefunction the 'chosen' pair 
are the two neutrons. 
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Interior Phase, T| 



Figure 6.7: Neutron-Deuteron scattering lengths that result from choice of LO three-body 

force, corresponding to a different interior phase, rj, for the trig-log behaviour of the wave- 
functions. The solid and dashed lines show the LO and NLO results respectively 




Interior Phase, t| 



Figure 6.8: Triton binding energies that result from different choices of the LO three-body 
force. The solid and dashed lines show the LO and NLO results respectively 
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Figure 6.9: The Phillips line as predicted by the pionless EFT at LO (bold line) and NLO 
(dashed line). The dots correspond to the predictions for the triton binding energy and 
neutron-deuteron scattering length in different models with the same two-body inputs [52]. 
The cross is the experimental result. 





LO 


NLO 


Physical 


Interior Phase of DWs, 77, where 
sin(5o In fiy + rf) 
and yu = lOOMeV. 


2.442 


2.419 




Triton Binding Energy 


8.08MeV 


8.68MeV 


8.48MeV 



Table 6.1: The interior phase of the DWs and the Triton binding energy at LO and NLO 
found by matching to the neutron-deuteron scattering length. 
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The highest value of the wavefunction occurs as a proton and neutron are paired with 
a second neutron around Ifm away, which occurs at a = njb, R = Ifm in both the singlet 
and triplet wavefunctions and also at a = nil, R = Ifm in the triplet wavefunction. This 
configuration also allows the three nucleons to be spread over the greatest volume, i.e the 
highest density at large values of R corresponds to the values a = n/6 in both singlet and 
triplet states and also to a = in the triplet state. For i? < Ifm the wavefunction does not 
distinguish much between the different configurations of the three particles (different a's) 
except for or ^ where the wavefunction vanishes. This value, or = 0, corresponds to the 
situation where the third nucleon is directly between the other two. 

Fig. l6.12l shows the predictions for the elastic neutron-deuteron phaseshift after match- 
ing the three-body force to the neutron-deuteron scattering length at LO (solid line) and 
NLO (dashed line). The circles show the most recent experimental results (circa 1967) [64] 
and the triangles show the results obtained from a combination of VI 8 and Urbana IX two 
and three-body forces [62]. Given the age of the experimental results, in makes sense to 
compare the data to that obtained from sophisticated potential models, rather than the ex- 
perimental data. These models are fitted to many different variables in two- and three-body 
systems and can be expected to produce a reasonable curve [24, 33, 66]. 

The EFT curve fits the potential model curve well at NLO. The computing effort in 
solving the EFT equations for the wavefunctions and the physical variables is minimal com- 
pared with the effort required for the more sophisticated approaches. 

To go to higher orders than NLO requires considerably more effort. A new three- 
body force is required, isospin symmetry violation must be considered and the computa- 
tional effort in calculating the four second-order perturbations is considerable. However, 
the agreement achieved with very few parameters and at the relatively low order of (2/Ao)^ 
is pleasing. 

6.2.5 Summary 

In this chapter we have derived and solved equations for the three body DWs for three 
Bosons and three nucleons at LO and NLO in the KSW EFT. In each case a 'phase-fixing' 
three-body force is required. 
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Figure 6.10: The triton wavefunction's projection 0(j) as found by matching the three-body 
force to the neutron-deuteron scattering length at LO (red) and NLO (black). The solid lines 
are the triplet wavefunction 0,(j), the dashed lines are the singlet wavefunction <!>siy) 
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Triplet state 'i'^R, a) 




12 3 4 



R [fm] 

Figure 6.11: Contour plots showing the triton wavefunction at LO. The highest particle 
density occurs at a = n/6, R ~ Ifm in both the triplet and singlet wavefunctions. This 
corresponds to a proton and a neutron on top of each other with the third neutron Ifm 
away. In the triplet wavefunction the a = n/2 'tail' is longer and peaks higher as this also 
corresponds to a proton neutron pair with a distinct neutron. 
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Figure 6.12: Physical phaseshifts for neutron-deuteron scattering (^o = V^/Sy^ - k^) with 
the LO three-body force fitted to the neutron-deuteron scattering length. The solid and 
dashed lines are the LO and NLO result respectively. The circles are the most recent exper- 
imental results from 1967 [64] and the triangles are the results obtained using the Argonne 
V18+Urbana IX two and three-body forces [62]. 
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The implications of the LO three-body force in the case of three nucleons is embodied 
in the Phillip's line, which shows the one parameter relation between the neutron-deuteron 
scattering length and the triton binding energy as predicted by pairwise nuclear potentials. 
The NLO results show convergence upon the potential model results in both the neutron- 
deuteron phaseshift and the Phillip's line. 

The equations, once the three-body force is determined by fitting to the neutron-deuteron 
scattering length, yield the triton wavefunction and all DWs describing elastic scattering of 
a neutron and a deuteron below breakup. The triton wavefunction result can be used in 
applications, such as electron triton scattering. 



Chapter 7 



Conclusions 



In this thesis we have looked at the use of the renoraialisation group in the development of 
effective field theories. Our three main results are the development of a tool, the DWRG 
and general methods of solving the DWRG equation; the use of the 3BDWRG equation in 
deriving the limit-cycle solution for the three-body force in the KSW EFT for short-range 
forces and the corresponding power-counting; the development and solution of equations 
that can be used to practically evaluate three body DWs in the KSW EFT for Bosons and 
the physically and mathematically interesting system of three nucleons in the ^5 1 channel. 

The DWRG is tool that enables us to resum some diagrams to all orders. The general 
method of solving the DWRG equation relies upon the basic loop integral, perturbations 
about which give an alternate power-counting scheme to naive dimensional analysis that 
corresponds to a trivial fixed point in the DWRG. What is pleasing about the results we have 
obtained, using the DWRG, is that established results such as the DWERE for well-behaved 
potentials share a common method with new results such as the limit-cycle solutions for 
short range forces with an attractive inverse square potential. 

The solution of the DWRG for the attractive inverse square problem required a little 
care. The system is not properly defined without the input of a scale, po, that defines a 
self -adjoint extension. The solutions of the DWRG are in the form of limit-cycles. The 
power-counting associated with the limit-cycle solutions have a LO marginal term, which 
corresponds to the degree of freedom offered by the choice of po. The solution of the 
DWRG answers many issues raised in refs. [30, 32] concerning the connection between 
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the short range effective force and the need for a self-adjoint extension, the possibility of 
multi-valued solutions of the RG, and the lack of a ground state. 

Our study of the KSW EFT three-body problem was motivated by the recent realisation 
that the three-body force is extremely important in some systems [22, 24]. Two systems 
that are particularly interesting are the three s-wave Bosons and three nucleons in the i 
channel. These systems have no net repulsion due to the pair-wise forces and result in 
singularities similar to the attractive inverse square potential. Because of this similarity, 
the 3BDWRG solutions exhibit limit-cycle behaviour and a three-body force that occurs at 
order (Q/Aq)'^. Our results confirm the power-counting proposed by Bedaque et al 
The power of our approach is that we are able to give a simple algebraic statement of that 
counting. 

Our results also illustrate the limit cycle behaviour of the three-body force that has been 
suggested based on analysis of the attractive inverse square potential [30, 32] and the STM 
equation [22,24]. 

The generalisation of the 3BDWRG to four or more bodies is straightforward, our 
method of solution involving the basic loop integral is easily generalised. The scale-free 
systems studied in this thesis and the novel power-counting schemes that correspond to 
them are likely to be important in N-body forces in the KSW EFT due to the scale-free 
nature of these systems at hyperradii far smaller than the two-body scattering length. 

The final chapter offers some physical results for the KSW EFT. The convergence 
of the results using a limit-cycle three-body force and their agreement with the results of 
sophisticated potential methods has already been shown, [24, 33]. To complement these 
results and to offer new possibilities of testing the EFT we have derived equations that 
allow calculation of the DWs for the three-body KSW EFT below threshold. Our equations 
are the coordinate space equivalent of the STM equations [56]. Taking advantage of what 
we had shown earlier, the LO three-body force can be chosen by simply fixing the phase of 
the DWs at small hyperradii. 

Our results for three nucleons in the i channel agree well with those of Bedaque et 
al (up to a difference in definition of NLO corrections). Our method allows us to calculate 
particle density of the triton wavefunction that may be tested by looking at electron-triton 
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scattering data. 

In summary, the DWRG has shown itself to be a versatile tool and has been used to 
shed light on several issues currently under debate. This thesis has produced several results 
that provide insight into current issues in the EFT community and open new doors for testing 
the KSW EFT for short-range forces. 



Appendix A 
Analyticity of J{p, k) 



We want to study the analyticity of the integral, 

q{e « - 1), 



/ = r dqCiklq) = InK f (A.l) 

-'o Jo a(eT - n. 



which we can write as, 

i = ink\ — — - = inAk - -nr-rv (^.2) 

Jink Ae^ - 1) J Ae"" - 1) 

where A is some constant. Since we are investigating the behaviour of the integral near 
/c ^ we can expand the integrand in the final term in a series in x and assume k to be 
within the radius of convergence to get, 

/Ink / 1 1 1 \ Tr'ii 
dx\ — + — + ... =1+ nk(lA + In)^ + ln2;r) + . . . . (A.3) 
\x^ 2x 12 I 6 
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Appendix B 

DWRG Analysis for the Attractive 
Coulomb Potential 



In this appendix we shall look quickly at the modifications to the DWRG analysis for an 
attractive Coulomb potential, i.e. when k is negative. 

Clearly, the power-counting around each fixed point is unchanged as these result from 
the form of the DWRG equation rather than details. The trivial fixed point is also unchanged, 
the expansion around it yields the distorted wave Bom expansion, eqn. (13.311) . 

The two differences that result from having an attractive Coulomb potential is the ex- 
istence of bound states that occur in the truncated Green's function and the definition of the 
non-trivial fixed point. In the repulsive case the starting point for the non-trivial fixed point 
solution was the basic loop integral. 



However, in the repulsive case this integral no longer converges because of the essential 



where c is the contour shown in fig. IB. 11 This integral is convergent because it approaches 
the essential singularity from the other direction. 




(B.l) 



singularity at § = 0. To remedy this, we now define J as 




(B.2) 
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Figure B. 1 : The contour of integration in the complex ^-plane for the non-trivial fixed point 
in the attractive Coulomb case. The bound states appear as poles at q = ik/n. 

The integral J can now be split up as before to isolate the singularities: 

J(p,k) = -Re r dqCik/q) - f dq^^^ + M(p,k), (B.3) 
Jc Jl P ~ <3 

M(p,k) = p'Rc f dq^^, (B.4) 

Jc' p - q 

where the contour c' extends the contour c along the real axis to infinity. As in the repulsive 
case, the second term is analytic, the term'M is not analytic but satisfies the homogeneous 
DWRG equation and the first term contains logarithmic dependence, which is still given 
by the analysis in the previous appendix. Hence we can write the non-trivial fixed point 
solution precisely as before, 

Vf = (J\pk) - M(p,k) - l(k,A)y\ (B.5) 

where X is defined in eqn. (13.391) . Now we have found a non-trivial fixed point solution the 
comments on the RG flow are the same as in the repulsive case. 

We will find that when writing down the distorted wave eff"ective range expansion there 
is the additional bound state terms in the truncated Green's function. Eqn. (13.471 ) should 
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now read, 
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where i^„(r) is the bound state with binding energy E„ = \k\ /n and is given by, 



>^n(r) = J-^2re-''l"L\_,{2Krln) ^ Iri-] as r ^ 0, (B.7) 



where L" is the associated Laguerre Polynomial. The term in braces in eqn. ( IB. 61 ) may be 
written as. 
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and so by Cauchy's theorem is equal to MK^J{p, k). The analysis follows as before, result- 
ing in the DWERE, 
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where. 
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Appendix C 

The Faddeev Equations for Three 
Identical Bosons 

For a system with three particles the Lippmann-Schwinger equation is inadequate. Although 
the three-body T-matrix satisfies the Lippmann-Schwinger equation, it is not possible to use 
this equation to find it because its kernel is not compact. This may be understood in several 
ways. Most simply in coordinate space it is a result of a non- vanishing potential as the mean 
separation of the three particles is taken to infinity but two of the particles remain close. This 
in turn is related to the possibility of incoming and outgoing boundary conditions involving 
bound states of particles and also to the possibility of so-called disconnected diagrams, in 
which one of the three-particles fails to interact. 

In short, a wave-function with an incoming two-body bound state and a single par- 
ticle satisfies the homogeneous version of the Lippmann-Schwinger equation resulting in 
non-unique solutions. The resolution to this quandary is well-known. We must ensure that 
disconnected diagrams are not allowed. One favoured method is to use the Faddeev equa- 
tions. 

Consider the case of three identical bosons interacting via a pairwise interaction V2b- 
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The wavefunction must be symmetric with respect to transposition of any two particle in- 
dices, hence the wavefunction, |*F^) with incoming boundary conditions of three free parti- 
cles (rather than a bound state of any two) can be written as 

\^^) = \ii^o)+{i+p)\r), (c.i) 

where \(ffo) is the 'in' -state and solves the free Schrodinger equation and P is a permutation 
operator, which permutes particle indices and has the matrix elements defined in eqn. (I5.3I) . 
The Schrodinger equation is written as, 

(Ho + a+P)V2B-EW^) = 0, (C.2) 
(1 + P)iHo + V2B - E)\ilj^) = -(1 + P)V2b(\>I^q) + P\il^^), (C.3) 

where to obtain the second line, we have used the fact that Hq commutes with the permuta- 
tion operator, P. Cancelling the (1 -i- P) on the left hand side and using G^^ = (E - V2B - 
Hq + iey^ we obtain, 

\r) = G\s{E)V2BW^) + G^^s{E)V^^P\r), (C.4) 
= Gl{E)t\E)\il,o) + Gl{E)t\E)P\ilj^), (C.5) 

where is the T-matrix given by the equation, f^iE) = V2B + VibG^b^^^^'^^- ^'l'^- 
is the Faddeev equation for the wavefunction component It is connected because of 
the P-operator which ensures that each interaction, given by the two-body T-matrix, is not 
between the same two particles as the preceding one. 

The corresponding equations for the full three-body Green's function and T-matrix 
now follow from similar decompositions. If we write the full T-matrix, T in terms of 
components, 

T^{E) = {l+P)T^{E), (C.6) 

then from the equation defining the full wavefunction in terms of the MoUer wave operator 
l+Gl{E)T\E), 

= (1 + Gl{E)T\E))\il,o) = 1^0) + (1 + P)G^o(E)T^\ME)), (C.7) 

and from the definition of in terms of its components, eqn. (15.11 ). we obtain the relation, 
= Gq{E)T'^{E)\iIjo), which when substituted into the Faddeev eqn. (15.41 ) yields the 
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Faddeev equation for the T-matrix component, 

T*{E) = f{E) + f{E)Gl{E)PT*{E). (C.8) 

We may find a similar relation for the full Green's function, g*{E) = G^(E) + (1 + P)G*{E), 
by using the equation g^{E) = G^(E) + G^(E)T*(E)G^(E). 



Appendix D 



Momentum Dependent Perturbations in 
the DWRG for Three-Bodies 



Momentum dependent perturbations about the limit cycle solution are found using the ansatz 



V3b(p, 7, K k'. A) = Vf^ip, r, A) + CA>(p, y. A), 



3B 



(f>ip, r. A) = k^"cPiip, y, A) + MP, % A), 



(D.l) 
(D.2) 



where cpi and (()2 only depend upon A logarithmically. Substituting this into the DWRG 
equation (I5.68I) . neglecting for now the discontinuities, then linearising gives a partial dif- 
ferential equation for 0i and 02, 
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where C{y, InA/p^) is given by, 
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From the first line we obtain the RG eigenvalue v = 2n and the solution 0i = . So that a 
with the further substitution. 



4.2 = [4',V<S-I>'"]V'S 



(D.5) 



t cr\ 
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we obtain a PDE for 03, 

oA op oy 1 - 

This equation should of course have discontinuities on the RHS. The solution of this PDE, 
including those discontinuities, is straight forward and can be reached in much the same 
way as the limit-cycle solution. The momentum perturbations occur at orders 2, 4, 6, . . . at 
the same orders of the energy perturbations. 
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Fourier Transforms 



We wish to perform the Fourier transform, 
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Doing the integral over all angles in /i gives, 

^1 = ^ Jfcfcsin(;t|y - y'\)e-^'^\ (E.3) 

27r2x|y-y'| Jo 

where we have put k = -2p/ V3. Using the result 

f dx xe-f^ sin Zjx = -f^^2(r ^I^^Th^) (E.4) 
Jo yS-^ + Z?^ 

we get, 

where (3o is given by eqn. (16.311) . The second integral, I2, is more involved, making the 
change of variable k" = k' + k/2 we obtain. 
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Performing the k" integral gives/2 = 1+ + I- where, 

/ - ^ r /k.(y+ly') -^|y'±x|VFT? 

--\r±x\J {inf • ^""-'^ 

Then the angular k integral can be done to give, 

h = ■ dkk^m{k\y + -y'\)e-^^^'^'-^'^' (E.8) 

- 2;r2|y+iy,||y, + x| Jo 2^'^ 

(E.9) 



2;r2|y+iy'||y'±x| Jo 



An'- Ql ' 

where we have again used result (IE.4I) . The result quoted in the text follows 



The second integral to be done is 

I6n r d^k r d\' 



Vo(y,y') = - 



\6n r d% r d'k' 
3M J (2^J (2^ 

r d\ r d^k' 



g!(k.y-k'.y') 



L2(;7,k,k',0)e 



/(k.y-k'.y') 



+ - 



16;r 



V3 



^ F + + k.k' + |i^2 
J (2;r)3 J (2;r)3''^ 



F + k'^ + k.k' + |/c2 



(E.IO) 
(E.ll) 



/3 but for an overall numerical factor is the same of the I2 with x = 0. Hence we have, 



h = 



(E.12) 



' V3;r2 ^2 ■ 

I4 can be done in much the same way as the I2 integral (with x = 0) all the way up to the 
final k integral. Hence we get. 



h = 



z dkk^^i^TK^sm(k\y+-y'\)e-'Ty''^^. (E. 

37r2|y + iy'ly Jo 2 



V37r2|y + iy'|/ 
Differentiating the result (IE.4I) on both sides w.r.t. /3 gives 



13) 



dxx^y^^e-P<^-^mbx = J^^^^fn -^^(7 ^W^') ' ^^2^2(7 ^J0^'), 
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which can be used to obtain, 
The quoted result follows. 



(E.15) 



The final integral is 

8;r r d% f d%' [ [3~~ 

_ 167ry2 r d^k r d^k 
~^~J (2^J On) 



g((k.y-k'.y') 



L2(;?,k,k',0)e'('^-^-'^'^') 



)3 + k'^ + k.k' + |/c2 
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\6ny 
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VP + /c2 



r d^k r d^k' 

r d\ r d^k' , , 



g/(k.y-k'.y') 



F + k'^ + k.k' + ^2 



g!(k.y-k'.y') 



F + + k.k' + |/f2 



rl2^3 + /4| + /5- 



(E.16) 



The results for 73 and I4 have already been established. I5, like I4, can be done in much the 
same way as the I2 integral (with x = 0) all the way up to the final k integral. Hence we get, 



1 1 

—— dkkie + l^)smik\y + -y'\)e- 

+ h'\y' Jo 2 



2;r21y+iy'|y 

Differentiating the result (IE. 141) on both sides w.r.t. /3 gives 

33 A, ,4 



(E.17) 



J^^ dx x(y^ + x^)e-^ sin bx = ^||^^4(r A^^) - Jf^^y^ K^iy ,^^), 

(E.18) 

which can be used to get, 

^ 3 V3/f4y'2 K4(kR) 3 V3/f3 



' 16;r2 ^4 4^2 ^3 
All the results brought together give the result in the text. 
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Appendix F 



Zero Energy Wavefunction for the Three 
Body EFT 



We wish to evaluate the integral 



QlQl 



where. 



/3 /3 3 

"R^ = -x^ + y'^ + 2yy', = ^1 -x^ + j2 + y2 + yy + _-^f^ 



In order to use residue calculus to evaluate the integral we write it as 



V3 



h = - 



4 



loo -^-^ nlnl ' 
QlQz 



d/Z^-Mog 



(F.l) 

(F.2) 

(F.3) 
(F.4) 

(F.5) 



V3;r(l - e'^'O 

The poles of the integrand of h are given by the zeroes of :/?± which may be written in terms 
of the hyperpolar coordinates as. 



y = +i—Re^"'. 
2 



(F6) 



Evaluating the integral using the two residues in the upper half -plane (recalling that < or < 
7r/2) we get, 

( V3/?V /7rs\ 

(F7) 
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Figure Fl: The branch structure of the integrand in I2, the labels show the argument of the 
branch point. On the left the case of < a < njb and on the right njb < a < n/2. 



I2 can be evaluated by closing the contour in the upper half plane then collapsing it onto the 
branch cuts there. The branch points are given by the zeroes of which are simply, 

y' = and y' = ^Re^'("^^-'(^K (F8) 

How the branch points up are joined up depends on the value of a. We do not want any of 
the branch cuts to cross the real axis. The branch structures are shown in Fig. IF. II 
Simply evaluating the integrals around the branch cuts gives, 



I7 = 



I 2 ) ,V3 i+2cos(f) ' wnenQf< ^, ^^^^ 

(j!MV whena>^ 
t\ 2 ) ,V3i+2cos(f) wnena > ^. 



Finally using the equation for s (16. 341) . and combining the two results we obtain the quoted 
result. 
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